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Abstract 

LAMINAR  AND  TURBULENT  BOUNDARY  LAYERS 
ADJACENT  TO  MELTING  VERTICAL  ICE 
WALLS  IN  SALT  WATER 

by  Edward  George  Josberger 

A laboratory  study  of  melting  vertical  ice  sheets  immersed  in 
sodium  chloride  solutions  of  oceanic  temperatures  and  salinities  found 
that  the  convective  motions  generated  by  the  melting  consist  of  two 
regions  confined  to  a boundary  layer  adjacent  to  the  ice.  The  first 
region,  at  the  bottom  of  the  ice,  is  laminar  and  bidirectional,  with  a 
2 to  3 nui  thick  upward  flowing  layer  next  to  the  ice  inside  of  a 10  mm 
thick  outer  downward  flowing  layer.  This  bidirectionality  results  from 
the  large  difference  between  the  thermal  and  saline  diffusivities  which 
contain  dilute  water  near  the  ice  and  allows  the  cooling  to  diffuse 
further  from  the  ice.  Near  the  ice,  dilute  cold  water  rises;  while  away 
from  the  ice,  cold  saline  water  sinks.  Second,  further  up  the  ice,  the 
inner  layer  goes  through  a region  of  transition  and  becomes  fully 
turbulent  when  the  saline  Grashof  number  reaches  2 x 108.  The  turbulent 
flow  is  upward,  unidirectional  and  the  dominant  flow  in  the  oceans. 

Ice-water  interface  temperature  measurements  over  the  ice  surface 
show  a vertically  varying  temperature  In  the  laminar  region  and  a uniform 
temperature  in  the  turbulent  region.  The  measured  melt  rates  were 
highest  in  the  transition  region,  lowest  in  the  laminar  region,  and 
had  a slow  vertical  variation  in  the  turbulent  region.  For  water 


temperatures  greater  than  2 5°C,  the  turbulent  flow  reversed  direction 
and  a laminar  bidirectional  region  formed  at  the  top  of  the  ice. 

In  a heoretical  study  of  the  turbulent  flow,  an  eddy  diffusivity 
models  the  turbulent  transport  of  momentum,  salt  and  heat  in  the 
conservation  equations  and  the  laboratory  measurements  determine  the 
model  parameters.  The  turbulent  model  consists  of  three  regions:  a 
laminar  inner  region,  an  intermediate  region  where  the  eddy  diffusivity 
increases  linearly  with  distance  from  the  ice  and  an  outer  region  of 
constant  eddy  diffusivity.  A similarity  transformation  reduces  the 
governing  partial  differential  equations  to  a system  of  ordinary 
differential  equations  which  are  then  numerically  integrated. 

The  numerical  results  show  the  following:  The  eddy  diffusivity 
depends  parabol leal ly  on  the  density  difference  across  the  boundary 
layer.  The  major  effect  of  the  turbulence  is  to  increase  the  transport 
of  dilute  water  away  from  the  interface.  The  dilution  effect  overwhelms 
the  thermal  effect  on  density  to  produce  a unidirectional  upward  flow. 
Finally,  the  injection  of  melt  water  into  the  boundary  layer  reduces 
the  salt  flux  to  the  ice  which  lowers  the  interface  salinity  as  the 
water  temperature  increases.  Use  of  the  results  to  form  a zero  net 
buoyancy  condition  on  the  buoyancy  integrated  across  the  boundary  layer 
predicts  the  reversal  of  the  turbulent  flow  which  accurately  fits  all  of 
the  known  data. 

Salinity  and  temperature  measurements  with  depth  in  the  vicinity  of 
an  iceberg  in  the  North  Atlantic  Ocean  show  variations  in  these  quanti- 
ties near  the  iceberg  that  may  be  the  result  of  the  upward  flowing 
turbulent  boundary  layer. 
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INTRODUCTION 


Ice,  in  the  form  of  glacier  and  sea  ice,  is  the  most  striking 
feature  of  the  polar  ocean'-.  The  discovery  of  off-shore  oil  deposits  in 
the  Arctic  Ocean,  the  possible  establishment  of  a krill  fishery  in  the 
Southern  Ocean,  and  icebergs  as  a possible  source  of  fresh  water  makes 
understanding  the  complex  ice-sea-water  interaction  phenomena  imperative. 
While  the  formation,  distribution  and  movement  of  pack  ice  in  the  ocean 
has  received  considerable  attention  in  AIDJEX,  the  melt  processes  have 
been  neglected.  In  particular,  the  melting  of  a vertical  ice  wall  in 
seawater  has  not  been  examined.  In  this  case,  the  melting  cools  and 
dilutes  the  seawater  adjacent  to  the  ice  to  produce  convection.  This 
convection  results  from  the  complex  diffusion  of  salt,  heat  and  momentum 
where  the  thermal  and  saline  effects  on  density  oppose  each  other. 

In  the  ocean  the  vertical  length  scale  for  sea  and  glacier  ice 
varies  from  250  m to  0.5  m.  Icebergs  have  the  largest  vertical  length 
scale  of  all  oceanic  ice,  as  large  as  250  m.  An  iceberg  forms  when  a 
tidewater  glacier  or  an  iceshelf  calves  or  breaks  and  sets  a large  piece 
of  glacier  ice  drifting  on  the  ocean.  Paterson  (1975)  defines  glacier 
ice  as  snow  that  has  metamorphosed  to  the  point  where  the  trapped  air 
pockets  are  sealed  off  from  one  another;  the  trapped  air  volume  varies 
from  0 to  7%.  The  vast  ice  shelves  of  the  Antarctic  Continent  produce 
tabular  icebergs  that  have  depths  up  to  250  m and  a typical  horizontal 
scale  of  1 km,  although  icebergs  as  large  as  50  km  on  a side  have  been 
reported  (McClain,  1978).  The  Amery,  Ronne-Filchner  and  Ross  iceshelves 
produce  the  majority  of  the  icebergs  found  in  the  Southern  Ocean  and 
Hult  and  Ostrander  (1973)  estimate  the  input  of  ice  to  be  6 x 1015  kg 
yr"1.  The  effect  of  this  ice  input  on  the  oceanography  of  the  Southern 
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Ocean  may  be  significant  and  only  a thorough  understanding  of  the  melt- 
ing process  will  determine  the  ice  melting  effects. 

Glaciers  flowing  from  the  Greenland  ice  cap  in  the  North  Atlantic 
produce  icebergs  that  are  smaller  and  fewer  in  number  than  the  icebergs 
in  the  Southern  Ocean.  Plate  XIII  shows  a typical  North  Atlantic  ice- 
berg; it  measures  approximately  120  m on  the  water  line,  rises  30  m out 
of  the  water,  and  has  a depth  of  100  m.  North  Atlantic  icebergs  are 
usually  unstable  with  rolling  occurring  frequently.  The  number  of  ice- 
bergs varies  from  year  to  year  but  in  a heavy  ice  year  over  a thousand 
icebergs  find  their  way  south  of  the  Grand  Banks  (Neuman  and  Pierson, 
1966)  where  they  constitute  a hazard  to  navigation.  In  1912  the  Titanic 
struck  an  iceberg  near  the  "Tail  of  the  Bank"  and  sank  with  the  loss  of 
1,502  lives  (Lord,  1976).  Today,  even  with  radar,  small  bergs  or 
growlers  may  not  be  detected  in  heavy  seas  or  pack  ice  and  they  still 
represent  a hazard  (Ruffman,  1976).  In  1977  a Canadian  ferry  bound  for 
Labrador  presumably  struck  a growler  while  operating  in  pack  ice  and 
sank  without  loss  of  life. 

Icebergs  are  not  only  a hazard  to  navigation,  but  they  may  also 
increase  biological  productivity  around  the  iceberg  through  upwelling 
driven  by  dilution  from  the  melting  ice.  In  the  Southern  Ocean  icebergs 
extend  below  the  main  pycnocline  at  50  m and  the  convection  next  to  the 
iceberg  may  transport  nutrients  from  the  deep  water  into  the  photic 
zone.  However,  controversy  exists  over  the  effect  of  the  melting  on  the 
surrounding  seawater.  Neshyba  (1977)  states  that  upwelling  is  important 
in  the  primary  production  of  the  Southern  Ocean,  while  Huppert  and 
Turner  (1978)  feel  that  stratification  in  the  Southern  Ocean  limits  the 
vertical  convection  and  no  significant  upwelling  occurs.  Neshyba  (1977) 
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estimates  that  the  convection  associated  with  the  sidewall  melting 
supplies  significant  amounts  of  nutrients  from  below  the  pycnocline  at 
50  m to  the  surface  layers  for  use  in  photosynthesis.  Neshyba  based  his 
calculations  on  turbulent  heat  transfer  coefficients  derived  from 
experiments  where  heat,  not  salt,  created  the  buoyant  driving,  a proce- 
dure that  may  not  be  correct.  In  a laboratory  experiment  Huppert  and 
Turner  (1978)  investigated  the  effect  of  stratification  on  the  convec- 
tion adjacent  to  a vertical  ice  wall  and  found  a cellular  convection 
pattern.  In  their  experiments,  with  0.25  m long  ice  slabs  submerged  in 
linearly  stratified  sodium  chloride  solutions,  the  cold  water  near  the 
ice  sank  a few  10' s of  mm  until  the  melt  water  reached  its  own  density 
level  where  it  flowed  horizontally  out  into  the  interior  of  the  tank. 

As  a result  of  this  cellular  convective  pattern,  horizontal  ridges 
formed  on  the  ice  surface  with  a vertical  length  scale  equal  to  the 
convective  cell  length  scale.  They  assert  that  the  vertical  step 
structure  in  temperature  and  salinity  found  in  the  Weddell  Sea  may  be 
due  to  melting  icebergs.  They  performed  their  experiments  at  20°C  with 
salinity  gradients  up  to  100  o/00  m-1,  conditions  which  do  not  represent 
oceanic  values.  Also,  for  uniform  conditions  I found  a reversal  in  the 
flow  at  20°C  and  oceanic  salinities.  These  factors  make  direct  applica- 
tion of  Turner  and  Huppert's  work  to  the  ocean  inappropriate;  the 
reversal  at  20°C  is  discussed  in  the  results  section. 

To  determine  the  effect  of  a melting  ice  wall  on  the  surrounding 
water  properties  in  the  ocean,  I performed  a field  program  adjacent  to 
an  iceberg  in  the  North  Atlantic  Ocean  in  cooperation  with  the 
International  Ice  Patrol,  United  States  Coast  Guard.  The  results,  given 
in  Appendix  III,  show  cooling  and  dilution  of  up  to  3°C  and  0.5  °/00  in 
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the  upper  30  m of  the  water  column  within  200  m of  the  iceberg  which 
supports  the  concept  of  an  upward  flowing  boundary  layer.  Plate  XIII 
shows  the  iceberg  and  a prominant  waterline  notch,  denoted  by  arrow  A, 
which  shows  that  the  iceberg  has  recently  rolled.  The  upper  horizontal 
edge  of  the  iceberg  was  once  a vertical  wall  and  this  polished  wall 
shows  no  sign  of  horizontal  ridging  on  the  scale  of  ten  meters  as  would 
be  expected  if  stratification  were  important. 

Since  icebergs  are  pure  water,  they  have  long  been  considered  as 
a source  of  water  for  drinking  and  irrigation.  For  example,  H.R.H. 
Prince  Mohamad  Al-Faisal,  a Saudi  Arabian  Prince,  has  proposed  the  tow- 
ing of  tabular  icebergs  from  Antarctica  to  the  Gulf  of  Arabia  to  provide 
fresh  water  for  his  arid  country  (Husseiny,  1978).  For  this  scheme  to 
be  successful,  the  melting  processes  must  be  thoroughly  understood  over 
a wide  range  of  temperatures,  from  near  the  freezing  point  in  the 
Southern  Ocean  to  greater  than  20°C  in  the  Indian  Ocean.  The  convective 
mixing  processes  must  also  be  understood  if  the  collection  of  the  sub- 
surface melt  water  is  used  to  provide  the  fresh  water. 

Not  all  ice  in  the  oceans  is  in  the  form  of  icebergs;  sea  ice 
covers  most  of  the  Arctic  Ocean  and  the  Southern  Ocean  during  their 
respective  polar  winters.  Pressure  ridge  keels  made  of  sea  ice  and  sea 
ice  adjacent  to  open  water  have  vertical  scales  of  0.5  to  10  meters. 

Cox  and  Weeks  (1974)  show  that  sea  ice  contains  salt;  first-year  ice 
may  have  a salinity  as  high  as  20  °/00  while  multiyear  ice  usually  has 
a salinity  of  2 °/00.  The  low  multiyear  Ice  salinity  allows  the  results 
of  this  study  to  be  applied  to  most  sea  ice  melt.ng  phenomena. 


Pressure  ridges  form  when  two  sheets  of  sea  ice  converge,  usually 
under  atmospheric  forcing,  and  undergo  mechanical  failure.  Blocks  of 
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sea  ice  pile  up  and  form  a ridge;  the  sail,  which  is  the  portion  of  the 
ridge  above  the  water,  is  usually  a few  meters  high,  but  it  is  supported 
by  a keel  tens  of  meters  below  the  water  (Parmerter  and  Coon,  1972). 
Pressure  ridge  keels  menace  submarine  navigation  and  are  important  in 
modeling  the  motions  of  the  Beaufort  ice  pack  (AIDJEX  Bulletin  No.  23). 
The  random  orientation  of  the  ice  blocks  in  the  keels  gives  the  melting 
of  tilted  ice  walls  geophysical  importance. 

Ice  with  vertical  length  scales  of  0.5  m to  1.0  m occurs  when  a 
lead  or  a polyna  opens  in  the  pack  ice.  Solar  radiation  warms  the  sea- 
water and  a fraction  of  this  heat  melts  the  ice  at  the  lead  edges  form- 
ing convective  boundary  layers.  In  a similar  situation,  the  spring 
runoff  from  rivers  flowing  into  the  Arctic  Ocean  produces  warm  low 
salinity  water  that  rapidly  melts  the  sea  ice.  In  this  case,  the 
salinity  can  be  below  25  °/00  and  the  temperature  between  -1.2°C  and  4°C. 
As  a result  of  the  nonlinear  equation  of  state  for  seawater,  in  the  lead 
case,  the  temperature  and  salinity  anomalies  due  to  melting  produce 
opposing  buoyancy  forces;  in  the  river  case,  the  anomalies  produce 
reinforcing  buoyancy  forces. 

Vertical  natural  convection  adjacent  to  ice  walls  occurs  in  the 
polar  oceans  over  a wide  range  of  length  scales  driven  by  a wide  range 
of  ambient  temperatures  and  salinities.  To  understand  fully  the  role 
of  ice  in  the  physical  oceanography  of  the  polar  oceans,  the  fluid 
dynamics  of  this  convection  must  be  studied. 

In  the  following  two  subsections  I discuss  first  the  equation  of 
state  and  the  freezing  point  dependence  of  seawater;  second,  the  dimen- 
sional analysis  applicable  to  the  flow  next  to  a melting  vertical  ice 
wall . 
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Freezing  Point  and  Density  of  Seawater 

The  convection  generated  by  a melting  vertical  ice  wall  is 
a complex  problem  in  fluid  dynamics  because  it  involves 
a phase  change,  a fluid  density  that  depends  on  both  temperature  and 
salinity  in  a nonlinear  fashion  and  the  diffusion  of  momentum,  heat  and 
salt.  The  difference  in  the  molecular  coefficients  of  diffusion  further 
complicates  the  problem. 

When  ice  melts  in  seawater,  it  cools  and  dilutes  the  adjacent 
water.  The  ice-water  interface  temperature  and  salinity  must  be  on  the 
freezing  curve  given  by  Fujino  et  aK  (1974)  as 


for 


Tf  ■ -0.036  - 0.0499S  - 0.000112S2  (1) 

17.7  °/00  < S < 35  o/00 


where  Tfp  is  the  freezing  point  and  S is  the  salinity.  For  S < 17.7  °/00 
I computed  T^p  using  the  results  given  in  Kaufman  (1960).  The  interface 
conditions  as  derived  in  the  boundary  conditions  chapter,  also  depend  on 
the  melting  rate  so  the  wall  temperature  Tw  and  the  wall  salinity  Sw 
are  not  known  a priori. 

The  ambient  far-field  temperature  T^  and  salinity  determine 
the  far-field  freezing  point  and  the  temperature  elevation  above  the 
freezing  point,  which  I call  the  thermal  driving,  Td> 


T . = T - Tfn(S  ) . 

d « fp  °° 


(2) 
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The  thermal  driving  is  a convenient  temperature  parameter  to  use  in  ice 
melting  problems  because  it  is  determined  solely  by  the  far-field 
conditions  and  is  independent  of  the  unknown  interface  conditions. 

When  fresh  ice  melts  in  seawater,  convection  results  from  the 
cooling  and  diluting  of  the  surrounding  water.  If  the  interior  ice 
temperature  equals  T , then  all  of  the  heat  transferred  to  the  ice 

W 

produces  a phase  change,  which  requires  a latent  heat  L of 

L = 3.3  x 105  J kg*1  . (3) 

To  compute  the  buoyancy  that  the  melting  generates  requires  an 
accurate  formulation  of  the  equation  of  state  for  seawater.  For 
compactness  and  ease  of  manipulation,  I use  the  formulation  of  the 
equation  of  state  given  by  Gebhart  and  Mollendorf  (1977)  at  1 atmosphere 
of  pressure.  Appendix  I gives  both  the  formulation  and  the  values  of 
the  coefficients  used  in  the  formulation.  Figure  1 shows  the  nonlinear 
behavior  of  the  density  of  salt  water  as  a function  of  S and  T.  The 
existence  of  a density  maximum  exemplifies  the  nonlinear  dependence. 

The  line  of  maximum  density  lies  between  S = 0 and  T = 4°C,  and 
S = 25  °/00  and  T = -1.4°C  where  it  intersects  the  freezing  point  line; 
for  temperatures  > -1.4°C  and  salinities  > 25  °/00  the  maximum  density 
lies  beneath  the  freezing  curve  in  the  frozen  zone.  The  figure  shows 
that  dilution  of  salt  water  always  produces  less  dense  water  while 
cooling  produces  either  less  dense  or  more  dense  water,  depending  on 
the  location  of  T and  S with  respect  to  the  maximum  density  line. 


TEMPERATURE  (°C) 


SALINITY  (%o) 

Figure  1.  A T-S  diagram  showing  the  location  of  the  freezing  curve, 

TfD,  and  maximum  density  curve,  Tpmax.  Far-field  conditions 
in  region  I produce  reinforcing  upward  buoyancy  forces;  far- 
field  conditions  in  region  II  produce  opposing  buoyancy 
forces. 


9 


At  the  maximum  density  line,  the  coefficient  of  thermal  expansion  a, 
where 


- 1 3p 
' 0 3T 


(4) 


changes  sign.  Therefore,  when  T and  S lie  in  region  II,  cooling  and 
dilution  produce  opposing  buoyancy  forces  because  a is  negative.  When 
T and  S lie  in  region  I,  cooling  and  dilution  produce  reinforcing 
upward  buoyancy  forces  because  a is  positive.  In  fresh  water  the 
density  maximum  at  4°C  gives  opposing  buoyancy  forces  when  lies 
between  4°C  and  8°C. 

For  most  oceanic  cases  the  buoyancy  forces  are  in  opposition  and 

dilution  has  a much  larger  effect  on  density  than  cooling.  Table  1 

gives  the  temperature  change,  aT,  necessary  to  offset  a 1 °/00 

dilution  of  30  °/00  seawater  such  that  the  density  remains  constant  for 

various  temperatures.  At  temperatures  less  than  6°C,  aT  is  so  large 

that  the  water  freezes  before  cooling  enough  to  offset  the  dilution. 

As  the  temperature  increases  aT  decreases,  which  means  that  at  low 

temperatures  the  salinity  dominates  the  density,  but  as  the  temperature 

increases  the  thermal  effect  becomes  more  important.  The  increasing 

effect  of  temperature  results  from  the  almost  parabolic  dependence  of 

density  on  the  temperature  elevation  above  the  temperature  of  maximum 

density,  T . 

m 


10 


Table  1 


S = 30  o/oo 

T(°C) 

0 

2 

4 

6 

8 

10 

20 


S - 1 °/00 
T(°C) 


-23.4 

-13.6 

-9.7 

-7.6 

-6.3 

-5.4 

-3.2 


Temperature  change  necessary  to  offset  the  density  change 
resulting  from  a 1 °/00  decrease  in  salinity  at  30  °/00. 


Dimens : orta  I A s t e 

The  flow  adjacent  to  a melting  vertical  ice  wall  falls  into  the 
category  of  boundary  layer  flow  because  the  flow  is  contained  in  a thin 
region  next  to  the  ice  whose  thickness  is  much  less  than  its  length. 
Schlichting  (1960)  scales  the  equations  governing  the  flow  to  obtain 
relationships  between  nondimensional  parameters  formed  by  combining 
the  molecular  properties  of  the  fluid  with  themselves  and  the  physical 
characteristics  of  the  flow. 

Three  nondimensional  quantities,  the  Schmidt,  Prandtl  and  Lewis 
numbers,  govern  the  distribution  of  sal t, heat, and  momentum  within  the 
boundary  layer.  The  ratio  of  the  thermal  diffusivity  k to  the  saline 
diffusivity  D defines  the  Lewis  number  Le; 


(5) 


11 


The  ratio  of  the  thermal  diffusivity  to  the  kinematic  viscosity  v 
defines  the  Prandtl  number  Pr; 


Pr  * — . 

K 


Finally,  the  ratio  of  the  saline  diffusivity  to  the  kinematic  viscosity 
defines  the  Schmidt  number  Sc; 


For  a laminar  boundary  layer  adjacent  to  a flat  plate,  Schlichting 
(1960,  p.  300)  shows  that  Pr,  Sc  and  Le  give  the  relative  thicknesses 
of  the  velocity,  thermal  and  saline  boundary  layers.  If  6 , 6U  and  6C 
represent  the  boundary  layer  thicknesses  for  velocity,  heat,  and  salt, 


■JT  3 , 


7^-  * *T7Le 

6h 


(8-a) 


( 8-b ) 


(8-c) 


The  Grashof  number,  Gr.and  the  Reynolds  number,  Re,  characterize 
the  nature  of  the  flow  within  the  boundary  layer.  The  Grashof  number 
is  the  ratio  of  buoyant  forces  to  viscous  forces,  defined  as 


g(pw  - pJ*5 
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where  g is  the  acceleration  of  gravity  9.8  m s'2t  pw  is  the  density  of 
water  at  the  ice-water  interface,  pb  is  the  density  of  water  in  the 
far  field,  and  l is  a vertical  length  scale.  For  purely  convective 
flows,  the  magnitude  of  Gr  determines  whether  the  flow  is  laminar  or 
turbulent.  Turner  (1973)  gives  109  as  the  critical  value  of  Gr  where 
the  laminar  flow  becomes  turbulent.  The  Reynolds  number  is  the  ratio 
of  inertial  forces  to  viscous  forces,  defined  as 

uZx 

Re  - ^ (10) 


where  u is  a velocity  characteristic  of  the  convective  motions. 

The  Nusselt  number,  Nu  and  the  Sherwood  number,  Sh,  are  the  non- 
dimensional  heat  or  salt  fluxes  to  the  ice.  Nu  and  Sh  are  the  ratios 
of  the  observed  heat  and  salt  fluxes  to  the  ice  at  a distance  downstream 
to  the  purely  diffusive  fluxes  based  on  the  same  downstream  distance: 


where  FH  is  the  observed  heat  to  the  ice,  F<-  is  the  observed  salt  to 
the  ice,  aT  is  the  temperature  difference  between  interface  and  far- 
field,  and  aS  is  the  salinity  difference  between  interface  and  far-field. 
The  heat  and  salt  transfer  to  the  ice  depends  on  the  flow  past  the  ice; 
buoyant  flow  and  flow  due  to  the  ice  moving  through  the  water  comprise 
the  total  flow. 

Next,  to  determine  if  free  or  forced  convection  dominates  the  flow 
adjacent  to  the  ice  for  oceanographic  conditions,  a scaling  of  the 


momentum  equations  following  Schlichtlng  (1960)  shows  that  free  convec- 
tion becomes  important  when 

Gr  = Re2  or  ^ ! . (12) 

Using  the  definitions  of  Re  and  Gr,  (12)  gives  an  estimate  of  the  speed 
necessary  for  forced  convection  effects  to  equal  free  convection  effects 
as 

u = (9^g1/2  • (13) 

For  oceanic  cases  I assume  that  fresh  water  is  generated  at  the 
ice-water  interface  so  that  a characteristic  value  of  g — is  0.1  m s-2. 
Substitution  of  this  value  and  a value  of  100  m for  Zx  into  Eq.  (13) 
gives 

u • 3 m s-1  . (14) 

In  June  1977  I observed  drogue  studies  next  to  an  iceberg  by  the  U.  S. 
Coast  Guard  that  show  the  iceberg  and  a 75  m deep  window  shade  drogue 
to  separate  at  a rate  of  several  kilometers  per  day  or  -.3  m s"1. 
Therefore,  natural  convection  dominates  the  flow  near  an  iceberg.  For 
the  length  scales  of  pressure  ridges  and  sea  ice  edges,  a flow  of 
1 m s"1  will  overwhelm  the  buoyant  flow  but  a 1 m s'1  flow  is  quite 
high  for  the  oceans.  Hence,  buoyant  motions  remain  important  at  all 
length  scales  and  determine  the  heat  and  salt  fluxes  to  the  ice. 

In  seawater  the  molecular  dlffusivities  of  salt,  heat,  and 
momentum  differ  by  three  orders  of  magnitude.  Appendix  II  gives  the 
values  of  the  molecular  diffusivities  used  in  this  study  and  shows  that 
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the  magnitudes  of  the  diffusi vities  of  salt,  heat,  and  momentum  are, 
respectively,  10"9,  1 0~ 7 , and  10'6  m2  s'1.  For  the  laminar  case, 
tqs.  (Ha,  b,  and  c)  give 


so  that 


-T-  - 0.07  . 
5H 


( 1 5 - a ) 

(15-b) 

(15-c) 


The  large  difference  between  and  6^  shows  that  the  opposing  buoyancy 
forces  act  in  different  regions  of  which  is  larger  than  either  ^ 
or  . The  upward  buoyancy  dominates  the  inner  region  while  the 
downward  buoyancy  dominates  the  outer  region;  the  present  study  found 
that  in  the  laminar  case,  this  buoyancy  distribution  produces  a 
bidirectional  flow. 

The  existence  of  opposing  buoyancy  forces  allows  for  the  defini- 
tion of  two  additional  Grashof  numbers  to  characterize  the  individual 
effects  of  thermal  and  saline  buoyancy.  The  saline  Grashof  number, 

Gr<.,  and  the  thermal  Grashof  number,  Gr^,  are  defined  as 


| 

i 

I 


\ 


' 
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( 16-a) 

* 


(16-b) 
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This  definition  of  Gr^  uses  |a|  because  a changes  sign  near  the  line  of 
maximum  density.  For  seawater 

»< 

Gr  = Gr$  (17) 

because  salinity  dominates  density  and  Gr$>>  Gr-p  As  an  example,  for 

a 200  m deep  iceberg  with  AS  = 35  °/00  and  AT  = 4°C,  Gr$  = 7 x 1017  and  k 

Gr^  = 9 x 101S.  Hence,  the  buoyancy  generated  by  changes  in  salinity 

should  dominate  the  flow.  The  present  study  shows  that  this  is  true 

at  T^  < 20°C  but  for  T^  > 20°C,  the  thermal  effects  dominate.  This 

reversal  results  from  the  increased  temperature  effect  on  density  at 

high  temperatures  and  the  large  value  of  Le. 

Finally,  Gr  exceeds  the  critical  value  of  10'1  for  oceanic  salinities 

i 


and  temperatures  for  vertical  length  scales  greater  than  0.5  m.  Thus, 
the  convection  is  turbulent  for  length  scales  greater  than  0.5  m and 
laminar  for  length  scales  less  than  0.5  m. 


Previous  Work 


This  chapter  reviews  the  previous  work  which  divides  into  four 
main  categories:  ice  melting  experiments,  laminar  boundary  layer  flows, 
transition  from  laminar  flow  to  turbulent  flow  and  turbulent  boundary 
layer  flows.  The  last  three  categories  further  divide  into  experimental 
and  theoretical  studies.  The  laminar  studies  include  both  single  and 
multiple  species  diffusion  while  the  transition  and  turbulent  studies 
consider  only  single  species  diffusion. 

Iae  Melting  Experiments 

The  literature  contains  only  three  experiments  where  natural 
convection  results  from  melting  vertical  ice  walls.  The  first  experi- 
ment by  Sandstrom  (1919)  took  place  in  salt  water,  the  second  took 
place  in  fresh  water. 

Sandstrom  (1919)  placed  large  blocks  of  ice  in  a tank  of  seawater 
to  model  the  circulation  driven  by  melting  ice  sheets  during  the  spring 
in  both  the  North  Atlantic  Ocean  and  Norwegian  fjords.  Figure  2 shows 
the  apparatus  and  the  observed  ice  ablation  and  circulation  in  the  tank. 
The  apparatus  bears  a striking  resemblance  to  the  present  apparatus  in 
that  he  recognized  the  need  to  control  stratification  in  the  tank  and  he 
also  used  dye  to  trace  the  fluid  motions.  Walfrid  Ekman,  a founder  of 
modern  oceanography,  made  miniature  current  meters  to  r.ieasure  the  fluid 
velocities.  At  100  mm  up  from  the  bottom  of  the  ice  block  the  water 
flowed  horizontally  towards  the  Ice,  then  divided,  with  some  water 
flowing  up  next  to  the  ice  and  the  rest  flowing  down  next  to  the  ice. 

A notch  formed  in  the  ice  at  the  jet  level  and  the  ice  above  the  notch 
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melted  faster  than  the  ice  below  the  notch.  A comparison  of  Figure  2 
to  Plate  VI  shows  these  results  are  Identical  to  the  present 
results.  Sandstrom  did  not  describe  an  inner  upward  flow  at  the  bottom 
of  the  ice,  but  since  this  layer  is  only  millimeters  thick,  he  may  not 
have  observed  it.  He  did  describe  cold  fresh  water  flowing  away  from 
the  ice  at  the  surface,  warm  salty  water  flowing  towards  the  ice  at 
depth  and  cold  salty  water  flowing  away  from  the  ice  at  the  bottom  of 
the  tank. 

The  far-field  conditions  for  the  1904  experiment  were  Ta  = +8"C, 

= 30  °/0o  for  3 thermal  driving  of  9.6°C.  Determining  the  average 
wall  salinity  at  = 9.6  from  Figure  8,  Sw  = 4.2  °/00,  and  taking 
the  critical  Gr^  = 1.6  x 108  yields  140  mm  for  the  length  of  the  laminar 
region.  The  additional  flow  of  low  salinity  water  up  the  side  of  the 
ice  from  melting  on  the  bottom  results  in  the  shorter  laminar  region 
because  the  additional  flow  hastens  the  transition  to  turbulence. 
Sandstrom  concluded  from  this  experiment  that  ice  melting  is  important 
in  driving  oceanic  currents. 

Second,  Bendell  and  Gebhart  (1976)  measured  the  heat  transfer  to  a 
vertical  ice  wall  submerged  in  fresh  water  to  Investigate  the  heat 
transfer  when  a density  maximum  exists  in  the  fluid.  When  a vertical 
ice  wall  melts  in  fresh  water,  the  existence  of  the  density  maximum  at 
4°C  provides  the  conditions  necessary  for  bidirectional  flow.  Consider 
the  following  cases:  If  the  ambient  water  temperature  is  greater  than 
8°C,  then  all  of  the  cooled  water  Is  more  dense  and  must  sink.  If  the 
ambient  temperature  is  less  than  4°C,  all  of  the  cooled  water  is  less 
dense  and  must  rise.  Finally,  when  the  ambient  temperature  lies 
between  4 and  8°C,  less  dense  water  is  found  near  the  ice  while  more 


dense  water  Is  found  further  away  from  the  ice,  a situation  which  may 
yield  a bidirectional  flow  next  to  the  ice.  Hence,  the  different 
convection  in  the  three  cases  will  modify  the  heat  transfer  to  the  ice. 

They  determined  the  bulk  heat  transfer  by  the  change  in  weight  of 
a counter-weighted  ice  block  over  the  course  of  an  experiment.  They 
found  the  heat  transfer  to  be  a minimum  at  5.5°C.  Two  thermocouples, 
one  at  each  end  of  the  ice  block  determined  the  direction  of  flow. 

When  the  upper  thermocouple  was  colder  than  the  lower  one,  they  assumed 
the  flow  to  be  upward;  for  the  opposite  situation,  they  assumed  the 
flow  to  be  downward.  They  found  upflow  when  the  ambient  temperature 
was  < 5.6°C  and  downflow  for  warmer  temperatures.  No  detailed 
description  of  the  flow  fields  was  given;  however,  they  assumed  in  all 
cases  that  the  flow  was  unidirectional. 

Third,  Huppert  and  Turner  (1978)  placed  0.25  m long  ice  slabs  in 
linearly  stratified  salt  water  at  20°C.  The  introduction  contains  a 
discussion  of  this  experiment. 

Uminai'  Studies 

This  review  begins  with  the  laminar  single  species  flow;  prior  to 
1954  [as  described  in  Gebhart  and  Pera  (1971)],  investigators  integrated 
the  boundary  layer  equations  with  assumed  velocity  and  species  profiles. 
These  integral  techniques  are  limited  to  Pr  » 1 ; when  Pr  is  not  near  1, 
the  boundary  layer  thickness  for  heat  (salt)  differs  from  the  momentum 
boundary  layer  thickness  and  these  integral  techniques  fail. 

Ostrach  (1953)  performed  a similarity  transformation  of  the  laminar 
boundary  layer  equations  and  then  numerically  integrated  the  resulting 
ordinary  differential  equations  with  Pr  varying  from  .01  to  1000  to 
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determine  the  effect  of  Pr  on  the  flow.  In  this  study  the  blowing 
velocity  at  the  wall  equaled  zero.  Ostrach  found  that  on  Increasing 
Pr  decreases  the  thickness  of  the  thermal  layer  when  compared  to  the 
thickness  of  the  velocity  layer  and  also  decreases  the  convective 
velocities.  The  velocities  decrease  because  the  buoyant  driving  is 
contained  close  to  the  wall  where  the  frictional  influence  of  the  wall 
is  greatest.  The  experimental  work  cited  in  Ostrach  agrees  with  the 
theory  for  Pr  up  to  100  but  no  experiments  have  been  performed  at  Pr  < 100. 

LaPadula  and  Mueller  (1970)  performed  a theoretical  investigation 
of  a solid  freezing  onto  a vertical  surface.  In  their  analysis  they 
included  the  effect  of  suction  or  blowing  at  the  interface.  They  found 
that  for  freezing,  the  mass  flux  towards  the  wall,  or  suction,  enhances 
the  heat  transfer  to  the  wall  while,  for  melting,  the  mass  flux  away 
from  the  wall,  or  blowing,  reduces  the  heat  transfer.  These  effects 
increased  with  increasing  Pr.  These  results  indicate  that  the  blowing 
due  to  ice  melting  must  be  included  in  our  analysis. 

When  multiple  species  diffuse  across  the  boundary  layer,  the 
resulting  buoyancy  forces  either  reinforce  or  oppose  each  other.  In 
multiple  species  diffusion  the  first  species  is  usually  heat,  and  the 
second  is  a dissolved  chemical  compound.  The  most  interesting  case, 
which  corresponds  to  the  oceanic  case,  occurs  when  the  buoyancy  forces 
oppose  each  other  and  Le  » 1 , which  will  both  separate  the  individual 
buoyancy  forces  within  the  boundary  layer  and  yield  bidirectional  flow. 

In  two  previous  experimental  studies  Adams  and  McFadden  and  Den 
Bouter  used  a heated  plate  to  generate  an  upward  force  while  sublimation 
of  a solid  from  the  plate  generated  the  downward  force.  When  water  was 
the  ambient  fluid,  electrochemical  means  generated  the  downward  force. 
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These  studies  concentrated  on  measuring  the  heat  and  mass  fluxes  at  the 
wall  for  a known  wall  temperature  or  measuring  the  wall  temperature  and 
mass  flux  for  a known  heat  flux.  In  both  of  these  studies  an  adequate 
description  of  the  flow  field  is  missing. 

First,  Adams  and  McFadden  (1966)  performed  an  experiment  with 
opposing  body  forces  in  air  by  using  a sublimating  plate  of  p-dichoro- 
benzene  to  produce  the  downward  force,  and  heat  to  produce  the  upward 
force.  In  their  case  Pr  = .73  and  Sc  = 2.23.  For  opposing  body  forces, 
the  heat  and  mass  transfer  results  were  10  to  15%  lower  than  the 
theoretical  predictions  given  by  integral  techniques  and  no  discussion 
of  the  flow  field  was  given. 

Second,  Den  Bouter  et  al.  (1968)  measured  the  heat  and  mass  trans- 
fer rates  in  an  aqueous  solution  of  copper  sulphate  to  verify  McFadden 's 
results.  For  body  forces  in  the  same  direction,  the  transfer  rates 
agreed  with  the  analytical  predictions.  However,  when  the  buoyancy 
effects  opposed  each  other,  the  observations  deviated  from  the  theory. 
The  flow  observations  were  stated  in  one  sentence:  "The  flow  was  either 
up,  down  or  oscillatory";  the  method  of  flow  visualization  was  not 
discussed. 

Previous  to  1971,  theoretical  work  on  the  effect  of  multiple  body 
forces  on  the  convection  received  black  box  treatment  in  that  the  best 
correlation  between  Nu  or  Sh  and  some  combination  of  Gr-p  Gr^,  Pr,  and 
Sc  was  sought.  Gebhart  and  Pera  (1971)  give  a thorough  review  of  the 
literature  and  discuss  the  shortcomings  of  previous  theories,  then  per- 
form a numerical  investigation  of  the  laminar  equations  on  the  assump- 
tion of  self-similar  flows.  They  present  solutions  for  a wide  range 
of  Prandtl  and  Schmidt  number  in  cases  where  the  buoyancy  effects 


reinforce  or  oppose  each  other.  In  the  case  where  Pr  = Sc  or  Le  = 1, 
the  various  species  diffuse  at  the  same  rate  so  the  buoyancy  results 
from  the  diffusion  of  a single  species  whose  effect  on  density  equals 
the  combined  species  effect  on  density.  More  interesting  situations 
arise  when  Pr  and  Sc  differ.  For  reinforcing  buoyancy  forces,  Gebhart 
and  Pera  show  that  the  boundary  layer  is  more  energetic  and  narrower 
than  a single  species  boundary  layer. 

For  opposing  buoyancy  forces,  the  boundary  layer  is  wider  and  less 
energetic  than  the  single  species  boundary  layer.  The  possibility  of  a 
bidirectional  flow  exists  when  Le  is  very  small  or  very  large;  however, 
Gebhart  and  Pera  do  not  present  solutions  for  oceanic  values  of  Sc  nor 
for  density  anomalies  that  would  be  found  when  seawater  is  in  contact 
with  fresh  ice.  Their  similarity  formulation  only  allows  the  boundary 
layers  to  thicken  and  accelerate  in  the  downstream  direction. 


Transition  to  Turbulence 

In  nature  laminar  flows  quickly  become  turbulent  through  the  growth 
of  infinitesimal  disturbances  within  the  laminar  flow.  A stability 
analysis  of  a laminar  flow  proceeds  by  perturbing  the  mean  laminar  flow 
and  then  determining  the  growth  rates  of  the  perturbations.  Because  no 
adequate  analysis  of  the  laminar  bidirectional  flow  exists,  the  stabil- 
ity of  this  flow  has  not  been  determined.  However,  the  stability 
properties  of  a single  species  boundary  layer  yields  insight  into  the 
transition  of  a bidirectional  boundary  layer.  Hieber  and  Gebhart  (1971, 
a,  b)  analyzed  the  stability  of  the  flow  generated  by  a uniform  heat 
flux  from  a flat  vertical  plate.  This  problem  Is  not  identical  to  the 
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present  experiments  but  the  physical  processes  are  similar  and  the 
authors  describe  a method  that  can  be  used  to  convert  from  a uniform 
flux  case  to  a uniform  wall  temperature  case.  They  identified  two  modes 
of  instability:  one  associated  with  the  instability  that  arises  from 
the  shape  of  the  velocity  profile,  an  Orr-Sommerfeld  instability,  and 
the  other  related  to  the  coupling  of  momentum  and  heat  through  the 
buoyancy  term  in  the  momentum  equation. 

For  small  Pr,  -0.01  < Pr  < 0.1,  the  coupling  effect  has  little  effect 
on  the  stability  or  the  position  of  the  neutral  curve.  The  neutral 
curve  is  the  locus  of  points  on  a plot  of  perturbation  frequency 
versus  a modified  Grashof  number 

• 08) 

such  that  the  imaginary  part  of  the  frequency  is  zero  hence,  the  waves 
do  not  grow  or  decay.  In  small  Prandtl  number  flows,  the  momentum 
layer  is  narrower  than  the  thermal  layer  so  a perturbation  in  the 
velocity  field  only  advects  a small  thermal  change  that  produces  an 
insignificant  buoyancy  perturbation. 

For  Pr  ^ 1 the  buoyant  coupling  dominates  the  instability  mecha- 
nism, and  the  neutral  curve  moves  towards  lower  modified  Grashof  numbers. 
For  large  Pr,  the  momentum  layer  is  much  wider  than  the  thermal  layer, 
so  that  a velocity  perturbation  will  advect  a significant  change  in 
temperature  because  the  local  temperature  gradient  is  large.  The 
changes  in  the  temperature  field  can  produce  momentum  in  two  ways: 

First,  if  the  buoyant  thermal  anomaly  is  advected  away  from  the  wall, 
it  will  accelerate  the  flow  since  it  is  in  a region  of  lower  buoyancy. 
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Second,  less  buoyant  water  that  Is  advected  towards  the  wall  will 
increase  its  buoyancy  since  It  has  moved  into  a region  of  lower 
density.  The  expected  wavelength  of  the  most  unstable  disturbance  for 
the  coupled  mode  is  of  order  xPr-1/2,  while  for  the  uncoupled  mode  the 
most  unstable  wavelength  is  of  order  xPr'1,  where  x is  the  distance 
downstream  from  the  tip  of  the  plate. 

In  an  experimental  study  Jalurla  and  Gebhart  (1974)  studied  the 
transition  of  the  natural  laminar  convective  boundary  layer  in  water 
for  a uniform  heat  flux  from  a flat  vertical  plate.  They  made  combined 
temperature  and  velocity  measurements  using  a copper-constantan  thermo- 
couple and  a hot  wire  anemometer  that  traversed  the  boundary  layer. 

The  results  show  that  the  velocity  and  thermal  boundary  layers  increase 
in  thickness  as  both  the  distance  from  the  tip  and  the  heat  flux 
increase.  The  changes  in  the  transition  region  correlate  with  G*/xn 
where  n is  approximately  0.5.  The  factor  of  l/xn  appears  because  the 
growth  of  the  disturbances  is  probably  related  to  the  local  values  of 
the  momentum  and/or  energy  fluxes,  quantities  that  do  not  correlate 
on  G*. 

The  authors  suggest  a new  parameter,  E,  to  characterize  the 
transition  to  turbulence  based  on  the  boundary  layer  kinetic  energy 
flux.  E is  defined  as  the  nondimensional  kinetic  energy  flux  of  the 
boundary  layer  based  on  the  laminar  flow: 


E 


2/15 


(19) 
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Using  the  conversion  given  by  the  authors  to  convert  from  uniform  flux 
to  isothermal  conditions,  the  value  of  E for  n\y  experiments  is  18, 
with  x equal  to  the  length  of  the  laminar  regime.  This  value  agrees 
remarkably  well  with  values  of  E from  13  to  24  from  experiments  cited 
in  Jaluria  and  Gebhart  performed  at  Prandtl  numbers  0.7  to  12.  This 
suggests  that  the  transition  to  turbulence  does  depend  on  the  available 
kinetic  energy  from  which  the  small  disturbances  can  draw  their  energy. 

These  past  studies  show  that  the  stability  properties  of  a single 
species  flow  are  quite  well  known  and  that  similar  stability 
characteristics  exist  for  the  bidirectional  case.  In  the  bidirectional 
case  the  inner  saline  layer,  with  its  large  buoyant  forcing  and  high 
Schmidt  number,  dominates  the  stability  of  the  flow  while  the  outer 
thermal  layer  modifies  the  stability  characteristics  of  the  flow  by 
changing  the  velocity  distribution  across  the  boundary  layer.  The 
effect  of  blowing  at  the  wall  on  the  stability  has  not  been  investigated 
and  it  may  have  a significant  effect,  especially  at  high  Schmidt  number. 

Turku  ten  t S tud  i cu- 
stodies of  the  turbulent  convective  boundary  layer  are  scarce;  the 
previous  experimental  studies  do  not  represent  the  complexities  of  ice 
melting  in  seawater  and  the  previous  theoretical  analyses  use  results 
from  forced  convection  theory  and  experiments.  Previous  experiments  had 
heat  diffusing  in  air  or  water,  for  Pr  = 0.73  or  10.0,  where  the  result- 
ing thermal  anomaly  creates  the  buoyancy.  An  electrically  heated  metaT 
plate  either  produces  a uniform  heat  flux  to  the  ambient  fluid  or  main- 
tains a uniform  wall  temperature.  Neither  theory  nor  experiment 
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Second,  less  buoyant  water  that  Is  advected  towards  the  wall  will 
increase  its  buoyancy  since  It  has  moved  into  a region  of  lower 
density.  The  expected  wavelength  of  the  most  unstable  disturbance  for 
the  coupled  mode  is  of  order  xPr'1/2,  while  for  the  uncoupled  mode  the 
most  unstable  wavelength  Is  of  order  xPr-1,  where  x is  the  distance 
downstream  from  the  tip  of  the  plate. 

In  an  experimental  study  Jaluria  and  Gebhart  (1974)  studied  the 
transition  of  the  natural  laminar  convective  boundary  layer  in  water 
for  a uniform  heat  flux  from  a flat  vertical  plate.  They  made  combined 
temperature  and  velocity  measurements  using  a copper-constantan  thermo- 
couple and  a hot  wire  anemometer  that  traversed  the  boundary  layer. 

The  results  show  that  the  velocity  and  thermal  boundary  layers  increase 
in  thickness  as  both  the  distance  from  the  tip  and  the  heat  flux 
increase.  The  changes  in  the  transition  region  correlate  with  G*/xn 
where  n is  approximately  0.5.  The  factor  of  l/xn  appears  because  the 
growth  of  the  disturbances  is  probably  related  to  the  local  values  of 
the  momentum  and/or  energy  fluxes,  quantities  that  do  not  correlate 
on  G*. 

The  authors  suggest  a new  parameter,  E,  to  characterize  the 
transition  to  turbulence  based  on  the  boundary  layer  kinetic  energy 
flux.  E is  defined  as  the  nondlmensional  kinetic  energy  flux  of  the 
boundary  layer  based  on  the  laminar  flow: 


E - G* 
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Using  the  conversion  given  by  the  authors  to  convert  from  uniform  flux 

to  isothennal  conditions,  the  value  of  E for  my  experiments  is  18,  W 

with  x equal  to  the  length  of  the  laminar  regime.  This  value  agrees 

! 

remarkably  well  with  values  of  E from  13  to  24  from  experiments  cited 
in  Jaluria  and  Gebhart  performed  at  Prandtl  numbers  0.7  to  12.  This 
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suggests  that  the  transition  to  turbulence  does  depend  on  the  available 
kinetic  energy  from  which  the  small  disturbances  can  draw  their  energy. 

These  past  studies  show  that  the  stability  properties  of  a single 
species  flow  are  quite  well  known  and  that  similar  stability 
characteristics  exist  for  the  bidirectional  case.  In  the  bidirectional 
case  the  inner  saline  layer,  with  its  large  buoyant  forcing  and  high 
Schmidt  number,  dominates  the  stability  of  the  flow  while  the  outer 
thermal  layer  modifies  the  stability  characteristics  of  the  flow  by 
changing  the  velocity  distribution  across  the  boundary  layer.  The 
effect  of  blowing  at  the  wall  on  the  stability  has  not  been  investigated 
and  it  may  have  a significant  effect,  especially  at  high  Schmidt  number. 

Turbulent  Studies 

Studies  of  the  turbulent  convective  boundary  layer  are  scarce;  the 
previous  experimental  studies  do  not  represent  the  complexities  of  ice 
melting  in  seawater  and  the  previous  theoretical  analyses  use  results 
from  forced  convection  theory  and  experiments.  Previous  experiments  had 
heat  diffusing  in  air  or  water,  for  Pr  = 0.73  or  10.0,  where  the  result- 
ing thermal  anomaly  creates  the  buoyancy.  An  electrically  heated  metal 
plate  either  produces  a uniform  heat  flux  to  the  ambient  fluid  or  main- 
tains a uniform  wall  temperature.  Neither  theory  nor  experiment 
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includes  the  effect  of  blowing  at  the  solid  boundary  or  opposing 
buoyancy  forces. 

Eckert  and  Jackson  (1951)  analyzed  the  turbulent  boundary  layer 
for  a uniform  wall  temperature  and  no  blowing  by  assuming  shapes  for 
the  velocity  and  temperature  profiles  and  then  applying  the  integral 
techniques  discussed  previously.  They  modeled  the  turbulence  by 
taking  the  wall  stress  from  measurements  in  forced  convection  because 
the  flow  near  the  wall  should  have  the  same  characteristics  in  each 
situation.  The  results  give  a boundary  layer  that  grows  as  x0,7  and  a 
maximum  velocity  that  increases  as  x1/,2t  where  x is  the  downstream 
coordinate.  The  average  Nu-dependence  on  Gr  and  Pr  was 
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Nu  - (GrPr)2/5  , (20) 

which  agrees  with  the  few  previous  experimental  studies  conducted  in  air. 

Kato  et_  alL  (1968)  also  analyzed  the  turbulent  free  convection  as 
a special  form  of  forced  convection.  They  took  eddy  diffusivitits  from 
forced  convection  results  and  assumed  that  the  eddy  diffusivity  of  heat 
equals  that  of  momentum.  The  eddy  diffusivity  was  zero  at  the  wall, 
remained  small  over  a viscous  sublayer  scale,  then  increased  linearly 
over  the  entire  boundary  layer  thickness.  Equation  (21)  gives  the  eddy 
diffusivity 

e/v  = .4y,[l-exp(-.0017  y'2)]  , (21) 


where  e is  the  eddy  diffusivity  and  y'  is  the  dimensionless  distance 
from  the  wall.  For  Pr  = 1,  y'  « 103  approximates  the  outer  edge  of  the 
boundary  layer  and  e/v  = 400.  The  integral  momentum  and  heat  equations 
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are  solved  at  several  downstream  positions,  corresponding  to  Gr  = 1010 
to  10lj,  by  a trial  and  error  method.  At  each  Gr,  an  assumed  tempera- 
ture distribution  determined  a velocity  profile;  the  new  velocity 
profile  then  produced  a new  temperature  distribution;  this  iterative 
process  continued  until  no  change  occurred  in  the  profiles.  The 
velocity  profiles  agree  with  those  of  Eckert  and  Jackson  and  with  the 
experimental  data  out  to  Gr  x Pr  equal  to  1012.  At  Pr  equal  to  10-  to 
10\ 

Nux  « Gr° • 31  , (22) 

where  Nux  is  the  local  Nusselt  number.  This  corresponds  to  a decrease 
in  the  heat  flux  to  the  wall  as  the  downstream  distance  increases.  The 
authors  only  give  velocity  and  temperature  profiles  for  Pr  = 1 and  omit 
the  effects  of  blowing  or  suction  at  the  wall  and  the  opposing  buoyancy 
forces,  effects  which  the  present  study  shows  to  be  important. 

Experimental  studies  of  the  turbulent  boundary  layer  divide  into 
two  categories,  uniform  vertical  wall  temperature  or  uniform  vertical 
heat  flux  at  the  wall.  The  results  show  that  the  associated  convection 
is  different  for  each  category.  In  either  case,  the  effects  of  blowing 
and  opposing  buoyancy  forces  are  not  present. 

Cheesewright  (1968)  measured  the  mean  velocity  and  temperature 
profiles  and  the  local  heat  transfer  rates  in  the  turbulent  flow 
produced  by  electrical  heating  of  an  aluminum  plate  2.75  m high  and 
0.6  m wide.  The  apparatus  maintained  vertically  uniform  wall  tempera- 
ture on  the  plate.  Platinum  resistance  thermometers  measured  both  the 
temperature  and  the  velocity  fields;  the  instrument  calibration  took 
place  in  the  laminar  boundary  layer  at  the  bottom  of  the  plate. 
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The  transition  from  laminar  to  turbulent  flow  occurred  at 
2 x 109  to  1010and  the  turbulent  measurements  were  carried  out  to 
Gr  = 1011.  The  measured  heat  transfer  agrees  with  the  theory  of  Eckert 
and  Jackson;  however,  the  assumed  shapes  for  velocity  and  temperature 
do  not  agree  with  the  experiment.  The  assumed  velocity  profile  is  too 
narrow  and  too  fast  when  compared  to  the  measured  velocity  profile  and 
the  assumed  temperature  anomaly  profile  decreases  too  slowly  as  the 
distance  from  the  wall  increases  when  compared  to  measured  profiles. 

When  flow  changes  from  laminar  to  turbulent,  the  turbulence  steepens 
the  temperature  gradient  near  the  wall  and  diffuses  the  temperature 
anomaly  further  from  the  wall.  The  data  near  the  wall  best  correlate  on 

n = y/x  Gr-1  (23) 

which  indicates  a different  turbulent  mechanism  acting  in  each  part  of 
the  boundary  layer. 

Vliet  and  Liu  (1969)  conducted  an  experiment  in  water  with  a 
vertical  plate  dissipating  a vertically  uniform  heat  flux.  They 
measured  the  mean  temperature  and  velocity  profiles  next  to  a 1.25  m 
long  plate  with  a copper-constantan  thermocouple  and  a hot  film 
anemometer.  The  water  temperature  varied  from  7°C  to  50°C. 

The  wall  temperature  distributions  as  measured  by  a series  of 
thermocouples  on  the  surface  of  the  plate  show  that  the  wall 
temperature  is  coldest  at  the  bottom  of  the  plate,  increases  in  the 
laminar  region  and  then  rapidly  decreases  in  the  transition  region. 
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In  the  turbulent  region  the  wall  temperature  increased  with  distance 
up  the  plate,  indicating  a decrease  in  the  heat  transfer  coefficient; 
the  decrease  was  proportional  to  x~*12.  Also  the  wall  temperatures  in 
the  turbulent  region  fluctuated  with  time,  a phenomenon  that  was  not 
observed  in  the  ice  interior  in  the  present  experiments. 

A comparison  of  the  laminar  and  turbulent  heat  fluxes  shows  that 
the  local  Nusselt  number  Nux  dependence  on  the  local  modified  Grashof 
number  Gr*,  where  Gr*  is  Gr^  multiplied  by  Nux  has  almost  the  same 
power  law  dependence  in  each  case,  0.2  versus  0.22.  This  indicates 
that  the  turbulent  flow  exhibits  strong  laminar  characteristics.  The 
transition  to  turbulence  occurs  in  a range  of  Gr*  equivalent  to  Grx  of 
109  to  1010,  1 to  2 orders  of  magnitude  greater  than  the  constant  wall 
temperature  case,  which  indicates  that  molecular  diffusion  processes 
remain  important. 

When  compared  to  the  data  of  Cheesewright,  the  measured  velocity 
and  temperature  profiles  are  in  good  agreement;  however,  the  transposed 
heat  transfer  data  do  not  agree.  The  authors  state  that  the  disagree- 
ment results  from  the  two  modes  of  heating  which  produce  slightly 
different  fluid  dynamics  in  each  case.  The  velocity  measurements  found 
the  maximum  velocity  within  7 mm  of  the  wall  and  a sharp  division 
between  the  outer  quiescent  fluid  and  the  inner  turbulent  fluid  that 
made  it  difficult  to  estimate  an  average  boundary  layer  thickness. 

From  their  measurements  the  authors  deduced  that  a one-parameter 
similarity  may  be  applicable  to  the  mean  fields. 
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In  summary  this  literature  review  shows  that  previous  work  on 
natural  convection  is  not  applicable  to  the  convection  generated  by  i ■ 

melting  vertical  ice  walls  in  seawater  because  these  studies,  except 
for  that  of  Sandstrom,  took  place  at  Gr  and  Sc  of  at  least 
two  orders  of  magnitude  less  than  the  expected  values  for  oceanic 
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situations  and  the  effect  of  blowing  is  omitted.  The  theoretical  ! 

laminar  and  turbulent  studies  suffer  from  the  same  drawbacks,  although  i 

LaPadula  did  investigate  the  effect  of  blowing  or  suction  at  the  wall 
for  laminar  single  species  convection.  Also  the  theoretical  turbulent 
studies  use  eddy  diffusivities  from  forced  convection  experiments. 

I 

k 

The  Present  Study  ] 

* 

I 

To  investigate  the  convection  that  occurs  when  vertical  ice  walls 
melt  in  seawater,  I performed  a series  of  laboratory  experiments  to 
characterize  the  convective  motions  and  determine  the  dependence  of 
T , S,  and  the  melt  rate  on  oceanic  values  of  T and  S . The  experi- 
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ments  show  that  the  convection  consists  of  both  laminar  and  turbulent 

I 

flows  and  oceanographically  the  most  important  flow  is  an  upward  flow- 
ing turbulent  boundary  layer  driven  by  low  salinity  water. 

I then  numerically  solved  the  equations  of  motion  governing  the 
turbulent  convection  under  the  assumption  of  self-similarity.  A 
spatially  varying  eddy  diffusivity  parameterized  the  turbulent  trans- 
port of  momentum,  heat,  and  salt  and  the  experimental  data  determined 
the  turbulent  parameters.  In  solving  these  equations  I retained  the 
effect  of  opposing  buoyancy  forces  and  blowing  or  mass  influx  at 
the  ice-water  interface, and  I found  that  these  effects  have  a profound 
influence  on  the  numerical  solutions. 
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In  an  application  of  the  laboratory  and  theoretical  results, 

I conducted  a field  study  adjacent  to  an  iceberg  in  the  North  Atlantic 

Ocean  in  cooperation  with  the  United  States  Coast  Guard  and  the 

International  Ice  Patrol.  The  field  study  consisted  of  a series  of 

STD  casts  around  an  iceberg  and  the  results  show  that  the  ice  melt 

modifies  the  salinity  and  temperature  of  the  water  in  the  upper  30  m 

within  200  m of  the  iceberg,  which  supports  the  concept  of  an  upward 

flowing  boundary  layer  adjacent  to  the  ice.  r 
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The  Laboratory  Experiment 

I 1 

Apparatus 

In  my  experimental  study  of  the  ablation  of  a vertical  ice  wall  in  < 

seawater  the  important  laboratory  parameters  are  the  salinity  and 

] 

temperature  of  the  ambient  water  and  the  length  of  the  ice  wall. 

| 

Oceanic  temperature  and  salinity  conditions,  as  well  as  1 m vertical 
length  scales,  are  easily  modeled  in  the  laboratory;  these  results  may 
be  extrapolated  to  vertical  length  scales  greater  than  1 m. 

The  experiments  took  place  in  a clear  acrylic  tank  that  measures 
0.45  m wide,  1.2  m long,  and  1.25  m deep  which,  when  filled,  contains 
approximately  650  kg  of  water.  I produced  water  at  oceanic  salinities 

l . 

by  mixing  food  grade  NaCl  into  the  water.  I then  determined  the 
resultant  salinity  with  either  an  American  Optical  or  Endeco  refracto- 

. 

meter,  respectively,  accurate  to  1.0  °/00  or  0.2  °/00.  I obtained 
oceanic  temperatures  by  placing  the  filled  tank  in  one  of  the  cold 
rooms  of  the  Sea  Ice  Laboratory,  Department  of  Oceanography,  University 
of  Washington,  and  setting  the  room  temperature  to  the  required  value. 

The  water  than  reached  thermal  equilibrium  with  the  air;  a process 
which  took  several  days  as  a result  of  the  large  mass  of  water.  The 
cold  rooms  are  large  walk-in  freezers  with  controlled  temperatures 
which  range  from  -40°C  to  +5°C.  The  temperature  within  each  room 
cycles  around  the  set  temperature  with  an  amplitude  of  1°C  and  a period 
of  40  minutes.  For  temperatures  higher  than  5°C,  I shut  off  the 
refrigeration  and  let  the  room  and  water  warm  to  the  required 
temperature. 

Before  placing  the  ice  in  the  tank,  1 insured  homogeneous  tempera- 


ture and  salinity  conditions  by  using  a stirring  motor  to  vigorously 


33 


agitate  the  water  and  then  letting  the  water  come  to  rest.  In  addition, 
I removed  the  ice  from  cold  storage  and  let  it  warm  to  a temperature  of 
approximately  -1°C  to  prevent  thermal  fracture  of  the  ice  slab  when  it 
was  immersed  in  the  water. 

In  the  experiments  I used  sheets  of  bubble-free  ice  with  an 
instrument  array  frozen  within  the  ice  to  model  the  melting  process. 

The  ice  sheets  measured  approximately  0.2  m wide,  0.1  m thick  and  varied 
in  length  from  0.5  m to  1.2  m.  I grew  bubble-free  ice  in  an  insulated 
tank,  0.4  m x 0.6  m x 0.4  m,  that  has  a cold  plate  on  the  bottom  and  two 
laboratory  stirring  motors  mounted  on  the  top.  As  the  ice  grew,  the  air 
bubbles  that  formed  at  the  ice-water  interface  were  swept  away  by  the 
water  motion  created  by  the  stirring  motors.  When  the  ice  thickness 
reached  half  the  desired  value  or  50  mm,  the  array  was  placed  in  the 
tank  and  the  ice  froze  around  the  array.  When  the  ice  reached  the 
desired  thickness,  I set  the  temperature  in  the  cold  plate  at  +2°C;  the 
subsequent  small  amount  of  melting  freed  the  ice  from  the  bottom  of  the 
tank  so  that  it  floated  to  the  surface  where  it  was  easily  removed.  As 
the  ice  grown  in  this  tank  measured  only  0.6  m long,  to  obtain  1.2  m 
long  ice  I joined  two  such  lengths  of  ice  by  butting  the  ends  together 
in  a cold  room  at  -15°C  and  seeping  water  into  the  joint.  For  these 
cases,  the  array  lies  only  in  the  upper  portion  of  the  ice.  I also 
performed  two  experiments  using  ice  that  had  a C0^  content  of  5%  or 
1.5T  by  volume. 

Before  storing  the  ice  or  at  the  beginning  of  an  experiment, 

I used  a laboratory  hot  plate  to  smooth  one  side  of  the  ice  slab  to  a 
surface  that  was  planar  to  within  2 mm/m.  As  Figure  3 shows,  I also 
cut  a bevel  on  the  ice  bottom  to  prevent  bottom  melting  from  interfering 


Figure  3.  The  experimental  apparatus. 
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with  side  melting  by  directing  positively  buoyant  water  formed  near  the 
bevel  up  the  back  side  of  the  ice.  To  hold  the  ice  in  place  a 12.7  mm 
diameter  aluminum  rod  screwed  into  the  plastic  array  and  then  was 
clamped  to  a similar  rod  across  the  top  of  the  tank.  A laboratory 
swivel  clamp  allowed  the  ice  face  to  be  aligned  parallel  to  gravity. 

The  instrument  array  frozen  in  the  ice  consisted  of  ten  precision 
thermistors  and  four  ablation  probes  mounted  on  an  acrylic  plastic 
support.  Figure  4 gives  the  arrangement  of  the  thermistors  and 
ablation  probes  on  the  support  which  was  0.5  m long  with  a cross  section 
of  12.7  mm  x 12.7  mm.  To  minimize  wake  effects  after  the  thermistors 
melt  through  the  ice,  no  two  thermistors  lie  on  the  same  vertical  line. 
Each  probe  protrudes  15  or  20  mm  beyond  the  plastic  support  so  that  the 
probes  melt  through  the  ice  at  different  times  during  an  experiment. 
Thermistors  3 and  4,  situated  at  the  same  horizontal  level,  were  used  to 
determine  the  cross-ice  variation  of  the  interface  temperature. 

The  thermistors  are  Yellow  Springs  Instruments  precision  thermi- 
stors, #44007,  with  a time  constant  of  1 s and  a dissipation  constant 
of  8 mw  deg-1  in  a well-stirred  oil  bath.  The  thermistors  are  roughly 
spherical  with  a diameter  of  2.4  mm.  To  make  the  thermistors  watertight 
I insulated  the  leads  with  a 15  to  20  mm  length  of  snug-fitting  shrink 
tubing  filled  with  self-curing  silicone  rubber.  I then  heated  the 
shrink  tubing  with  a hot  air  gun  which  caused  the  tubing  to  shrink  and 
extrude  the  excess  rubber.  To  protect  the  thermistor  from  temperature 
extremes  during  the  heating  process,  I immersed  the  exposed  thermistor 
tip  to  a depth  of  1 mm  in  fresh  water.  After  curing  for  several  days 
the  insulation  became  watertight.  I then  mounted  each  sensor  on  a 5 mm 
wide,  1 mm  thick  plastic  strip  of  varying  length  and  fastened  this  strip 
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® THERMISTORS 
■ ABLATION  PROBES 


Figure  4.  The  thermistor 


array. 


37 


to  the  main  support  at  the  required  level. 

The  four  ablation  probes  mounted  on  the  array  are  strips  of  plastic 
rule  measuring  35mmx3mmxlirm  mounted  in  a vertical  plane  perpendi- 
cular to  the  ice  surface  at  the  positions  shown  in  Figure  4.  Within 
the  vertical  plane,  each  probe  tilted  downward  at  an  angle  of  10°  so  that 
the  real  image  and  the  reflected  Image  of  the  probe  extending  out  from 
the  ice  formed  an  inverted  "V".  The  probes  were  viewed  at  an  angle  of  5° 
from  the  ice  plane  because  the  large  change  in  index  of  refraction  result- 
ing from  the  low  salinity  water  at  the  interface  hindered  sighting 
directly  along  the  ice  plane.  To  determine  the  ablation  rate  I periodi- 
cally measured  the  horizontal  distance  from  the  tip  of  the  probe  to  the 
apex  of  the  "V"  with  a Keuffel  and  Esser  measuring  transit  with  an  over- 
all range  of  5 mm  that  is  accurate  to  0.05  mm.  The  graduations  on  the 
plastic  rule  served  as  reference  points  for  measuring  lengths  greater 
than  5 rnn. 

Because  the  far-field  temperature  and  salinity  determine  the  ice 
melting  processes,  I measured  the  far-field  condition  with  an  array 
consisting  of  three  thermistors  and  three  surgical  tubes.  The  three 
thermistors  are  the  same  as  used  on  the  ice  array,  and  they  are  sealed 
in  the  same  manner.  I mounted  them  on  a 12.7  mn  x 12.7  itm  acrylic 
plastic  rod  1.2  m long;  0.25  m separates  each  sensor  with  the  bottom- 
most thermistor  at  the  bottom  of  the  rod.  Each  surgical  tube,  which 
was  taped  to  the  support,  terminated  at  one  of  the  thermistors  and  with 
a syringe  I drew  water  samples  for  later  salinity  determination  by 
either  refractometer  previously  described.  Figure  3 shows  the  place- 


ment of  the  ice  and  far-field  array  in  the  tank.  I placed  this  array 
0.5  m from  the  ice  surface  such  that  the  far-field  array  and  the  array 
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within  the  ice  spanned  the  same  vertical  distance  inside  the  tank. 

During  an  experiment  I recorded  the  thermistor  resistances  with 
a Non-Linear  Systems  data  acquisition  system  that  scanned  each  of  the 
thermistors  once  every  minute,  measured  the  resistances  with  a 
Keithly  Model  160  Digital  Multimeter  and  recorded  the  data  on  magnetic 
tape.  The  Keithly  meter  has  an  accuracy  of  * 0 . 01  k-Ohm,  which 
corresponds  to  an  accuracy  in  temperature  of  approximately  *0.01CC. 
Figure  5 gives  a characteristic  curve  of  the  resistance  data  for  a 
typical  experiment.  The  flat  section  labeled  "A"  gives  the  ice-water 
interface  temperature.  The  ice  melted  past  the  sensor  at  the  point 
labeled  "B"  and  the  outer  edge  of  boundary  layer  passed  the  sensor 
at  the  point  labeled  "C".  The  interval  of  time  from  the  thermistor 
meltout  to  the  outer  edge  of  the  boundary  layer  along  with  the  melting 
rate  gives  an  estimate  of  the  thermal  boundary  layer  thickness. 

To  visualize  the  convective  motions  generated  by  the  melting  ice 

1 illuminated  small  reflecting  flakes  suspended  in  the  salt  water  with 
a vertically  oriented  plane  of  light.  The  flakes  are  finely  ground 
mother-of-pearl,  obtained  from  the  Mearl  Corporation,  41  East  42nd 
Street,  New  York  City,  New  York  10017,  with  the  trade  name  "AA-Mearl 
Maid  Flakes".  The  flakes  are  in  paste  form  and  1 ml  of  paste  emulsified 
in  water  insured  a dense  enough  distribution  of  particles  in  the  tank  to 
trace  the  flow  field.  Because  the  settling  time  of  the  flakes  is  of 
order  days,  the  flakes  accurately  followed  the  water  motions.  To  pre- 
vent dissolution  of  these  calcium  carbonate  flakes,  I made  the  water 
slightly  basic  by  the  addition  of  0.2  kg  of  sodium  bicarbonate.  To  view 
the  flakes, I illuminated  them  with  a plane  of  light  from  a standard 

2 inch  by  2 inch  slide  projector.  I generated  a plane  of  light  by 


ohms 
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projecting  the  light  through  a slide  mount  which  had  two  razor  blades 
taped  to  the  mount  to  form  a slit  2 mm  wide.  The  light  entered 
the  tank  from  the  end  near  the  ice  slab  and  passed  through  the  bubble- 
free  ice  to  illuminate  the  convective  boundary  layer  next  to  th'e  ice. 

I viewed  and  photographed  the  flow  field  through  the  side  of  the  tank 
normal  to  the  light  plane.  I photographed  the  flow  using  a 35  mm  Nikon 
camera  with  a 50  um  lens  and  Kodak  Panatomic-X  film;  the  exposures  were 
usually  at  f-8  for  a period  of  30  s.  Plate  I is  an  example  of  a 
streak  photograph;  the  particle  motions  appear  on  the  photographs  as 
white  streaks  against' a black  background.  The  streak  photographs  give 
only  a qualitative  picture  of  the  flow  because  the  flakes  rotate  during 
the  exposure  and  thus  may  not  be  visible  during  the  entire  exposure. 

To  estimate  the  velocities  associated  with  the  convective  motions, 

1 injected  a neutrally  buoyant  mixture  of  India  ink  and  salt  water 
into  the  boundary  layer.  I then  measured  with  a stop  watch  the  time  it 
took  for  the  leading  edge  of  the  dye  to  pass  between  thermistors  in  the 
ice.  Also,  the  dispersion  of  the  dye  across  the  boundary  layer  gives 
an  estimate  of  the  boundary  layer  thickness.  I obtained  another  esti- 
mate of  the  boundary  layer  thickness  by  viewing  a Moire  grid  through  the 
boundary  layer  such  that  the  line  of  sight  was  parallel  to  the  ice 
surface.  The  large  change  in  the  index  of  refraction  between  the  water 
near  the  ice  and  the  ambient  water  distorts  the  Moire  grid  and  makes  the 
outer  limit  of  convective  motions  visible. 

To  maintain  homogeneous  far-field  conditions,  I relied  on  the  large 
mass  of  water  in  the  tank  and  a pumping  and  heating  system.  In  the 
pumping  system  shown  in  Figure  3 the  colder  less  saline  water 
generated  by  the  melting  ice  was  pumped  from  the  surface  of  the  tank 
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through  a constant  temperature  bath  and  returned  to  the  bottom  of  the  j, 

tank  behind  a perforated  diffusing  screen.  This  screen  allowed  the 

return  flow  to  mix  in  the  tank  without  affecting  the  convective  motions 

generated  by  the  melting  ice.  The  pumping  rate  was  of  order  100  mm3s_1 

which  produced  an  upward  velocity  in  the  tank  of  0.15  mm  s_1 , or  much 

less  than  the  10  mm  s'1  observed  velocities. 

To  summarize,  a typical  experiment  proceeded  in  the  following 
fashion:  I removed  the  ice  from  cold  storage  at  -15°C  and  allowed  the 
ice  to  warm  to  -1°C  to  prevent  thermal  fracture  during  placement  in  the 

I 

tank.  This  procedure  took  approximately  30  minutes  during  which  time  the 

. • 

tank  was  vigorously  stirred  and  allowed  to  come  to  rest  to  insure 

homogeneous  far-field  conditions.  I then  connected  the  thermistors  to  s 

the  data  acquisition  system  and  began  recording  just  prior  to  submerging 

• j. 

the  ice.  Once  the  ice  was  submerged,  a plumb  bob  assured  vertical 
alignment  of  the  ice  surface  and  the  experiment  began. 

♦ f 

During  the  experiment  I measured  the  lengths  of  the  melt  probes 
every  30  minutes  and  drew  far-field  salinity  samples.  I observed  and 
photographed  the  flow  using  reflecting  flakes  and  dye  injections. 

Towards  the  end  of  an  experiment,  photographs  recorded  the  final  ice 
shape.  The  experiment  ended  when  the  ice  melted  past  all  of  the  thermi- 
stors frozen  in  the  ice.  After  removing  the  array  and  any  residual  ice 
from  the  tank,  I rinsed  the  array  with  fresh  water  and  placed  it  in  a 


crushed  ice  and  water  bath  to  determine  the  0°C  thermistor  resistances. 
After  this  calibration,  data  recording  stopped  and  the  magnetic  tape 
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Introduction 

This  chapter  presents  and  discusses  the  nature  of  the  flow 
generated  by  a melting  vertical  ice  wall,  the  vertical  variation  of 
the  ice-water  interface  temperature,  and  the  ice  melt  rates.  To 
summarize,  the  laboratory  experiments  for  ambient  temperatures  less 
than  20"C  show  that  the  flow  next  to  the  bottom  portion  of  the  ice 
consists  of  a laminar  bidirectional  flow  with  a vertical  extent  limited 
by  a critical  Grashof  number  of  1.6  x 10M.  Above  this  laminar  region, 

I found  an  upward  rcr.ing  turbulent  boundary  layer.  Between  the  laminar 
and  turbulent  regions,  a horizontal  jet  of  ambient  water  flows  towards 
the  ice  that  melts  a notch  in  the  ice  wall.  Additional  experiments  show 
that  high  ambient  temperatures,  bubbles  in  the  ice,  and  tilting  the  ice 
wall  modify  the  flow,  temperature  and  melt  characteristics  from  those 
observed  in  the  bubble-free  case  for  T^  < 20°C. 

In  the  laboratory  study  most  of  the  experiments  concentrate  on  the 
case  of  vertically  oriented  bubble-free  ice  at  T^  < 20°C.  The  discus- 
sion continues  with  the  specific  cases  of  Td  > 20°C,  bubbly  ice  and 
ti lted  wal Is. 


Tho  (\jjtc-  of  Pubblc-Fz'co  Too,  T ^ < 20°C  ,md  Vertical  Wall. 

The  laminar  bidirectional  and  unidirectional  flow. 

At  the  bottom  portion  of  the  ice  I observe  laminar  flow,  either 
bidirectional  or  upward,  whose  characteristics  depend  on  the  far-field 
conditions.  When  Sro  and  Tt  lie  in  region  II  of  Figure  1,  the  region 
where  the  salinity  effect  on  density  opposes  the  temperature  effect, 
the  flow  is  upward  near  the  ice  and  downward  away  from  the  ice. 
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When  Su  and  Tu  lie  in  region  I of  Figure  1,  the  region  where  the 
salinity  and  temperature  effects  reinforce,  the  flow  is  everywhere 
upward. 

For  the  first  case,  the  velocity  field  consists  of  a thin  upward- 
flowing  region  next  to  the  ice  inside  of  a wider  downward  flowinq  region. 
Figure  6,  a sketch  of  the  observed  flow,  and  Plates  I and  II  illustrate 
the  bidirectional  nature.  These  two  plates  are  60  s streak  photographs 
of  the  flow  field  for  Su  = 33  °/oo  and  T„  = -1.50C;  f.ii(  iop  uf  Plate  II 
joins  Plate  I at  the  level  marked  C.  Immediately  next  to  the  ice 
Plates  I and  II  show  an  upward-flowing  region  of  order  1 mm  thick  which 
is  thinnest  at  the  bottom  of  the  ice  and  increases  in  thickness  with 
distance  from  the  ice  bottom;  dye  injections  show  the  velocities  in  this 
region  to  be  of  order  1 mm  s-1.  Two  features  on  each  plate  mark  the 
outer  edge  of  the  upward  flowing  layer:  the  apexes  of  the  "V"  shaped 
streak  trajectories,  and  the  region  of  short  streaks  and  dots  a few 
millimeters  from  the  ice  located  between  two  regions  of  longer  streaks, 
lhe  arrows  on  the  plates  point  to  several  of  the  "V"  shaped  trajectories 
and  the  arrow  shafts  lie  along  the  region  of  short  streaks  and  dots. 

The  trajectories  show  a parcel  of  water  moving  out  of  the  downward- 
flowing outer  region,  through  a region  of  no  vertical  motion,  and  into 
the  upward-flowing  inner  region. 

The  downward-flowing  outer  thermal  region  is  of  order  10  mm  thick 
and  grows  in  thickness  as  it  flows  downward  until  it  approaches  the 
bottom  ice  edge.  As  the  outer  flow  approaches  the  bottom  of  the  ice, 
the  flow  both  constricts  and  curves  slightly  under  the  ice.  Both  the 
flow  constriction  and  curvature  below  the  ice  bottom  results  from  the 
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TURBULENT  REGION 


Figure  6.  A sketch  of  the  observed  flow  field  next  to  the  ice. 
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Plate  I 


Bidirectional  flow  for  = -1.15°C  and  = 33.0  °/00 

upper  portion  of  the  Ice,  see  text  for  explanation, 

60  s exposure. 
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Plate  II.  Bidirectional  flow  next  to  lower  portion  of  the  ice, 
same  experiment  as  Plate  I,  60  s exposure. 
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absence  of  the  wall  which  causes  the  flow  to  accelerate;  this  accelera- 
tion, through  Bernoulli's  law,  reduces  the  pressure  within  the  boundary 
layer.  Hence,  the  flow  narrows  and  curves  under  the  ice  tip. 

The  large  difference  between  the  molecular  diffusivities  of  heat 
and  salt,  characterized  by  Le  = 200,  produces  the  observed  bidirectional 
flow.  Adjacent  to  the  ice,  the  melting  dilutes  and  cools  the  salt  water. 
As  a result  of  the  large  difference  in  diffusivities,  the  dilution  pro- 
duces positively  buoyant  water  near  the  ice  that  rises  while  the  cooling 
produces  negatively  buoyant  water  farther  from  the  ice  that  sinks.  Even 
though  the  dilution  buoyancy  near  the  wall  is  much  larger  than  the 
opposed  buoyancy  in  the  thermal  layer,  the  frictional  influence  of  the 
wall  is  largest  near  the  wall  so  that  the  velocities  in  each  region  are 
of  the  same  order. 

For  comparison  when  the  far-field  conditions  lie  between  the  freez- 
ing curve  and  the  maximum  density  curve  in  Figure  1,  I observe  an 
upward-flowing  unidirectional  laminar  boundary  layer.  In  this  case  both 
dilution  and  cooling  produce  less  dense  water  within  the  boundary  layer; 
cooling  decreases  the  density  because  the  sign  of  a is  positive  in 
Region  I of  Figure  1.  I performed  two  experiments  with  far-field 
conditions  T = 0.05°C,  S = 14.2  °/00,  and  T = 1.8°C,  S = 8.0  °/00  to 
determine  the  characteristics  of  the  flow  in  this  situation.  In  both 
cases  I observed  an  upward-flowing  unidirectional  boundary  layer  of 
order  10  mm  thick.  Plate  III,  a streak  photograph  from  the  first  experi- 
ment, displays  the  unidirectional  velocity  field  and  a comparison  with 
streak  Plates  I and  II  showing  the  bidirectional  flow  clearly  illustrates 
the  difference  between  the  two  cases. 
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Plate  in. 


Unidirectional  upward  laminar  flow,  T = -0.4°C  S = i«  0/ 

OO  * oc-  10  ' 00  » 

60  s exposure. 
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Transition  to  turbulence. 

The  turbulence  studies  concentrate  on  the  case  with  opposing 
buoyancy  forces  because  common  oceanic  salinities  and  temperatures  lie 
in  region  II  of  Figure  1,  or  the  region  of  opposing  buoyancy  forces. 

For  the  bidirectional  case,  the  onset  of  turbulence  in  the  inner  saline 
layer  determines  the  vertical  extent  of  the  laminar  bidirectional  flow. 

I observe  the  transition  from  laminar  to  turbulent  flow  in  two  ways;  by 
viewing  a Moire  grid  through  the  boundary  layer  and  by  injecting  dye 
into  the  boundary  layer. 

The  use  of  both  these  techniques  reveals  small  amplitude  waves  with 
horizontal  crests  parallel  to  the  ice  surface  and  perpendicular  to 
gravity  at  the  outer  edge  of  the  saline  layer  in  the  upper  third  of  the 
laminar  region.  These  waves  propagate  upward  with  a phase  speed  of 
order  1 mm  s'1  and  a wavelength  of  order  5 mm.  The  wave  amplitude  grows 
to  approximately  3 nm  at  which  point  the  waves  "break";  about  30  mm 
above  the  point  at  which  the  waves  break,  the  inner  flow  becomes  turbu- 
lent. When  compared  to  the  theoretical  stability  study  of  Heiber  and 
Gebhart  (1971),  these  laboratory  observations  show  that  the  stability  of 
the  inner  saline  layer  is  governed  by  a coupled  buoyancy-momentum 
mechanism. 

Vertical  grooves  which  form  in  the  ice  in  the  upper  third  of  the 
laminar  region  give  indirect  evidence  of  an  additional  transitional 
mechanism.  These  grooves  have  a width  and  depth  of  order  1 rrm  and  a 
separation  of  order  25  mm.  Gebhart  (personal  communication)  attributes 
these  grooves  to  a secondary  flow  feature  consisting  of  longitudinal 
rolls  that  flow  in  a circular  fashion  perpendicular  to  the  direction  of 
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the  mean  flow.  The  increased  melting  results  from  the  advection  of  heat 
towards  the  ice  by  these  rolls  in  the  flow. 

When  the  inner  flow  becomes  turbulent,  the  turbulence  diffuses  the 
dilute  water  near  the  ice  further  outward  where  the  dilution  buoyancy 
overwhelms  the  thermal  buoyancy  to  produce  an  upward  unidirectional  flow. 
At  the  level  of  transition  to  turbulence,  the  divergence  between  the 
upward-flowing  turbulent  boundary  layer  and  the  downward-flowing  laminar 
thermal  boundary  layer  produces  a horizontal  jet  of  ambient  water 
flowing  towards  the  ice.  Plates  IV  and  V show  the  jet  which  appears  as 
a broad  band  of  horizontal  streaks.  The  point  labeled  "B"  at  the  ice- 
water  interface  on  each  photograph  marks  the  stagnation  stream  line  of 
the  jet  and  will  be  called  the  bifurcation  point.  Water  in  the  jet 
above  the  bifurcation  point  enters  the  upward-flowing  boundary  layer 
while  water  below  this  point  enters  the  downward-flowing  thermal 
boundary  layer.  Because  the  streaks  in  the  jet  are  much  longer  than 
the  horizontal  streaks  outside  the  jet,  the  velocities  in  the  jet  are 
much  greater  than  the  suction  velocity  associated  with  the  boundary 
layer.  The  speeds  in  the  jet  are  of  order  1 mm  s_1  and  the  jet  has  a 
thickness  of  approximately  100  mm. 

As  a -esult  of  the  jet,  the  ice  melts  fastest  in  the  region  of  the 
jet  and  a notch  forms  in  the  ice.  Plate  VI  shows  the  ice  after  being 
immersed  in  water  with  a salinity  of  29  °/00  and  a temperature  of  +1.5°C, 
Td  = 3.2°C,  for  a period  of  5 hr.  17  min.  Point  B is  the  bifurcation 
point  taken  from  Plate  IV  which  shows  the  initial  ice  shape  for  this 
experiment.  The  experiment  ended  when  the  notch  melted  through  the  ice. 
Once  the  notch  begins  to  form,  it  grows  rapidly  and  stays  at  the  same 
level . 


r 


Plate  VI 
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Increased  melting  at  the  bifurcation  point,  same  experiment 
as  Plate  IV,  after  melting  for  5 hr  17  min. 
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The  length  of  the  laminar  regime,  l7,  is  the  distance  from  the 
bottom  tip  of  the  ice  to  the  bifurcation  point.  Table  2 gives  the 

1 

observed  values  of  L^,  the  average  wall  salinity  and  temperature,  and  L 

the  saline  and  thermal  Grashof  numbers  for  each  of  the  laminar  experi- 
ments. With  increased  thermal  driving,  decreases  from  300  mm  at 
Td  = 1.65°C  to  110  mm  at  T^  = 8.71°C.  The  average  is  1.6  x 108, 
while  the  average  is  3.8  x 106.  The  coefficient  of  saline  expansion 
used  is  8.1  x 10"4  °/0o”1  from  Appendix  I,  while  the  coefficient  of 

1 1 

thermal  expansion  used  to  calculate  G-j.  is  evaluated  at  the  far-field 

conditions.  Because  Gr^  is  two  orders  of  magnitude  larger  than  Gr^,  the 

stability  of  the  saline  layer  determines  L7 , which  is  a convenient  i 

£ t | 

length  scale  for  the  presentation  of  the  laminar  wall  tei  ture  data. 


Laminar  wall  temperatures. 


In  the  bidirectional  case  the  average  wall  temperature,  T , 

w 

decreases  as  T.  increases  with  the  following  vertical  variation  of  T : 

d 3 w 

T is  coldest  at  the  bottom  of  the  ice  and  slowly  warms  with  increasing 

W 

distance  from  the  ice  tip.  Figure  7 exhibits  the  vertical  variation 


of  Tw  for  the  laminar  experiments.  The  normalized  interface  temperature 


r is 


Figure  7 shows  r as  a function  of  the  distance  from  the  bottom  tip 
of  the  ice  divided  by  the  length  of  the  laminar  regime,  L^.  Where  the 
L^  was  not  measured  directly,  I computed  it  using  a critical  Gr^  = 

1.6  x 10®  and  the  measured  salinity  difference  between  the  wall  and  the 
far  field.  The  figure  shows  that  r varies  by  as  much  as  20%  over  the 
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Table  2 


Experiment 

number 

T. 

(°c) 

<°/oo> 

Td 

(°C) 

l 

(mm) 

Tw 

(°c) 

1 

-0.5 

29.0 

1.08 

-1.53 

2 

-0.72 

33.7 

1.13 

260 

-1.81 

3 

0.0 

29.0 

1.65 

300 

-1.49 

4 

1.55 

29.0 

3.13 

270 

-1.11 

5 

2.20 

30.0 

3.83 

180 

-0.97 

6 

2.70 

29.0 

4.28 

-0.80 

7 

2.66 

34.4 

4.55 

155 

-1.03 

8 

6.85 

34.0 

8.71 

110 

-0.32 

Suninary  of  laminar 

bidirectional  experiments. 

l not  observed 

for  experiments  1 

and  6. 

Table  3 

T. 

(°c) 

(°/oo) 

Td 

(°C) 

Tw 

(°C) 

Sw 

(°/oo) 

r 

+0.05 

14.2 

0.78 

-0.66 

12.8 

0.90 

+1.80 

8.0 

2.50 

-0.29 

5.6 

0.70 

Far-field  conditions  for  the  two  unidirectional  experiments 
and  the  resulting  interface  conditions. 
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length  of  the  laminar  layer,  which  corresponds  to  a temperature  change 
of  approximately  0.35°C  for  the  experimental  conditions.  The  maximum 
variation  occurs  in  the  range  of  Td  from  3°C  to  5°C.  For  > 5°C  the 
variation  of  r decreases  because  T at  any  level  begins  to  approach  the 

W 

fresh  water  1 i mi t . 

Figure  8 gives  T , obtained  by  averaging  all  of  the  temperature 
w 

measurements  in  the  laminar  region  from  an  experiment  and  then  normali- 
zing with  Tf-p(Soo)  for  that  experiment  as  a function  of  Td-  The  norma- 
ized  average  interface  temperature,  r,  where 


gives  a measure  of  the  dilution  occurring  at  the  interface,  r varies 

from  1 , when  the  ice  is  in  equilibrium  with  the  far-field  and  no 

melting  occurs,  to  0 when  the  melt  rate  is  great  enough  to  prevent 

salt  from  reaching  the  interface.  The  interface  temperature  and  hence 

the  interface  salinity  both  approach  the  fresh  water  limit,  T = 0°C  and 

w 

S = 0 °/00,  as  the  thermal  driving  increases  with  90%  of  the  change 
occurring  for  Td  1 10°C.  For  Td  > 10°C  the  interface  conditions 
asymptotically  approach  the  fresh  water  limit. 

For  the  unidirectional  laminar  case,  I found  vertically  uniform 


wall  conditions  in  each  of  the  experiments.  Table  3 gives  the  far-fieli 
conditions,  TH  and  T for  the  two  experiments.  The  uniform  S and  T 

u w W W 

and  unidirectional  flow  suggests  that  the  analysis  of  Gebhart  and  Pera 
(1971)  is  applicable  to  this  case;  however,  they  do  not  include  the 
effect  of  blowing  at  the  interface. 
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Observations  of  melting  in  the  laminar  bidirectional  region. 

In  the  laminar  bidirectional  case,  the  ice  ablation  is  greatest  at 
the  ice  bottom,  then  decreases  with  distance  from  the  bottom  until  the 
transition  region  is  reached  where  the  ablation  abruptly  increases. 

Plate  VII  is  a profile  of  the  ice  from  the  same  experiment  as  Plates  IV 
and  VI  taken  220  min  after  the  start  of  the  experiment.  The  original 
ice  thickness  was  55  mm;  the  dashed  line  on  Plate  VII  represents 
the  approximate  location  of  the  ice-water  interface  at  the  beginning  of 
the  experiment.  The  large  melt  rate  at  the  ice  bottom,  which  rapidly 
changes  the  ice  shape  in  this  region,  results  from  the  influx  of  ambient 
water  into  the  upward-flowing  saline  boundary  layer  from  below  the  ice. 
The  melt  rate  then  decreases  as  the  saline  boundary  layer  thickness 
increases  and  the  thermal  boundary  layer  thickness  decreases.  At  the 
level  where  the  inner  saline  flow  becomes  turbulent,  the  horizontal  jet 
of  ambient  water  flowing  towards  the  ice  increases  the  melt  rate  to 
produce  the  observed  notch.  Above  the  notch,  in  the  turbulent  regime, 
the  melt  rate  is  less  than  the  melt  rate  in  the  notch,  but  greater  than 
the  melt  rate  in  the  laminar  regime.  Figure  9 for  the  same  experi- 
ment as  Plate  VII  gives  the  ablation  rates  measured  with  the  transit; 
the  melt  rates  are  of  order  100  mm  day-1. 

For  all  experiments  with  bubble-free  ice  for  Td  < 20°C,  uniform  ice 
ablation  takes  place  at  each  horizontal  level  and  during  an  experiment 
the  vertical  edges  of  the  ice  become  rounded.  However,  the  edge  radius 
of  curvature  is  less  than  10  mm,  or  much  smaller  than  the  cross  ice 
dimension  of  200  mm;  therefore,  edge  rounding  also  has  little  effect  on 
the  overall  flow  and  the  observed  melt  rates. 
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Figure  10  gives  the  melt  rate  measured  at  the  vertical  and 
horizontal  mid-point  of  the  laminar  layer  for  the  experiments  from  the 
transit  measurements.  A power  law  regression  of  the  transit  data  gives 

M - 0.34  y-»  - U . (24) 

Figures  11  and  12  display  the  melting  data  in  the  conventional 
manner;  in  terms  of  the  Nusselt  number,  Nu,  and  Sherwood  number,  Sh,  as 
a function  of  the  local  Grashof  number.  In  Figure  11  the  distance  from 
the  ice  bottom  to  the  point  of  measurement  was  used  to  compute  Gr<-,  Nu 
and  Sh.  In  Figure  12  the  distance  from  the  bifurcation  point  to  the 
point  of  measurement  was  used  to  compute  Gr^,  Nu  and  Sh.  Also,  I evalu- 
ated at  the  far-field  conditions  to  give  temperature  anomalies 
maximum  weight  in  the  calculation  of  Gr^. 

Because  of  poor  estimates  of  the  ice  thickness  between  a thermistor 
and  the  outer  edge  of  the  ice,  the  melt  rates  obtained  from  the  melt- 
through  time  are  only  approximate  and  subject  to  large  errors. 
Consideration  of  the  transit  data  alone  reduces  the  scatter  in  the  data. 
The  resulting  correlation  of  Nu  on  Gr<-  and  Grj  indicates  that  a 1/4 
power  law  dependence  is  not  appropriate;  for  correlations  on  Gr^,  the 
data  indicate  a 1/8  power  law,  while  for  correlations  on  Gr^,  the  data 
indicate  a 1/2  power  law.  The  failure  of  a 1/4  power  law  for  Sh  results 
from  the  omission  of  blowing  from  previous  studies  that  give  the  1/4 
power  law  correlation.  LaPadula  and  Mueller  (1970)  show  that  the  effect 
of  blowing  increases  with  increasing  Pr  or  Sc  which  explains  the  observed 
correlation  differences  between  Nu  and  Sh.  The  effect  of  blowing  on  the 
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Figure  11.  Correlation  of  Nu  and  Sh  on  Gr$;  • corresponds  to  Nu, 
■ corresponds  to  Sh.  The  dashed  lines  each  have  a 
slope  of  1/4. 


i 

i 


66 


heat  transfer  is  small  because  of  the  low  value  of  Pr  while  the  effect 
of  blowing  on  the  salt  transfer  is  large  because  Sc  is  large,  200  times 
greater  than  Pr. 

T he  turbulent  boundary  layer . 

Above  the  bidirectional  laminar  flow,  1 observe  an  upward-flowing 
turbulent  boundary  layer.  In  this  boundary  layer  the  turbulent  motions 
disperse  the  dilute  upward-flowing  water  adjacent  to  the  ice  across  the 
thickness  of  the  boundary  layer,  with  the  result  that  the  upward 
buoyancy  and  momentum  overwhelm  the  downward  thermal  buoyancy.  The 
large  difference  in  the  index-of-refract ion  between  the  low  salinity 
water  near  the  wall  and  the  higher  ambient  salinity  water  allows  for 
direct  observation  of  the  turbulent  flow.  When  I sighted  across  the 
boundary  layer  at  a Moire  grid,  the  low  salinity  water  appeared  as  1 nm 
wide  wisps  that  were  randomly  ejected  from  the  flow  near  the  wall.  The 
turbulence  then  rapidly  mixes  the  low  salinity  water  to  the  point  where 
it  is  no  longer  visible. 

Dye  injections  show  that  the  outer  edge  of  the  boundary  layer 
fluctuates  toward  and  away  from  the  ice  with  the  passage  of  turbulent 
eddies  and  these  eddies  entrain  ambient  fluid  into  the  boundary  layer. 
Townsend  (1976)  shows  that  the  boundary  between  a turbulent  flow  and 
the  quiescent  fluid  remains  sharp  while  wavering  back  and  forth.  The 
temperature  records  give  further  evidence  of  this  phenomena.  Figure  13 
is  an  example  of  the  temperature  records  from  the  turbulent  region  and 
these  records  are  discussed  fully  below.  The  boundary  layer  thickness 
is  approximately  10-30  nin  and  becomes  thicker  with  increasing  downstream 


TEMPERATURE  (#C) 
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distance.  Dye  injections  show  the  maximum  upward  velocity  to  be  of 
order  10  mu  s'1  located  at  2-3  mm  from  the  ice.  The  upward  velocities 
increase  with  both  increasing  downstream  distance  and  thermal  driving. 

Ice-water  interface  temperatures  in  the  turbulent  region. 

In  the  turbulent  flow,  I find  not  only  unidirectional  upward  flow, 

but  also  uniform  ice-water  interface  temperatures  over  the  entire  ice 

surface.  The  length  of  the  turbulent  flow  regime  varies  from  0.3  m to 

1.0  m and  in  all  cases  the  measured  values  of  T are  within  0.02°C  of 

w 

each  other  for  any  particular  experiment.  Table  4 summarizes  the  far- 
field  conditions  and  resulting  interface  temperatures  and  melt  rates 
for  the  turbulence  experiments.  Figure  14  gives  the  observed  values  of 
r for  0°C  £ Trf  1 20°C.  This  figure  shows  that  r decreases  linearly  with 
increasing  T^  for  T^  * 7°C  and  then  r approaches  the  fresh  water  limit 
of  0.0  as  T^  increases  to  20°C.  The  temperature  line  on  Figure  14  is 
the  curve 

r = 1.0  - tanh(0. 155  x Trf)  . (25) 

This  function  is  especially  suited  to  represent  the  data  because  they 
both  linearly  decrease  for  small  arguments  and  approach  0.0  for  large 
arguments. 

The  temperature  records  yield  not  only  the  interface  temperature, 

but  also  provide  a qualitative  description  of  the  cross  stream  turbulent 

temperature  fluctuations.  Figure  13  shows  the  temperature  as  measured 

by  a thermistor  moving  across  the  boundary  layer  at  a speed  equal  to  the 

ablation  rate.  In  this  case  T = -0.64°C,  T = +2.7°C,  Tj  * 4.4°C  and 

w ® d 
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Table  4 

T„ 

s 

oo 

Td 

Tw 

x M x 10" 3 mm  s"1 

(°C) 

O 

o 

o 

(°C) 

(°C) 

(mm) 
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the  melting  rate,  M,  is  7.3  mm  hr-1.  This  record  shows  several  features 
of  the  turbulent  flow  as  the  thermistor  traverses  the  boundary  layer. 
First,  within  2 mm  of  the  ice,  no  significant  temperature  fluctuations 
occur  and  the  temperature  increases  monotonical ly.  Second,  in  a region 
from  2 mm  to  10  mm  away  from  the  ice,  the  temperature  fluctuations  occur 
as  large  as  0.5°C  and  the  magnitude  of  the  fluctuations  increase  with 
increasing  distance  from  the  ice.  Third,  the  large  temperature  fluctu- 
ations at  the  outer  edge  of  the  boundary  layer  result  from  the  wavering 
motion  of  the  boundary  between  turbulent  and  quiescent  fluid.  Last,  the 
thermal  boundary  layer  thickness  of  10  mm  roughly  equals  the  momentum 
boundary  layer  thickness  observed  by  dye  injections.  The  sampling  rate 
of  once-per-minute  provides  only  a qualitative  picture  because  fluctu- 
ations with  frequencies  greater  than  0.02  s-1  remain  undetected  and  both 
dye  injections  and  the  Moire  grid  show  fluctuations  in  the  velocity 
field  at  frequencies  of  at  least  1 s-1. 


Melt  rates  in  the  turbulent  region. 

For  the  ice  adjacent  to  the  turbulent  flow,  I found  that  the  melt 
rate  dependence  on  T^  changes  from  a power  law  dependence  to  a linear 
dependence  at  Td  = 9°C.  Also,  at  a particular  Td,  I found  the  melt  rate 
to  have  the  following  vertical  variation:  The  melt  rate  is  highest  at 
the  bifurcation  point  and  decreases  with  increasing  distance  from  the 
bifurcation  point.  For  all  experiments  at  Td  < 20°C  the  ice  melted  in 
a smooth  fashion;  however,  the  wakes  from  the  thermistors  and  melt 
probes  produce  1 mm  wide,  1 mm  deep  and  50  mm  long  perturbations  on  the 
ice  surface. 


i 

i 


y 


f 


72 


Figure  14  gives  the  observed  melt  rates,  M^,  for  up  to  30°C, 

where  M is  the  melt  rate  measured  by  transit  at  a distance  x above  the 
x 

bifurcation  point  multiplied  by  x1/4.  For  all  cases,  x lies  between 
0.2  m and  0.5  m.  I chose  this  formulation  for  because  the  uniform 
wall  temperatures  strongly  suggest  an  x1/4  similarity  transformation 
where  the  melt  rate  decreases; as  x-1/4.  Therefore,  represents  a 
characteristic  melt  rate  for  a given  Td  that  is  independent  of  x;  this 
convenient  formulation  aids  in  the  similarity  development  discussed  in 
the  Similarity  section. 

Figure  15  gives  the  average  observed  melt  rate  M without  the  x1//m 
dependence  and  M is  the  average  of  all  transit  r lsured  melt  rates  taken 
during  each  experiment.  The  solid  line  on  Figure  15  is  a power  law 
regression  fit  of  all  the  data  given  by 

M = 2.85  x 10‘4T  i-3*  mm  s-1r  2 . (26) 

a 

However,  as  seen  for  a linear  dependence  best  fits  the  data  for 

Td  - rc- 

Returning  to  Figure  14,  the  data  show  that  the  dependence  of 
on  T^  changes  in  the  vicinity  of  9°C.  For  Td  < 9°C  a power  law 
regression  gives 

M = 7.54  T j1  • e 3 x 10-4  mm^4  s“'  (27) 

x d 

while  for  9°C  < Td  < 20°C  the  data  best  fit  a linear  regression  given  by 
Mx  = [38.06  Td  - 7.29]  x 10*4  mm3/4  s’1  . (28) 
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Neither  the  linear  nor  the  power  law  fit  of  the  data  is  satisfactory 

over  the  entire  range  of  T^.  The  change  in  melt  rate  dependence  results 

from  the  change  in  buoyant  forcing  determined  by  the  difference  between 

the  interface  and  far-field  conditions.  For  T.  < 9°C,  both  (S  - S,  ) and 

(Ta,  - Tw)  increase  linearly  with  increasing  T^;  but  for  > 9°C,  the 

saline  difference  remains  nearly  constant  while  the  thermal  difference 

continues  to  increase.  Because  salinity  determines  the  density,  the 

buoyant  forcing  increases  with  for  < 9°C  and  the  increasing 

buoyant  forcing  and  temperature  difference  across  the  boundary  layer 

combine  to  give  a melt  rate  that  increases  as  some  power  ’of  T^.  Above 

Td  * 9°C  the  buoyant  forcing  remains  constant  while  the  temperature 

difference  increases  with  T^,  therefore  the  melt  rate  increases  linearly 

with  T..  For  T,  > 20°C  the  thermal  effect  overwhelms  the  saline  effect 
d d 

on  density  and  produces  the  flow  reversal  discussed  in  the  next  section. 
In  the  linear  regression  I omitted  the  data  points  at  = 27 . 7°C  and 
26.6°C  because  the  flow  in  these  experiments  had  reversed  and  as  such  no 
longer  represents  the  flow  conditions  found  at  < 20°C. 
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The  High  Temperature  Flow  Reversal 

Experiments  at  Td  > 25°C  show  that  the  flow  next  to  the  ice 
consists  of  a laminar  bidirectional  boundary  layer  at  the  tO£  of  the 
ice  and  a downward  flowing  turbulent  boundary  layer  beneath  the  laminar 
flow.  In  addition,  the  ice  ablates  smoothly  in  the  laminar  region  and 
irregularly  in  the  turbulent  region  where  the  ablation  produces  cusps  in 
the  ice  surface. 

I performed  two  experiments  at  > 25°C,  Td  = 27.7°C,  and  26.6°C 
and  Figure  16,  a sketch  of  the  observed  flow,  illustrates  the  fluid 
motions  and  the  resulting  ice  shape.  At  the  top  of  the  ice,  I find  a 
laminar  bidirectional  layer,  60  mm  long,  with  the  inner  layer  flowing 
up  and  the  outer  layer  flowing  down.  The  inner  layer  is  1-2  mm  thick 
while  the  outer  layer  is  10-20  mm  thick  and  the  ice  adjacent  to  this 
flow  melted  smoothly.  Below  this  laminar  flow  1 observe  a 20-30  mm 
thick  very  turbulent  downward  flowing  boundary  layer  that  grows  in 
thickness  with  distance  down  the  ice.  The  ice  next  to  this  layer  melts 
irregularly  to  produce  cusps  on  the  ice  surface.  Plate  VIII  shows  the 
smooth  region  of  the  top  of  the  ice  and  the  cusped  region 
over  the  rest  of  the  ice.  At  the  beginning  of  the  experiment,  this 
piece  of  ice  was  planed  smooth;  I took  the  photograph  after  the  ice  had 
been  submerged  in  salt  water  at  24.8°C  and  33.6  °/0o»  Td  = 26.6CC,  for 
23  minutes. 

The  cusps  have  a downstream  length  scale  that  becomes  longer  as  the 
experiment  progresses.  Plates  IX  and  X taken  10  min  and  23  min,  respec- 
tively, after  the  start  of  the  experiment,  show  the  evolution  of  the 
cusps.  Initially,  well-ordered  cusps  form  within  12  min  from  the  begin- 
ning of  an  experiment  with  a downstream  length  of  10  im  and  a cross 
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Plate  VIII.  Ice  shape  after  melting  for  23  min,  = 24.8  C, 
S =33.6  °/00-  Note  smooth  upper  portion  and 

^oo 

cusped  lower  portion. 


Plate  IX.  Cusps  in  ice  surface  formed  for  Tt  = 26.0°C,  S =31.0  °/00, 
after  melting  for  10  min. 
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min,  same  experiment  as 
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stream  length  of  20  mm.  After  30  min  the  cusps  are  disorganized  with 
equal  downstream  and  cross  stream  lengths  of  approximately  30  mm.  In 
both  cases  the  cusps  are  larger  where  the  boundary  layer  is  thicker. 

The  depth  of  the  cusps  increases  with  time  and  the  final  depth  of  the 
cusps  is  approximately  5 mm.  Dye  observations  showed  that  the  physical 
effect  of  the  cusps  is  to  create  a back  eddy  in  the  flow  behind  each  ice 
crest,  as  illustrated  in  Figure  16  by  the  small  arrows,  where  the 
buoyant  water  produced  at  the  ice  surface  reinforces  the  upward  flow 
near  the  ice.  When  dye  was  injected  into  a cusp,  the  dye  flowed  upward 
and  then  outward  at  the  crest  of  the  cusp,  then  downward  in  the  mean 
flow.  Occasionally,  some  dye  flowed  over  a crest  and  into  a higher  cusp 
before  being  swept  downward.  I observed  no  net  upward  flow  adjacent  to 
the  ice  in  the  turbulent  portion  of  the  boundary  layer. 

I estimated  the  melting  rates  in  these  experiments  by  measuring  the 

ice  thickness  before  and  after  the  experiment.  Because  the  ice  surface 

is  cusped  and  the  ice  thickness  at  the  end  of  the  experiment  is  measured 
to  the  top  of  the  cusps,  this  technique  gives  a lower  bound  to  the  ice 
ablation.  The  melt  rate  for  Tw  = 24.8°C  and  Sv  = 33.6  °/0o*  ^ = 26.6V'C, 
was  175  mm  hr-1;  the  melt  rate  for  Tv  = 26.0°C  and  Sv  = 31.0  °/00, 

Tj  - 27.7°C,  was  114  mm  hr"1.  The  uncertainty  in  measuring  the  final 

ice  thickness  produces  the  disparity  in  melt  rates.  The  measured 
interface  temperatures  for  the  second  experiment  were  all  -.02°C  so 
that  fresh  water  is  produced  at  the  wall. 

The  separation  of  salinity  and  temperature  effect  on  density  across 
the  boundary  layer  produces  the  high  temperature  reversal.  Next  to  the 
ice  there  is  a small  amount  of  very  positively  buoyant  water  inside  of 
a large  amount  of  slightly  less  buoyant  water  that  at  T^  * 20°C 
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overwhelms  the  upward  buoyancy.  The  high  temperature  reversal  chapter 
gives  a complete  discussion  of  this  phenomena. 

Bubble  Ice  Ejq  cvinenis 

Paterson  (1975)  shows  that  glacial  ice,  the  source  of  icebergs, 
has  an  air  content  as  high  as  7"  by  volume.  The  formation  of  glacial 
ice,  which  results  from  the  metamorphosis  of  snow  under  the  weight  of 
additional  snow,  traps  air  in  the  ice.  To  investigate  the  effect  of 
released  air  from  the  ice  on  the  convective  processes,  I performed  two 
experiments  with  ice  that  had  a CO;  content  of  1.5*  and  51.  I used  CO; 
because  of  its  high  solubility  in  cold  water,  a property  that  enabled  me 
to  easily  produce  bubble  ice. 

To  produce  ice  with  a suitable  bubble  content,  I let  1 kg  of  dry 
ice  sublimate  in  water  at  0°C  in  the  bubble-free  ice  tank,  then  lowered 
the  room  temperature  to  -20°C.  Because  the  sides  of  the  tank  are 
insulated,  the  freezing  takes  place  at  the  top  of  the  tank  and  progresses 
downward.  As  the  ice  grows,  the  dissolved  CO;  is  excluded  from  the 
crystalline  ice  structure  and  forms  bubbles  at  the  ice-water  interface. 
These  bubbles  then  freeze  into  the  ice  block  as  the  ice  continues  to 
form.  The  bubbles  have  a characteristic  diameter  of  1 mm  and  a uniform 
distribution  within  the  ice.  I determine  the  CO;  content  of  the  ice  by 
melting  a known  mass  of  ice  in  heated  paraffin  oil  and  collecting  the 
CO;  in  a burrette. 

I chose  far-field  conditions  of  T^  = +2.6°C,  Sa,  = 32.7  °/00  for  the 
experiments  so  that  I could  compare  the  results  to  b"bble-free  ice 
experiments  performed  with  similar  far-field  conditions.  For  the  experi- 
ment with  5*  CO;  at  T.  = 4.06,  the  turbulent  interface  temperature  is 
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-1.35°C  which  corresponds  to  a value  of  r = 0.75.  If  the  ice  were 
bubble-free,  r would  equal  0.44  which  corresponds  to  a warmer  interface 
temperature.  Interface  temperatures  from  the  other  experiment  are  not 
available  due  to  the  failure  of  the  data-logger. 

The  melting  rate  estimates  for  both  experiments  fall  within  10'T 
of  the  values  of  the  bubble-free  experiments;  however,  the  ice  ablation 
does  not  proceed  uniformly  over  the  ice  surface.  Plate  XI,  from  the  5", 
CO;  case,  shows  the  ice  form  after  melting  for  3 hours;  cusps,  similar 
to  those  found  in  the  high  temperature  bubble-free  experiments,  form 
on  the  ice  surface.  These  cusps  have  a length  of  30  mm  and  a depth  of 
5mn.  Because  of  the  irregular  nature  cf  the  ablation,  I could  only 
estimate  the  melt  rate.  The  melting  ice  continuously  releases  CO; 
bubbles  into  the  boundary  layer  where  these  bubbles  quickly  rise.  The 
bubbles  always  remain  within  5 mm  of  the  ice  and  as  they  rise  they  drag 
water  with  them,  imparting  upward  momentum  to  the  water.  Dye  injections 
show  the  upward  velocity  30  cm  from  the  tip  to  be  8 mm  s'1.  Compared 
to  the  bubble-free  experiments,  these  experiments  give  a colder  inter- 
face temperature  and  show  that  the  ice  no  longer  melts  in  a smooth 
fashion.  Also,  many  features  remain  the  same  for  both  experiments;  at 
the  bottom  of  the  ice  the  flow  is  bidirectional  and  when  the  inner  saline 
layer  becomes  turbulent,  the  flow  becomes  unidirectional  but  the  charac- 
teristic notch  at  the  transition  point  does  not  form.  I also  observe 
uniform  wall  temperatures  in  the  turbulent  regime  and  the  upward 
velocities  in  the  turbulent  boundary  layer  are  approximately  10  mm  s'1. 

The  effect  of  the  bubbles  appears  to  be  threefold.  First,  the 
bubbles  in  the  ice  roughen  the  ice  surface  which  previously  was  smooth 
which  enhances  the  salt  transport  to  the  ice  and  lowers  Tw  from  the 
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Plate  XI.  Bubble  ice  experiment,  55.  CO.  by  volume,,  ice  shape  after 
melting  for  2 hr  54  min;  Tv  = 2.26°C,  Sv  =32.8  °/00. 
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corresponding  bubble-free  case.  Second,  the  ice  melts  in  an  irregular 
fashion  that  produces  cusps  in  the  ice.  Third,  the  bubble  release 
increases  the  production  of  positive  momentum  in  the  boundary  layer. 

The  Tilted  Wall  Experiment 

All  ice  walls  in  nature  are  not  vertical;  for  example,  pressure 
ridge  keels  are  composed  of  randomly  oriented  ice  blocks  and  icebergs 
have  slanted  walls.  Usually,  the  seawater  is  underneath  the  ice  but 
occasionally  this  configuration  is  reversed  as  in  the  case  of  iceberg 
rams.  To  determine  the  effect  of  tilting  the  ice  wall  on  the  convective 
motion  generated  by  the  melting  ice,  I performed  an  experiment  where  the 
angle  of  the  ice  with  respect  to  the  vertical  varied  from  0°  to  60°.  In 
this  experiment  I used  a slab  of  bubble-free  ice  that  measured  0.11  m 
thick,  0.6  m long,  and  0.2  m wide  which  had  the  thermistor  array  frozen 
in  the  ice  such  that  the  thermistors  were  18  mm  from  one  ice  face. 
Initially,  I placed  the  ice  vertically  in  the  tank,  then  tilted  the  ice 
at  15°,  32°,  and  60°  to  the  vertical.  The  ice  spent  40  min  at  0°,  26 
min  at  15°,  34  min  at  32°,  and  30  min  at  60°;  the  side  of  the  ice 
nearest  the  thermistors  faced  downward.  The  far-field  conditions  for 
this  experiment  are  T^  = 3.7°C  and  = 29.5  °/00  which  gives  T^  = 5.31°C. 

The  flow  visualization  methods  are  the  same  as  those  used  in  the 
vertical  wall  experiments;  namely,  dye  injections  and  streak  photographs. 
For  0°  tilt,  these  techniques  show  that  the  flow  is  the  same  on  each 
side  of  the  ice  and  as  before  consists  of  a bidirectional  laminar  layer 
at  the  bottom  of  the  ice  and  an  upward-flowing  turbulent  boundary  layer 
over  the  rest  of  the  ice. 
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As  the  tilt  angle  increases  the  inner  saline  flow  under  the  ice 
becomes  thinner  and  less  energetic.  The  outer  laminar  thermal  flow 
under  the  ice  becomes  unsteady  in  that  the  colder  more  dense  water 
forms  plumes  that  sink  into  the  interior  of  the  tank  and  the  upward 
flowing  turbulent  boundary  layer  disappears.  The  disappearance  of  this 
flow  results  from  the  density  stratification  under  the  ice  that 
stabilizes  the  laminar  saline  flow.  For  60°  tilt  the  outer  flow  takes 
on  the  appearance  of  convection  from  a horizontal  plate. 

For  the  water  above  the  ice,  the  convection  above  the  ice  is  more 
energetic  than  that  below  the  ice  because  the  dilute  melt  water  is  now 
free  to  rise.  For  the  three  tilt  angles,  I found  no  laminar  flow  along 
the  ice  surface,  only  convection  that  resembles  convection  from  a 
horizontal  plate. 

Plate  XII,  a 15  s streak  photograph  at  a 60°  tilt,  depicts  the 
difference  between  the  convection  on  each  side  of  the  ice.  The  motions 
in  the  upper  portion  of  the  picture  are  disorganized  and  severely 
influenced  by  the  proximity  of  the  surface  which  is  the  white  line  at 
the  top  of  the  picture.  Below  the  ice,  the  fluid  motion  is  less 
dramatic;  however,  several  sinking  plumes  are  shown. 

Plate  XII  also  shows  the  proximity  of  the  thermistors  to  the 
lower  ice  face  which  enabled  me  to  determine  the  lower  interface 
temperature  for  the  different  angles.  Because  the  small  amount  of  ice 
between  the  thermistor  and  the  interface  has  an  e-folding  thermal 
response  time  of  7 min  and  each  portion  of  the  experiment  lasts  several 
response  times,  the  thermistors  accurately  measure  T provided  T is 
the  same  on  each  side  of  the  ice.  I computed  an  accuracy  of  ± 0.07°C 
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for  these  measurements  under  the  assumption  of  equal  wall  temperatures 
by  assuming  the  worst  case.  For  this  case,  T = 0°C  at  the  upper  ice 

W 

surface  and  a linear  temperature  gradient  exists  through  the  ice. 

I chose  0°C  at  the  upper  side  as  a worst  case  because  the  wall  tempera- 
ture is  usually  warmer  when  the  melting  rate  is  higher. 

Figure  17  shows  the  distribution  of  wall  temperature  along  the 
lower  ice  surface  for  the  four  cases.  At  0°  tilt  the  wall  temperature 
has  the  expected  distribution,  coldest  at  the  bottom  and  then  warmest 
and  uniform  in  the  turbulent  region.  The  pattern  disappeared  when  the 
ice  was  tilted  and  the  most  striking  feature  is  the  overall  decrease  in 
interface  temperature  as  the  tilt  angle  increased.  At  the  mid-point  of 
the  ice  the  interface  temperature  decreased  from  -.41°C  at  0°  tilt  to 
-.82°C  at  60°  tilt.  The  spatial  variation  at  a given  angle  changed  to 
give  the  coldest  temperatures  at  the  ends  of  the  ice  and  the  warmest 
temperatures  in  the  middle  of  the  ice.  This  trend  was  evident  at 
tilts  of  15°  and  60°  but  not  at  32°  tilt. 

Plate  Xll  also  shows  the  melt  rate  is  higher  on  the  upper  side  of 
the  ice  than  on  the  lower  side  of  the  ice.  At  the  start  of  the  experi- 
ment the  ice  thickness  was  symmetrical  about  the  thermistor  array. 
Detailed  measurements  of  the  melt  rate  were  not  made,  as  this  experiment 
was  undertaken  only  to  determine  the  change  in  the  nature  of  the 
convection  as  the  angle  of  the  ice  block  changed. 
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The  tick  marks  denote  the  location  of  the  thermistors. 
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Concluding  Remarks 

These  laboratory  experiments  show  that  the  flow  generated  by  a 
melting  ice  wall  consists  of  both  laminar  and  turbulent  flows.  For 
oceanic  conditions,  the  vertical  extent  of  the  laminar  flow  next  to 
an  ice  wall  is  less  than  1.0  m for  Td  > 1°C.  Hence,  the  most  important 
flow  in  the  ocean  is  the  upward  flowing  turbulent  boundary  layer.  In 
the  following  chapters  I develop  the  fundamental  equations  and  boundary 
conditions  that  govern  the  turbulent  flow.  I then  simplify  the  equations 
with  a similarity  transformation  and  numerically  solve  the  resulting 
ordinary  differential  equations. 
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Fundamental  Equations 


The  governing  equations  of  turbulent  vertical  natural  convection 
are  those  of  a Boussinesq  fluid  subject  to  buoyancy  forces  that  arise 
from  the  cooling  and  dilution  at  the  ice-water  interface.  To  simplify 
the  governing  equations,  I parameterized  the  turbulent  transfer  of  salt, 
heat,  and  momentum  with  a spatially  varying  eddy  diffusivity  and  made 
the  boundary  layer  approximation  to  the  equations. 

Batchelor  (1967)  gives  the  mass  conservation  equation  for  a fluid 


as 


1 Dp. 
p Dt 


v 


u = 0 


(29a) 


Tennekes  and  Lumley  (1972)  give  the  mean  momentum  equations  of  a 
Boussinesq  fluid  in  two  dimensions  for  the  coordinate  system  illustrated 
in  Figure  18  as 

+ 3U  = _ 1 3£  + M)  + v/£u  + 3iu  ) 

3t  3x  3y  p 3x  prt  a (sx^  3y2  / 
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(29  b) 


IV  + ^ + v 3V 
3t  “sx  v3y 


where  u,  v are  mean  velocity  components;  u1,  v'  are  turbulent  velocity 
fluctuations;  v is  the  kinematic  viscosity;  p is  pressure;  g is  the 
gravitational  acceleration,  pQ  is  the  undisturbed  density;  and  p is  the 
density  perturbation  due  to  changes  in  the  salinity  or  temperature  of 
the  fluid. 
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Tennekes  and  Lumley  (1972)  show  that  the  terms  u 'u 1 and  v 'v ' 
represent  normal  stresses  or  pressures  and  have  little  effect  on  the 
transport  of  mean  momentum.  Hence,  I will  neglect  these  terms  in  Eqs. 
(29b,  c).  The  terms  u V and  v'uT  are  equal  and  represent  shear 
stresses.  I use  the  simplest  model  for  these  terms:  The  flux  of  mean 
momentum  due  to  the  turbulent  motions  equals  a spatially  varying  eddy 
diffusivity  times  the  mean  velocity  gradient  across  the  boundary  layer, 

uV  = TV  = A (30) 

where  A is  the  eddy  diffusivity  of  momentum,  a function  of  x and  y. 

The  turbulent  model  section  gives  a complete  discussion  of  the  spatial 
variation  of  A.  Substitution  of  Eq.  (30)  into  (29b,  c)  gives 
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The  equations  for  the  diffusion  of  heat  and  salt  complete  the  set 
of  fundamental  equations.  These  equations  also  contain  flux  terms  that 
are  modeled  with  an  eddy  diffusivity.  Following  Tennekes  and  Lumley 
(1972),  I take  the  eddy  diffusivity  for  momentum  and  a diffusing  element 
(heat  or  salt)  as  equal  because  turbulence  is  a flow  property  and  as 
such  transports  heat,  salt,  and  momentum  equally  and  the  transport  is 
independent  of  the  respective  molecular  diffusion  coefficients.  The 
mean  equations  for  the  diffusion  of  heat  and  salt  with  an  eddy 


diffusivity  are 
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T'.v  Souru.Lu'y  Layer  Apyroximaticn 

The  observed  laboratory  flows  fall  into  the  category  of  boundary 
layer  flow  and  the  boundary  layer  approximation  further  simplifies  the 
governing  equations.  In  the  boundary  layer  approximation,  scaling  the 
governing  equations  shows  that  downstream  gradient  terms  are  always 
smaller  than  cross-stream  gradient  terms  and  are  neglected  unless  they 
are  multiplied  by  the  downstream  velocity.  In  addition,  the  u-momentum 
equation  dominates  the  diffusion  of  momentum  and  the  v-momentum  equation 
is  neglected.  Also,  the  pressure  within  the  boundary  layer  equals  the 
pressure  outside  of  the  boundary  layer;  in  this  case  the  pressure  is 
hydrostatic.  To  arrive  at  the  boundary  layer  approximation,  scale  the 
continuity  equation  with  the  following  relationships 
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(33a) 
(33b) 
(33c) 
(33  d) 


where  - are  quantities  of  order  one  and  u*.  v*,  , and  ' are  the 
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velocity  and  distance  scales.  Substitution  of  (33)  into 

(29a),  eliminating  the  order  1 quantities,  and  omitting  the  partial 

time  derivative  gives 


1 

P 

A coordinate  transformation  to  a reference  frame  moving  with  the  ice- 
water  interface  eliminates  the  time  dependence  from  the  governing 
'•equations  and  a later  section  gives  a complete  discussion  of  the  trans- 
formation and  its  ramifications.  The  maximum  change  in  density  affecting 
the  left-hand  side  of  Eq.  (34)  occurs  when  the  interface  salinity  is 
zero  which  gives 

^ = 2 x 10-2  . (35) 
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Substitution  of  (35)  into  (34)  scales  the  density  variations 
to  give 


P 


(36) 


Because  of  the  small  value  of  , the  change  in  density  has  little 
influence  on  the  flow  and  the  left-hand  side  of  (29a)  is  neglected. 
The  mass  conservation  equation  then  becomes 
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and  the  flow  is  incompressible.  Substitution  of  Eqs.  (33)  into  (37) 


. | . _ ■ H I | 


95 


j 


gives 


where  o means  of  same  order.  Rearranging  (37)  shows  that 


(39) 


which  means  that  the  ratio  of  the  downstream  to  cross-stream  velocity 
must  be  of  the  same  order  as  the  ratio  of  the  downstream  to  cross- 
stream length  scales.  The  experiments  show  that  l = 1000  mm  and  ; 

x y 

= 20  mm;  therefore  the  downstream  velocity  is  50  to  100  times  the  cross- 
stream velocity. 

Substitution  of  the  scaled  variables,  (33),  into  the 
x-momentum  equation,  (31a),  results  in 


A term-by-term  comparison  shows  that  the  frictional  portion  of  the 
equation  can  be  simplified  because 


(41) 


but  neither  the  molecular  nor  turbulent  diffusive  terms  can  be  eliminated 
since  the  magnitude  of  the  eddy  diffusivity  is  unknown.  Both  advective 
terms  must  be  retained  as  they  are  the  same  order,  u;/^,  and  keeping 
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all  terms  of  this  order,  which  is  the  same  order  as  the  buoyant  term, 
the  x-momentum  equation  becomes 


li±+l3Uv^  * - 11£  + 

3t  3x  3y  p 3X 


(42) 


Schlichting  (1960)  carries  this  procedure  out  in  detail  and  shows  that 
the  scale  for  the  v-momentum  equation  is  v*/Z£  which  is  much  less  than 
u*/‘y  allowing  for  the  neglect  of  the  v-momentum  equation.  Schlichting 
also  shows  that  the  pressure  within  the  boundary  layer  equals  the 
pressure  outside  the  boundary  layer;  for  this  case  the  pressure  is 
hydrostatic,  and  given  by 

If  *'09  ' <«> 


Substitution  of  (43)  into  (42)  and  using  the  scaling  arguments 
to  simplify  the  heat  and  salt  diffusion  equations  gives  the 
simplest  form  of  the  governing  equations 


and 


3u 

3X 


+ 


(44a) 

(44b) 

(44c) 

(44d) 
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In  this  study,  a coordinate  transformation  to  a reference  frame 
moving  with  the  ice-water  interface  eliminates  the  time  dependence  and 
then  a similarity  transformation  converts  the  governing  partial 
differential  equations  to  ordinary  differential  equations.  But  first, 
to  complete  the  governing  equations,  the  boundary  conditions  must  be 
derived. 


Boundary  Conditions 


The  boundary  conditions  for  a melting  or  freezing  liquid-solid 
interface  follow  from  the  application  of  the  no-slip  condition  and 
conservation  of  mass,  heat,  and  salt  at  the  interface.  Figure  18  shows 
the  two-dimensional  geometry  used  in  the  derivation  of  the  boundary 
conditions  for  a melting  two-dimensional  curved  boundary,  o(x,  t). 

At  the  interface  the  no-slip  condition  requires  that  the  tangential 
velocity  be  zero.  The  velocity  components,  u and  v,  written  in  terms  of 
the  velocity  components  normal  and  parallel  to  the  interface,  up  and  vr> 
are 

u = ur  cos  Y - vr  sin  y , (45a) 

and 


v = ur  sin  y + vr  cos  y » 


(45b) 


where  y equals  the  angle  between  the  x-axis  and  the  tangent  to  the 
interface.  At  the  interface,  the  component  parallel  to  the  interface, 
ur,  is  zero  which  gives 


and 


u = -v  sin  y . 
v = vp  cos  y at  y = o(x,  t) 


Elimination  of  vr  from  (46)  gives 


(46a) 

(46b) 


u + v tan  y = 0 at  y s a(x,  t)  , (47) 
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and  recognizing  that  tan  > = |2-,  the  no-slip  boundary  condition  becomes 

o X 


u + v|^=0aty  = o(x.  t)  . 


(48) 


An  alternative  derivation  states  that  the  sum  of  the  projections  of  u 
and  v on  to  the  tangent  of  the  curve  y = o(x,t)  must  be  zero. 

I next  derive  the  boundary  conditions  which  follow  respectively 
from  conserving  mass,  heat,  and  salt  for  a control  volume  that  contains 
the  solid  liquid  interface.  The  cross-hatched  region  in  Figure  18  is  the 
control  volume  used  in  the  derivations  and  letting  the  horizontal  bounds 
of  this  control  volume  approach  the  interface  from  each  side  gives  the 
appropriate  boundary  condition. 

First,  mass  conservation  provides  an  additional  condition  on  the 
velocity  at  the  interface;  mass  conservation  is  written  as 


a_ 

at 


j jpdV  + 
V 


ndA  = 0 


(49) 


where  the  volume  and  surface  integrals  are  over  the  control  volume  shown 
in  Figure  A. 

Application  of  Gauss's  law  changes  the  surface  integral  to  a volume 
integral  and  noting  that  the  choice  of  x and  x + dx  is  arbitrary  and  not 
a function  of  time  gives 


x+dx 

dx' 

x 


3_ 

at 


b 

Pdy  + 


x+dy 
f dx' 


a x 


b 

’ 

v -(pu)dy  = 0 . 
a 


(50) 


For  brevity,  the  first  and  second  multiple  Integrals  in  (50)  a**e 


(00 


denoted  as  I and  II. 

Integration  of  I from  a to  0 and  then  from  0 to  b gives 

0 b 

1 = ft  {j  pdy  + f pdy)  ‘ (51) 

a 0 


Performance  of  the  integration  and  noting  that  a and  b are  not  functions 


of  time  yields 


x+dx 


I = 


/ 3o  3o\  . 

(pi  at  ' pw  3t)dx  • 


(52) 


Integral  II,  after  the  application  of  Leibnitz's  rule  with  the 
condition  that  the  velocities  in  the  ice  are  zero  yields 


x+dx  b 

\h\ pud*  -(pu>bk+  (pv)Udx 

x a 


(53) 


Takinq  the  limit  of  a and  b approaching  0 gives 

x+dx 

" ■ j !x)d*  • <M> 

x 

b 

The  integral  term,  ^ pudy,  vanishes  in  the  limit  because  a decomposition 
of  the  integral  a into  two  parts,  one  on  each  side  of  the  inter- 
face, shows  that  the  integral  within  the  ice  is  identically  zero  and  the 
integral  in  the  water  equals  zero  in  the  limit. 

The  sum  of  I and  II  gives  the  complete  boundary  condition.  The 
integrand  of  this  integral  must  equal  zero  because  the  bounds  of  the 
integral  are  arbitrary  and  after  some  rearrangement  the  boundary 
condition  becomes 
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(55) 


at  y = o(x,  t). 

The  physical  meaning  of  this  boundary  condition  is  most  easily 
understood  for  a flat  plate,  = 0.  For  this  case  the  velocity  normal 
to  the  wall  is  equal  to  the  normalized  density  difference  between  the 
liquid  and  the  solid  times  the  melting  or  freezing  rate.  For  an  ablat- 
ing ice  wall,  the  melting  of  a unit  mass  of  ice  yields  a unit  mass  of 
water  that  occupies  a 10T.  smaller  volume  than  the  ice.  This  volume 
decrease  is  exactly  balanced  by  a suction  velocity  toward  the  ice.  For 
freezing,  the  volume  expansion  caused  by  liquid  water  changing  to  ice 
expels  water  from  the  interface. 

Second,  the  derivation  of  the  thermal  boundary  condition  proceeds 
in  a similar  manner  using  heat  conservation  instead  of  mass  conservation 
at  the  interface.  I write  heat  conservation  as 


; 

i 

\ 


f 


a 

at  j 

’ 

pedV  + 

peu  • rid  A + 

i t 

l A 

F.  • ndA  = 0 
h 


(56) 


where  e is  the  internal  energy  of  the  liquid  or  solid,  F^  is  the  flux  of 
heat  due  to  diffusion,  u is  the  vector  velocity,  p is  the  density,  and  n 
is  the  unit  outward  normal.  I have  omitted  the  mechanical  production  of 
heat  because  Schlichting  (1960)  shows  that  this  effect  is  small  for  the 
slow  velocities  observed  in  these  flows.  Application  of  Gauss's  law  to 
the  control  volume  in  Figure  A gives 


X*dX  V JF  , 

dxl(-ri+^-dy*0  . 
V 3X  9y  ' 


( 57) 


For  simplicity,  the  three  integrals  are  called  in  order: 
1 1 1 ^ • Expansion  of  the  first  integral,  1^,  gives 


Ih,  1 1 h , and 


x+dx 

'h  ■ | (°1ei  * pwew)  If  dx  • 


(58) 


where  the  subscripts  i and  w refer  to  ice  and  water,  respectively.  This 
integral  represents  the  net  change  of  internal  energy  as  the  ice  melts. 

Next,  expansion  of  I Ih , which  represents  the  advection  of  heat  by 
the  velocity  field,  in  a manner  similar  to  the  derivation  of  (53) 

9,ves  x*dx  b 

nh  - \ (h  I ‘'budy  + peu|a  U - reu|b  §£  * pe»lb-  ,-evla)  dx  . 
x a 

(59) 

In  the  limit  of  a and  b approaching  o with  zero  velocity  in  the  ice, 

x+dx 

(60) 


11^  becomes 


Hh  = | pwew(v‘u  lf)y  = odx  • 


Finally,  the  expansion  of  1 1 1^ , which  represents  the  diffusive  heat 
flux,  gives 


1 
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x+dx 


+ (F> 


dx' 


(61) 


The  sum  of  1^,  1 1 ^ , and  1 1 1 ^ gives  the  thermal  boundary  condition; 
with  arbitrary  integration  limits  the  sum  of  the  integrands  must  equal 
zero;  hence,  the  thermal  boundary  condition  with 


F = -KvT 


(62) 


becomes 


(ciei  • Vw)itVw(v-uH)t*i5r-|<w5r 


3T 


3T. 


K 1?  ' Ki  1j)!f  = 0 at  y = o . 


(63) 


At  the  interface,  the  internal  energy  of  the  water  equals  the 
internal  energy  of  the  ice  plus  the  latent  heat,«£, 

ew  = ei  + . (64) 


Substitution  of  (55)  and  (64)  into  (63)  gives  the  complete 
thermal  boundary  condition 


‘ p1*!t  + Ki  3y 


3T4  9T 

- ^ — + 


3y 


/ 3T 

K 3T  ‘ K 


3T.  \ , 

Iff  = 0 at  y = O .(65) 


i 3x 


( 


This  condition  states  that  the  difference  in  the  amount  of  heat  entering 
one  side  of  the  interface  and  leaving  the  other  side  equals  the  product 
of  the  latent  heat  and  the  freezing  or  melting  rate.  Carslaw  and  Jaeger 
(1959)  call  this  condition  the  Stefan  boundary  condition. 
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The  third  boundary  condition  follows  from  salt  conservation  at  the 
interface  and  the  derivation  proceeds  in  the  same  manner  as  the  deriva- 
tion of  the  thermal  boundary  condition.  The  salt  conservation  equation 
contains  three  terms:  the  time  rate-of-change  of  salt  within  the 
control  volume,  and  both  the  advection  and  diffusion  of  salt  in  or  out 
of  the  control  volume.  Formally,  this  becomes 


3_ 

3t 


’ * 

pSdV  + ( 

Spu  • ndA  + | 

V 

J 

t 

l A 

ndA  = 0 


(66) 


As  before,  Gauss's  law  converts  the  surface  integrals  to  volume 
integrals,  then  the  integrals  are  expanded  using  the  identical  control 
volume  giving 


x+dx 


b 

pSdy  + 


x+dx 

dx 


b 

v 


x a x a 


•( 


pSu  + F. 


;)dy  - 


(67) 


Equation  (67)  is  the  saline  analog  to  (57),  the  expansion  of 
the  three  integrals  proceeds  in  a similar  fashion  with  the  additional 
assumption  that  the  ice  contains  no  salt.  For  icebergs  this  is  an 
excellent  assumption  because  glacier  ice  is  metamorphosed  snow  or  pure 
water.  For  multiyear  sea  ice,  either  pressure  ridge  keels  or  large 
floes,  observations  by  Cox  (1974)  show  that  the  ice  salinity  is  less 
than  4 °/00.  In  the  limit,  the  boundary  condition  becomes 


-S 


Sv  - Su  - D 


aS  + D 3S  ao 
ay  u 3x  ax 


0 at  y 


o.  (68) 


The  flat  wall  case  facilitates  the  interpretation  of  this 
boundary  condition;  in  this  case  the  condition  is 

s(|f  - v)  + = 0 at  y = O . (69) 

Equation  (69)  states  that  the  salt  flux  to  the  wall  equals  the 
product  of  the  wall  salinity  and  the  melting  rate.  Comparison  of 
Eqs.  (69)  and  (65)  shows  that  the  wall  salinity  is  analogous  to 
the  latent  heat  and  can  be  thought  of  as  a latent  salinity.  Unfortu- 
nately, the  wall  salinity  varies  and  depends  both  on  the  salt  and  heat 
fluxes  to  the  ice,  unlike  the  latent  heat  which  is  constant. 

Because  Eq.  (55)  shows  that  the  suction  velocity  at  the  ice  wall 
is  one-tenth  the  melting  rate,  I approximate  the  saline  boundary  condi- 
tion as 

s H ■ - 0 § at  * ■ ° • <70> 

I used  the  boundary  conditions  for  a flat  wall  to  discuss  and  model 
the  convective  motions  observed  adjacent  to  the  melting  ice  wall.  The 
curvature  does  not  affect  the  dynamics  of  the  fluid  motions  because  the 
radius  of  curvature  of  the  ice  surface  is  much  larger  than  the  length 
scale  of  the  flow.  Inspection  of  Plate  VII  which  illustrates  the 
ice  surface  after  being  submerged  for  220  minutes,  verifies  the  assump- 
tion of  an  almost  flat  plate.  Two  regions  where  this  approximation 
breaks  down  are  at  the  ice  bottom  and  in  the  'notch'  region  where  the 
laminar  layer  becomes  turbulent.  Both  the  notch  and  the  ice  bottom 
represent  the  leading  edge  of  an  upward  flowing  boundary  layer;  in  this 
case  the  turbulent  flow  and  the  inner  laminar  flow.  Therefore,  the 
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following  analysis  is  invalid  in  these  regions  but  valid  everywhere  else 
in  the  turbulent  flow. 

In  summary,  the  boundary  conditions  for  a flat  wall  are  at  y = o 

( 1 


T = - mS  , 

w w 


u = 0 , 


v = 


pw~p i Jo 

P at  ’ 
w 


i jT  ' " \ 


JT 

sy  ’ 


$ l£  /i  _ J?! i \ . o — 

w Jt  V c I u ay  • 


(71a) 

(71b) 

(71c) 

( 71  d) 

( 7 1 e ) 


and  at  y = 


u * 0 , 

( 7 1 f ) 

s = s , 

iV 

( 7 1 g ) 

T = T . 

«» 

( 71  h) 

The  conditions  at  y 3 o add  to  the  complexity  of  the  problem.  Only  the 
condition  u = 0 is  known  explicitly,  the  remaining  conditions  are  rela- 
tionships between  the  independent  variables  of  the  problem.  With 
boundary  conditions  of  this  form  the  solution  of  the  governing  equations 
becomes  difficult,  either  analytically  or  numerically. 


1 


i 
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Transformation  to  a Moving  Coordinate  System 


The  experiments  show  that  at  any  height  the  melting  proceeds  at  a 
constant  rate  and  the  flow  next  to  the  ice  is  in  steady  state  so  that 
a transformation  to  a coordinate  system  fixed  to  the  ice-water  interface 
removes  the  time  dependence  from  Eq.  (44).  The  constant  melt  rate  then 
appears  as  a fluid  source  at  the  interface. 

Writing  the  v-component  of  the  velocity  as  v = vr  + vt>  where  vr 
is  the  velocity  in  the  moving  reference  frame  and  vt  is  the  transla- 
tional velocity,  gives  for  the  left-hand  side  of  (44b) 


_3 

3t 


+ u 


3X 


IV 


Vt^  3y 


(72) 


If  the  temporal  changes  are  solely  due  to  advection  by  the  translational 
velocity,  then 

S*«tfs0  • <”> 

and  (72)  reduces  to 


3 , 3 

(i  + v 

ax  r ay 


(74) 


where  the  coordinate  y is  the  distance  from  the  ice  wall. 

In  the  moving  frame  the  normal  velocity  boundary  conditions  becomes 


vr  + vt 


pw  ~ pi  3o 

pw  ’ 


(75) 


4 


>4 
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and  since 


v 


r 


fj_  3o 
»w  ^ 


at  y = 0. 


(76) 


In  the  geometry  shown  in  Figure  18  melting  corresponds  to  a negative 

which  means  is  positive,  so  that  there  is  blowing  at  the  wall.  For 

freezing,  ^ is  positive,  vr  is  negative,  and  there  is  suction  at  the 

wall.  Since  ice  is  less  dense  than  water,  the  blowing  or  suction 

pi 

velocity  is  reduced  by  a factor  of  — . 

pw 

Another  way  to  conceptualize  the  effect  of  the  moving  coordinate 

system  is  to  fix  the  boundary  layer  in  space  and  then  translate  the  ice 

toward  the  boundary  layer  as  fast  as  it  melts.  This  introduces  mass 

into  the  boundary  layer  at  a rate  equal  to  p..  . Since  mass  cannot 

accumulate,  there  must  be  a velocity  at  the  interface  carrying  mass  away 

from  the  wall  at  a rate  of  p v . The  sum  of  these  two  mass  fluxes  must 

w r 

be  zero  so  that 


or 


3c 

at 


p v 
*Vr 


0 , 


y 

r 


_1  i£. 
pw  ^ 


at  y = 0 


(77a) 

(77b) 


which  is  again  (76). 

The  saline  boundary  condition,  (71e), 


as 

ay  * 


(78) 


also  changes  in  the  moving  coordinate  system.  Substitution  of 


into  (78)  gives 


s»(v  * H - !?)  = D f • 


Since 


vr(0)  = - 


3a 

»w  3t 


(79)  becomes 


or 


3o  _ D 3S 
'w  Pw  3t  ' u 3y 


3o  _ _L  IS 
n ~ ’ pi  ^ 3y 


(79) 


(80) 


The  remaining  boundary  conditions,  both  at  the  interface  and  at  a large 
distance  away  from  the  wall,  transform  unchanged.  Further  discussion 
will  take  place  in  the  moving  frame  and  for  convenience  the  subscript  r 
will  be  dropped. 


The  Similarity  Transformation 


To  simplify  the  governing  equations,  incompressibility  and  the  two 
dimensional  nature  of  the  flow  allows  for  the  replacement  of  u and  v 
with  derivatives  of  the  stream  function.  Then,  nondimensional  salinity 
and  temperature  variables  replace  S and  T in  the  equations.  Finally,  a 
similarity  transformation  reduces  the  governing  partial  differential 
equations  to  a system  of  ordinary  differential  equations  and  transforma 
tion  of  the  boundary  conditions  completes  the  similarity  development. 

The  stream  function,  v,  is  defined  as 


it 

ay 


and  v 


= . ill 


3X 


(81) 


Normalization  of  T and  S on  the  salinity  and  temperature  differences 
across  the  boundary  layer  gives  the  following  convenient  nondimensional 
formulation: 


S-S 


* = 


s -s 

w » 


(82a) 


T-T 


e = 


VTw 


(82b) 


Upon  substitution  of  (81)  and  (82)  into  (44),  the  governing  equations 
become 

it  i£t_  . it  iii  = ..  (p'po)  + 1 r \ a2 


ay  axay  ax  ay* 


ay  • 


(83a) 


iiii.itii=i_r  / n+a \ ili 
ay  ax  ax  ay  ay  ayJ 


(83b) 


Ill 


To  complete  the  set  of  equations,  the  equation  of  state,  described  in 
Appendix  I,  provides  the  final  relationship  between  p and  4 and  e.  Thus, 
(83)  and  the  equation  of  state  are  the  four  equations  that  govern  the 
variables  41,  4,  0 and  p. 

The  analysis  of  boundary  layer  flows  comnonly  employs  a similarity 
transformation  that  assumes  the  downstream  variation  of  the  flow  depends 
only  on  a scale  factor  determined  by  the  downstream  coordinate.  A 
transformation  in  this  case  is 


4 * axpf(n) 

(84a) 

A = xrN(n) 

(84b) 

<p  = 4>(n) 

(84c) 

6 = 0(n) 

(84d) 

- 

(84e) 

x 


where  f,  N,  4 and  0 are  unknown  functions,  n is  the  similarity  variable, 
a,  b,  p,  q,  and  r are  unknown  constants  to  be  determined.  In  this 
formulation,  I assume  4 and  e are  only  functions  of  n for  the  following 
reasons:  First,  my  laboratory  observations  show  that  for  uniform  far- 
field  conditions  S and  T have  no  measurable  x-variation.  Second,  the 

W W 

experiments  of  Vliet  and  Liu  (1969)  show  that  the  temperature  distribu- 
tion across  the  flow  depends  only  on  the  similarity  variable.  Finally, 
other  investigators  make  this  assumption  [cf.  Gebhart  and  Pera  (1971) 
and  Turner  (1973)]  with  reasonable  success. 
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Substitution  of  (84)  into  (83a)  yields 


a2b2x  2p-2q-l £(p_q)-f,2  -pff  ] 


p(<j>,6)-p  _ 

= -g  + ab3x  p-,3q  [(v+xrN)f,M  +xrN  'f  "]  (85) 


where  primes  denote  differentiation  with  respect  to  n. 

To  transform  (85)  into  an  ordinary  differential  equation,  the 
.equation  must  be  independent  of  x.  Inspection  of  the  buoyancy  term  in 
(85)  shows  that  this  term  is  only  a function  of  n because  neither  <j>  nor 
e are  explicit  functions  of  x by  the  equation  of  state.  Thus  by 
inspection  three  conditions  must  be  met: 


and 


2p  - 2q  - 1 = 0 , 

p - 3q  = 0 , 

p - 3q  + r = 0 . 


The  solution  of  (86)  yields 


p = 3/4,  q = 1/4,  and  r = 0; 


(86a) 

(86b) 

(86c) 


(87) 


as  a result,  the  eddy  diffusivity  is  a function  of  n only.  Then,  (85) 
becomes 

^ - [2f ' 2 - 3ff "] 

p($.e)-p0 


+ ab3[f ' ' A ‘ + ( v+A ) f 1 1 1 ] . 


(88) 
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To  simplify  (88)  define  a and  b by  the  equations 


a2b2  (pw'pJ 
4 " "9 


and 


vb3a  = -g 


Poo 

(pw"pj 


(89a) 

(89b) 


which  gives 


a = 4vb  = 4v 


and 


l 4v2  . 


l/4 


b = 


'9(pw-pJ 


4v2 


1/4 


(90a) 

(90b) 


With  these  formulations  for  a and  b and  letting  pq  = p (88) 
becomes 


pU>0)-p„  «i  . 

2f ' 2 - 3ff ' ' = + f ■ •'L  + n + ^)f  ■ ■ < 

Pw-Po  V ^ 


Substitution  of  (84)  into  the  thermal  and  saline  diffusion 
equations,  (83b)  and  (83c)  yields 

- pab  xP’9'1^'  = b2x"2q[(D+A)*' ' + /A'] 

- pab  xp"q_1fe'  = b2  x"2q[(,c+A)e ' ' + e'A’]. 


(91) 


(92a) 

(92b) 


Using  the  definitions  of  a and  b from  (90)  and  taking  p and  q from 


the  diffusion  equations  become 

“3f  4> 1 = (Sc  + £)*"  + 

(93a) 

“3fQ*  = (Pr  + *)e”  + e'^L  . 

(93b) 

and 
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Equations  (91),  (93a)  and  (93b)  are  the  transformed  diffusion  equations 
of  momentum,  salt,  and  heat.  If  there  is  no  turbulence,  then  these 
equations  reduce  to  the  laminar  equations  given  by  Schlichting  (1960). 

Transformation  of  the  boundary  conditions  and  confirmation  that 
these  conditions  satisfy  the  similarity  functions  completes  the  equations 
governing  the  turbulent  flow.  At  a large  distance  from  the  ice,  the 
vertical  velocity  is  zero  and  the  temperature  and  salinity  equal  their 
far-field  values,  or 

u = 0 , (94a) 

S = . * (94b) 

and 

T = at  y = - . (94c) 

In  n coordinates,  because  x is  held  fixed  as  y approaches  infinity,  the 
boundary  conditions  evaluated  at  y = oo  also  hold  at  n = =».  Transforma- 
tion of  the  far-field  boundary  conditions  (94)  gives  the  following: 

First  since 

, 1/2  (95) 

u = 4vb2  x1'  f 1 (n) 

and  because  v,  b,  and  x are  non-zero,  the  boundary  condition 
(94a)  becomes 

f = 0 at  n = oo  . (96) 

The  boundary  conditions  (94b  and  c)  when  used  with  the  definitions  of 
4>  and  e , (82)  yield 


4>  = e = 0atn  = 


oo 


(97) 
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The  Ice-water  Interface,  y = 0,  transforms  by  (84e)  to  n = 0; 
so  that  the  boundary  conditions  at  y = 0 are  applied  at  n x 0.  The 
no-slip  condition,  u = 0 at  y = 0,  using  (95),  gives 

f ' (0)  = 0 at  n = 0 . (98) 

The  salinity  and  temperature  at  the  interface  is,  respecti vely , and 

T which  gives 
w 3 

♦ = 0=1  at  ii  = 0 (99) 

from  (82).  Equation  (99)  appears  to  have  removed  the 
indeterminant  nature  of  the  freezing  curve  boundary  condition  from  the 
problem;  however,  this  condition  reappears  in  the  heat  and  salt  flux 
boundary  condition  which  is  discussed  next. 

The  thermal  and  saline  flux  boundary  conditions  in  the  moving 
reference  frame,  (80)  and  ( 71 d ) , combine  to  give 

1,001 


Use  of  definitions  Eqs.  (82)  and  changing  the  y-differentiation 
to  an  n-differentiation  through  (84e),  (100)  becomes 


♦ ■ "w  WT-> 
r - DLpw  T5-TJ 


at  y 


0 , 


(101) 


where  primes  denote  differentiation  with  respect  to  n.  The  freezing 
condition  now  appears  as  (101). 
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The  transformation  of  the  blowing  boundary  condition,  (76), 
completes  the  set  of  boundary  conditions  at  the  wall.  The  similarity 
transformation  gives  the  v-component  of  velocity  as 

v = -3vbx-1/4  f + vb2x-1/2  yf'  ; (102) 

at  the  interface  y = 0,  which  along  with  (76)  gives 


v(0)  = 


3o  _ 

at 


3vbx-1/4  f (o)  , 


(103) 


which  is  a function  of  x.  Solution  of  (103)  for  f(0)  yields 


f = 


Do  1/4 

pi  nx 


3vbp 


at  n = 0 


w 


(104) 


For  melting  ^ < 0 and  f ( 0 ) is  < 0;  for  freezing  > 0 and  f ( 0 ) is  > 0. 
The  assumption  that  the  spatial  variation  of  the  melt  rat&  is  propor- 
tional to  x'1/4  is  required  to  remove  the  x-dependence  and  the  similarity 
transformation  is  complete.  Letting 


|| . mox->A 


(105) 


gives  for  (104) 


f(0)  = 


(106) 


where  Mq  is  a characteristic  melting  rate  for  the  entire  ice  wall.  This 
formulation  of  the  melt  rate  restricts  the  spatial  variation  of  the  eddy 
diffusivity,  which  is  discussed  in  the  following  chapter. 


r 
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The  Turbulent  Model 

In  this  section  I model  the  turbulent  transports  processes  with  an 
eddy  diffusivity  that  depends  both  on  the  similarity  variable  n and  the 
buoyant  forcing.  The  bulk  measurements  made  in  these  experiments  permit 
the  evaluation  of  this  eddy  diffusivity.  This  simple  model  permits  the 
investigation  of  the  complex  effects  of  blowing  and  opposed  buoyancy 
forces. 

Townsend  (1976,  p.  259)  shows  that  the  turbulent  structure  of  a 
wall  jet,  a boundary  layer  flow,  consists  of  an  inner  and  an  outer 
region,  separated  by  the  location  of  the  mean  naximum  upward  velocity. 
For  my  case,  Figure  19  shows  these  two  regions  and  the  variation  of 
the  eddy  diffusivity  across  the  flow.  By  analogy  with  Townsend,  in  the 
inner  region,  which  lies  between  the  wall  and  the  position  of  the 
maximum  velocity,  the  turbulence  has  the  characteristics  of  wall  turbu- 
lence; the  eddy  diffusivity  increases  linearly  with  distance  from  the 
wall.  In  the  outer  region,  which  lies  outside  the  position  of  maximum 
velocity,  the  turbulence  has  the  characteristics  of  free  turbulence  and 
the  eddy  diffusivity  remains  constant. 

Further,  Townsend  shows  that  for  a smooth  wall  the  inner  region 
divides  into  a viscous  sublayer  near  the  wall  and  a region  where  the 
eddy  diffusivity  increases  linearly  with  distance  from  the  wall.  The 
viscous  sublayer  exists  because  the  ice  melts  in  a smooth  fashion;  if 
the  ice  melted  to  produce  a rough  surface  then  the  sublayer  would  be 
destroyed.  In  support  of  the  existence  of  the  viscous  sublayer  the 
temperature  records  from  thermistors  located  in  the  turbulent  flow 
(Figure  13)  show  a smooth  monotonic  increase  in  the  temperature 
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immediately  after  the  ice  melts  past  a thermistor,  then  as  the  ice 
melts  further  past  the  thermistor,  the  record  shows  turbulent  tempera- 
ture fluctuations.  Also,  particle  trajectories  and  dye  injections  show 
very  low  turbulent  activity  within  2 ntn  of  the  ice. 

The  outer  edge  of  the  viscous  sublayer  represents  the  region  where 
the  turbulent  transfer  of  momentum  approximately  equals  the  molecular 
transfer  of  momentum;  however,  turbulent  velocity  fluctuations  still 
persist  in  the  sublayer.  Townsend  (1976,  p.  354)  shows  that  for  larqe 
Pr  or  Sc  the  heat  or  salt  transfer  in  the  outer  region  of  the  viscous 
layer  results  from  turbulent  convection,  while  close  to  the  wall  the 
heat  or  salt  transfer  results  from  conduction.  Townsend  calls  this 
inner  region  the  'visco-conductive1  layer  and  the  position  of  the  outer 
edge  of  this  layer  occurs  where  the  turbulent  transfer  approximately 
equals  the  molecular  transfer  of  heat  or  salt.  Because  the  molecular 
thermal  diffusivity,  k,  is  much  larger  than  molecular  saline  diffusivity, 
D,  the  thermal  visco-conductive  layer  will  be  larger  than  the  saline 
visco-conductive  layer. 

To  incorporate  all  of  these  effects,  1 modeled  the  turbulence  with 
an  eddy  diffusivity  that  depends  only  on  the  similarity  variable  n. 

Near  the  wall  the  eddy  diffusivity,  A,  equals  zero,  corresponding  to  the 
region  of  no  turbulence.  In  the  wall  layer  the  eddy  diffusivity 
increases  linearly  with  increasing  distance  from  the  wall.  Beyond  the 
point  of  the  mean  maximum  velocity,  an  eddy  diffusivity  is  constant. 

The  mathematical  formulation  of  this  model  is 

A ( n ) = 0 for  0 < n £ An  , (107a) 

A(n)  = Ao[1-sech(a(n-An))] 


(107b) 
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illustrated  by  Figure  19.  The  region  of  no  turbulence  is  An,  a is  a 
shape  factor  and  sech,  the  hyperbolic  secant,  is  defined  as 

sech(y)  = — - . (108) 

eY-e'Y 

In  this  model  the  eddy  diffusivity  remains  constant  as  n increases; 
ideally  A should  approach  zero  as  the  velocity  decreases  with  increasing 
n.  Two  factors  allow  this  ambiguity  to  exist:  First,  the  outer  edge  of 
the  turbulent  boundary  layer  is  not  clearly  defined;  visual  observations 
of  the  entrainment  of  quiescent  fluid  by  the  turbulent  flow  show  the 
position  of  the  outer  edge  to  fluctuate.  Keeping  A constant  over  the 
region  of  entrainment  accurately  models  the  turbulent  processes  in  this 
region.  Second,  outside  the  region  of  entrainment  the  upward  velocity 
and  the  gradients  of  momentum  salt  and  heat  are  zero.  Because  the 
diffusive  flux  of  a quantity  equals  the  product  of  its  gradient  and  its 
total  diffusion  coefficient,  the  solution  to  the  governing  equations 
remains  unaffected  by  the  value  of  the  diffusion  coefficient  as  long  as 
the  gradient  equals  zero.  The  numerical  scheme  takes  advantage  of  this 
result  by  using  a larger  integration  step  size  away  from  the  wall  than 
near  the  wall.  If  A approached  zero  for  large  n,  then  the  large  Sc  of 
the  flow  would  require  a small  step  size  to  prevent  numerical  errors 
from  accumulating.  Betchov  and  Criminale  (1964)  used  a similar  eddy 
diffusivity  model  to  study  the  oscillations  of  a turbulent  flow  which 
they  state  as  applicable  to  turbulent  jets;  in  their  model  they  reduced 
the  eddy  diffusivity  in  the  far  field  to  1/3  its  maximum  value  to  avoid 
numerical  difficulties. 
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To  estimate  the  order  of  magnitude  of  A,  I next  compare  the 
observed  turbulent  flow  characteristics  with  the  theoretical  solutions 
of  Ostrach  (1954)  for  the  laminar  flow  generated  by  the  diffusion  of  a 
single  species  at  Sc  = 1000  and  the  results  show  that  A is  of  the  same 
order  as  v.  For  the  laminar  case  taking  - Sw  equal  to  15  o/OQ  gives 
the  following  boundary  layer  characteristics  at  a point  0.3  m downstream 
from  the  origin  of  the  flow;  the  melting  rate  equals  1.03  x 10'3mm  s'1, 
the  maximum  upward  velocity  equals  6 mm  s'1  and  boundary  layer  thickness 

. 

equals  30  mm.  In  my  laboratory  turbulent  case,  a value  of  15  °/00  for 
Soo  ' sw  corresponds  to  a thermal  driving  of  4°C,  and  the  boundary  layer 
has  the  following  characteristics  at  0.3  m downstream  from  the  bifurca- 
tion point:  The  melting  rate  equals  1.74  x 10"3  mm  s*1,  the  maximum 
upward  velocity  is  of  order  10  mm  s'1  and  the  boundary  layer  thickness 
equals  approximately  20  to  30  mm.  The  similarity  between  these  two 
cases  shows  that  A is  the  same  as  v and  molecular  effects  must  be 
included  in  a realistic  flow  model.  These  low  turbulent  levels  result 
from  the  large  Sc.  Even  in  the  turbulent  case  the  low  molecular  diffu- 
sivity  of  salt  contains  the  majority  of  the  saline  anomaly  near  the  ice 
wall  well  within  the  viscous  sublayer  where  the  frictional  effect  of  the 
wall  can  counter  the  large  buoyant  forces.  Also,  the  small  value  of  the 
Reynolds  number  for  the  observed  flows  suggests  weak  turbulent  mixing; 
the  laboratory  Reynolds  number  was  of  order  102. 

Because  of  these  low  turbulence  levels,  in  my  model  I retained  the 
molecular  effects  and  used  an  eddy  diffusivity  which  is  only  a function 
of  the  similarity  variable  n.  This  formulation 
points  in  the  boundary  layer  the  turbulence  has 
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This  assumption  may  not  adequately  represent  the  turbulence  at  all 
points,  but  the  low  turbulent  levels  and  the  large  Influence  of 
molecular  effects  allows  the  model  to  simulate  the  laboratory  flow. 
Modeling  the  turbulence  In  this  fashion  admits  a similarity  solution 
in  which  the  effects  of  opposing  buoyancy  forces  and  blowing  can  be 
readily  Investigated.  The  results  of  the  model  are  then  valid  for 
large  Sc  and  Gr  up  to  1013;  however,  for  an  iceberg  Gr  equals  10lp 
and  direct  application  of  the  model  to  icebergs  may  not  be  appropriate. 


Numerical  Scheme 


Because  the  governing  equations  and  the  boundary  conditions  are 
nonlinear  with  variable  coefficients,  I solved  them  numerically.  To 
solve  the  equations  I used  the  quasi-linearization  scheme  of  Radbill 
and  McCue  (1970)  programmed  on  the  CDC  6400  computer  at  the  University 
of  Washington.  Quasi-linearization  allows  for  the  linear  superposition 
of  the  particular  and  homogeneous  solutions  to  satisfy  the  boundary 
conditions. 

In  this  case,  specification  of  the  boundary  conditions  at  n = 0 
and  n = 00  gives  a two  point  boundary  value  problem  and  I applied  the 
boundary  conditions  at  n = 00  at  n = 20  which  correspond  to  approximately 
three  boundary  layer  thicknesses.  A guess  of  the  unknown  boundary 
conditions  at  n = 0 allows  for  the  integration  of  the  equations  over  the 
entire  interval.  At  n = 20  the  difference  between  the  calculated  and 
given  boundary  conditions  gives  the  disparity  between  the  iterated  and 
exact  solution.  With  linear  equations  these  differences  would  determine 
a new  set  of  coefficients  for  the  next  iteration  through  a linear  combi- 
nation of  the  homogeneous  and  particular  solutions.  In  a nonlinear  case 
Radbill  and  McCue  use  a Newton-Raphson  algorithm  to  generate  the  next 
iteration  which  is  a linear  function  of  the  previous  iteration  and  then 
use  linear  methods  to  generate  a new  set  of  coefficients  for  the  next 
iteration.  The  iteration  process  continues  until  the  difference  between 
the  iterations  is  less  than  a specified  value,  in  my  case  5 x 10"3. 

To  check  the  sensitivity  of  the  solution  to  the  convergence  parameter, 

I ran  two  duplicate  cases  with  the  convergence  parameters  equal  to 
5 x 10" 3 and  10" 3 respectively  and  found  no  discernable  differences. 
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I next  describe  the  Newton-Raphson  method  used  in  the  numerical 
integration  for  the  case  of  a nonlinear  scalar  equation 

9(x)  = 0 . (109) 

If  the  function  g(x)  decreases  monotonically  and  is  strictly  convex  in 
a sufficiently  large  neighborhood  of  x,  the  exact  root  g(x)  is  linearly 
approximated  at  x^  by  g[x^]  and  its  tangent  g'[x^], 

g(x)  a g[x(o)]  + [x-x(o)]g'[x(o)]  , (110) 


where  the  superscript  in  parentheses  denotes  the  iteration  number. 

f ‘ f ' 

Figure  20  shows  the  successive  approximations  of  g(x)  and  the  rapid 
approach  of  the  iterated  root,  x^,  to  the  actual  root.  Equation  (110) 
gives  a further  approximation  to  r when  g(x)  equals  zero  which  gives  the 
general  recurrence  relation: 


(n+1)  „ (n)  . g[x(n)l 

g:ixfnT] 


(111) 


Radbill  and  McCue  show  that  this  method  converges  both  quadratically  and 
asymptotically;  each  additional  iteration  doubles  the  number  of  correct 
digits  in  the  approximation.  They  also  generalize  the  Newton-Raphson 
method  to  a system  of  n nonlinear  ordinary  differential  equations  and 
give  several  detailed  examples  of  the  application  of  the  quasi- 
linearization technique  to  problems  in  fluid  mechanics. 

A useful  feature  of  this  scheme  is  a variable  integration  step  size 
which  I adjusted  to  minimize  the  computing  time  required  to  obtain 
convergence.  Table  5 gives  the  division  of  the  integration  Interval 
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and  the  size  of  the  steps  in  each  division;  each  step  was  further 
divided  by  15  In  the  actual  Integration.  The  high  gradients  of  velocity, 
salinity  and  temperature  near  the  wall  which  result  from  the  lack  of 
turbulence  and  the  high  Schmidt  number  require  a small  step  size  near 
the  wall  to  prevent  error  accumulation.  Further  from  the  wall  the 
gradients  relax  and  a larger  step  size  is  used. 


Table  5 


Interval  Division  Step  Size 
O.0  to  0.5  0.05 

0.5  to  4.0  0.25 

4.0  to  7.0  0.50 

7.0  to  20.0  1.00 


Table  5.  Step  size  in  each  division  of  the  Integration 
interval . 


Numerical  Results 


This  section  shows  how  the  numerical  scheme  permits  me  to  evaluate 
the  turbulent  parameters  which  determine  the  boundary  layer  character- 
istics. The  results  also  show  that  both  blowing  at  the  interface  and 
opposing  buoyancy  forces  must  be  included  in  the  analysis  to  model 
adequately  the  flow. 

Appendix  IV  contains  the  computer  printouts  of  f11,  f1,  f,  i', 
o',  o and  the  buoyancy  for  0 s n i 20  and  2°C  i Tj  i 12"C.  This  appen- 
dix contains  three  cases;  the  first  case  includes  the  effect  of  blowing 
and  opposed  buoyancy  forces.  The  second  case  omits  blowing  at  the 
interface  and  the  third  case  omits  the  temperature  effects  on  density. 

I evaluated  the  turbulent  parameters  by  varying  these  parameters 
until  the  calculated  values  of  ^'(0)  and  o'(0)  equaled  the  experimental 
values  calculated  from  the  smoothed  laboratory  data.  Table  6 gives 

4>jab(0),  °iab^’  dnc)  fHb^  for  1 Td  1 ^°C,  whore  t,u'  subscript 
lab  refers  to  the  laboratory  results.  The  calculations  at  each  T^  used 
fj  b(0)  at  that  T^  for  the  blowing  boundary  condition.  Also,  the 
resulting  flow  had  to  satisfy  two  other  less  exact  but  equally  important 
conditions:  The  entire  flow  must  be  upward  and  the  boundary  layer  must 
be  approximately  20  mm  or  5 n units  thick. 

Once  this  method  succeeded  at  T^  = 'TC,  1 found  that  a parabolic 
dependence  of  the  magnitude  of  the  eddy  diffusivity  on  the  density 
difference  between  the  wall  and  the  far-field  accurately  modeled  the 
flow  from  2°C  > T^  i 12°C.  For  T^  3 14°C  the  solution  showed  a signi- 
ficant downward  velocity  and  the  calculated  values  of  ^‘(0)  and  o'(0) 
disagreed  with  the  laboratory  values.  For  T^  I 16°C  the  numerical 
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Table  6 


a b 


Tdro 

♦'Ub<°) 

f'lal)(0)xl°-' 

( mm 5 / ‘*  s_1) 

(mm-3/1 

2 

-5.02 

-0.556 

-0.26 

11.0 

1.53 

4 

-4.67 

-0.720 

-0.68 

12.8 

1.78 

6 

-3.83 

-0.837 

-1.24 

13.7 

1.90 

8 

-2.93 

-0.937 

-1.92 

14.2 

1.97 

10 

-2.04 

-0.982 

-2.59 

14.4 

2.00 

12 

-1.33 

-0.989 

-3.22 

14.5 

2.02 

Table  6. 

^’lab^'  0 

' lab(O)*  ^ lab^’  a anc*  ^ comPutecl  from  The 

smoothed  laboratory  data. 
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scheme  failed  to  converge,  which  shows  that  at  high  Td  the  flow 
characteristics  change  radically.  This  result  is  in  fair  agreement 
with  the  observed  flow  reversal  found  at  Td  > 20°C.  In  the  following 
section  I use  the  model  results  at  Td  = 12°C  to  form  zero  net  buoyancy 
condition  that  accurately  predicts  the  turbulent  flow  reversal. 

To  evaluate  the  turbulent  parameters,  which  are  A,  the  maximum 
value  of  the  eddy  diffusivity;  An,  the  thickness  of  the  laminar  region; 
and  a,  the  shape  factor,  I set  a equal  to  1.5  and  then  varied  A and  An 
in  a systematic  fashion  until  the  calculate  values  of  4>'(0)  and  o'(0) 
equaled  the  laboratory  values  shown  in  Table  6.  This  technique 
rapidly  yielded  values  for  A and  An  because  changes  in  each  parameter 
affected  the  salinity  and  temperature  distributions  differently. 

I chose  a = 1.5  to  allow  the  eddy  diffusivity  to  reach  its  maximum 
value  in  the  region  of  the  maximum  velocity,  as  required  by  Townsend 
(1976).  Figure  21  gives  the  f‘  profiles  from  the  final  turbulent 
model  with  « = 1.5  and  shows  that  the  maximum  velocity  occurs  for 
0.7  ini  2.0.  The  position  of  maximum  velocity  remained  virtually 
unaffected  when  I varied  a from  0.5  to  2.0,  however  1 could  not 
satisfy  the  thermal  and  saline  boundary  conditions.  The  eddy 
diffusivity  attains  90%  of  its  maximum  value  when  an  = 3.0;  therefore, 
this  value  of  a places  the  maximum  eddy  diffusivity  at  the  required 
location  for  Td  > 2°C.  For  Td  = 2°C,  the  eddy  diffusivity  in  the  region 
of  maximum  velocity  is  much  less  than  the  molecular  viscosity  and  hence, 
the  eddy  diffusivity  has  little  direct  effect  on  the  distribution  of 
momentum.  To  simplify  the  calculations  over  the  entire  temperature 
range,  I kept  a constant  for  all  Trf. 
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Figure  21 . 


131 


With  a fixed  at  1.5,  I then  varied  A and  An  until  <j>'(0)  = i)>.jab(0) 
and  0 ' (0)  = e-jg^O)  at  = 4°C.  Figure  22  shows  the  dependence  of 

4>'(0)  and  0 ' ( 0 ) on  A and  An  with  a = 1.5  at  = 4°C,  where  a more 
negative  value  of  either  <t> ' ( 0)  or  e 1 ( 0 ) corresponds  to  a higher  salinity 
or  temperature  gradient  at  the  interface.  Variations  of  A have  a larger 
effect  on  e'(0)  than  on  $'(0)  while  the  opposite  holds  true  for  varia- 
tions of  An;  this  is  a result  of  the  large  Lewis  number  that  character- 
izes these  flows.  The  thermal  anomaly  with  its  high  diffusivity, 
diffuses  through  the  laminar  region  into  the  turbulent  region  where  the 
eddy  diffusivity  affects  the  heat  transfer.  The  low  value  of  the  saline 
molecular  diffusivity  contains  the  saline  anomaly  close  to  the  ice; 
hence,  changes  in  the  thickness  of  the  laminar  region  greatly  affect  the 
salt  distribution  within  the  boundary  layer.  A change  in  An  also 
affects  the  temperature  distribution  for  two  reasons:  First,  the 
temperature  distribution  is  affected  in  the  same  manner  as  the  salt 
distribution  only  to  a lesser  extent  because  of  the  large  Lewis  number. 
Second,  the  density  dependence  on  salinity  produces  large  changes  in  the 
buoyant  forcing  which  affect  the  flow  and  its  transfer  characteristics. 

Figure  22  shows  that  A = 0.3  mm2  s"1  and  An  = 0.08  almost  give 
the  required  values  of  <j>'(0)  and  e ' (0)  for  Trf  = 4°C.  For  the  model  at 
Tj  = 4°C  I chose  A and  An  equal  to  0.28  mm2  s_1  and  0.078,  respectively, 
for  the  turbulent  parameters  because  these  values  gave  an  average  error 
in  the  salinity  and  temperature  gradients  of  2.5%,  and  an  average  error 
of  less  than  5%  for  an  extension  of  the  turbulent  model  over  the  range 
of  Td  from  2°C  to  12°C. 

To  extend  the  turbulent  model  over  the  temperature  range,  I took 
the  following  formulation  for  A: 
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/P  ~ p ,n 

A = An(— 

o\Pco  - Pf/ 


(112) 


where  pw  is  the  interfacial  water  density,  p^  is  the  fresh  water  density 
at  T = 0,  and  n is  an  undetermined  exponent.  I choose  (112)  because 


the  experimental  melt  rates  indicate  a turbulent  transfer 


mechanism  that  increases  with  increasing  density  difference  between  the 


far-field  and  the  interface.  The  depencence  of  pw  on  makes  A a 


function  of  Td  and  the  experimental  results  show  that 


p„  - p 


p - p 


- = tanh(0.15  Td) 


(113) 


because  density  depends  linearly  on  salinity.  This  formulation  gives  a 


rapidly  increasing  A for  Td  up  to  9°C  and  an  almost  constant  A for 
Td  > 9°C,  a variation  suggested  by  the  observed  melt  rates.  Figure  14 
gives  the  experimental  melt  rate  dependence  on  Td>  Below  9°C  the  melt 


rate  increases  as  T^*6  which  indicates  that  transfer  of  heat  and  salt 


not  only  increases  as  the  result  of  higher  temperature  and  saline 


differences  across  the  boundary  layer,  but  also  the  turbulent  transfer 


mechanism  becomes  stronger.  Above  9°C  the  melt  rate  increases  linearly 


which  indicates  that  the  turbulent  transfer  mechanism  remains  constant 


while  the  temperature  difference  across  the  boundary  increases  linearly. 


To  determine  the  coefficient  n,  I numerically  solved  the  governing 


equations  at  Td  equal  to  2°C,  4°C,  and  8°C  for  three  values  of  n:  1 


1.5,  and  2.  with  An  held  constant  at  0.09.  For  each  value  of  n,  I held 


An  constant  at  0.09  and  I adjusted  AQ  to  give  A at  4°C  equal  to 


0.28  mm2  s"1.  Figure  23  compares  the  smoothed  experimental  data  with 


1 


a‘  (o) 


I 

l 


points  for  ♦'(0)  and  e'(0)  from  the  experiments, 
and  from  four  models  of  the  eddy  diffusivity; 
for  discussion. 
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the  results  from  each  set  of  calculations.  The  heavy  solid  line  Is  the 
locus  of  experimental  points.  Each  tick  mark  on  this  line  represents  the 
boundary  conditions  for  consecutively  increasing  T^,  beginning  at  2°C 
at  the  lower  right  and  ending  at  12°C  at  the  upper  left.  The  arbitrary 

division  of  melt  rate  dependence  on  Td  produces  the  kink  at  Td  = 9°C. 

a 

The  solutions  for  n = 1,  denoted  by  x,  give  temperature  and  salinity 

gradients  that  are  too  low  at  T^  = 8°C  and  too  high  at  2°C.  The 

solutions  for  n = 3/2,  denoted  by  A,  are  in  better  agreement  with  the 

laboratory  results.  The  solutions  for  n = 2,  denoted  by  Q,  give 

results  that  all  lie  above  the  experimental  line  which  indicates  a 

change  in  An  would  give  better  overall  agreement.  Values  of  0.078  for 

An  and  0.922  mm-'  s"1  for  Aq  give  the  results  denoted  by  0 where  the 

numbers  inside  the  0’s  are  the  T.  for  each  calculation.  These  are  the 

d 

final  parameters  of  the  model  and  they  give  an  average  error  of  less 

than  5%  over  the  ranqe  of  T.. 

d 

The  final  form  of  the  turbulent  model  is 

a = 1.5  , 

An  = 0.078  , 

and  A = 0.922[tanh(0.15  T^)]2  mm;  s"1  (114) 

which  gives  an  average  error  of  less  than  5$  for  2°C  < T^  < ^"C. 

Using  these  values,  I next  discuss  the  profiles  of  salinity, 
temperature  and  velocity  across  the  boundary  layer  and  the  variations  in 
these  profiles  that  result  from  a change  in  the  thermal  driving. 

Figure  24  gives  the  profiles  of  * and  e for  a thermal  driving  of 
4°C,  8°C,  and  12°C.  The  two  most  important  features  are  the  difference  in 
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the  saline  and  thermal  boundary  layer  thicknesses  and  the  thickening  of 
each  boundary  layer  as  the  thermal  driving  increases.  The  difference  in 
the  boundary  layer  thicknesses  results  from  the  low  turbulence  levels 
and  the  large  Lewis  number.  The  laminar  region  adjacent  to  the  ice  and 
the  small  value  of  D act  to  contain  95%  of  the  saline  anomaly  to  within 
1 n unit  of  the  ice  wall.  Although  most  of  the  saline  anomaly  is 
contained  near  the  ice,  the  eddy  diffusivity, which  increases  with 
distance  from  the  wall,  allows  the  dilution  effect  to  diffuse  out  to  4 
or  5 n-  The  eddy  diffusivity  reaches  the  same  order  of  magnitude  as  D 
at  a distance  of  approximately  0.01  n from  the  ice  and  this  effect 
produces  the  sharp  bend  in  the  <t>  profiles  at  0.5  n.  The  thermal 
boundary  layer  thickness  is  4 to  5 n units  and  the  e profiles  do  not 
exhibit  the  same  sharp  bend  found  in  the  $ profiles.  This  behavior 
results  from  a maximum  eddy  diffusivity  that  is  less  than  or  of  the  same 
order  as  the  molecular  thermal  diffusivity. 

Next,  Figure  21  gives  the  nondimensional  velocity  or  f‘  profiles 

for  Td  = 2°C,  4°C,  and  8°C.  The  profile  for  Td  = 12°C,  which  is  almost 
identical  to  that  of  Td  = 8°C,  has  been  omitted  but  Appendix  IV 
contains  values  of  f'(n)  for  Td  = 12°C.  The  similarity  in  profiles 
results  from  the  almost  identical  buoyant  forcing  in  each  case.  The 
difference  between  the  far-field  and  interfacial  conditions  determines 
the  overall  buoyant  forcing  at  each  Td>  When  Td  increases  from  8°C  to 
12°C,  the  laboratory  results  show  that  the  saline  difference  across  the 
boundary  layer  increases  by  10%  while  the  thermal  difference  increases 
by  50%.  This  large  change  in  temperature  difference  along  with  the 
increasing  effect  of  temperature  on  density  creates  enough  downward 
buoyancy  to  counter  the  increased  upward  buoyancy  that  results  from 
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increased  salinity  difference.  Hence,  as  increases  from  8°C  to  12°C, 
the  velocity  profile  undergoes  little  change. 

Figure  21  shows  that  the  momentum  boundary  layer  thickness, 

6 , is  constant  for  all  T..  definition  of  $ as  the  distance  from 
m d m 

the  ice  where  f‘  drops  to  90%  of  its  maximum  value  gives  <5  approxi- 
mately 11  n.  This  value  for  6 of  11  n corresponds  to  approximately 
25  mm  which  agrees  with  the  laboratory  observations;  however,  the 
long  tails  of  the  f'-profiles  make  determination  of  the  boundary 
layer  thickness  arbitrary.  The  constant  6m  for  varying  Td  results 
from  the  low  values  of  A which  are  always  less  than  the  molecular 
diffusivity,  v.  At  all  Trf  v determines  the  boundary  layer  thickness 
and  the  turbulence  has  only  an  indirect  effect  on  the  distribution  of 
momentum.  The  turbulence  acts  to  transport  the  dilute  water  near  the 
ice  away  from  the  ice  where  this  buoyant  water  generates  upward  momentum 
away  from  the  strong  frictional  effect  of  the  wall.  This  dilute  water 
is  approximately  0.5  °/00  less  saline  than  the  ambient  water  but  the 
large  effect  of  salinity  on  density  produces  enough  buoyancy  to  over- 
whelm the  negative  thermally-produced  buoyancy.  The  numerical  solutions, 
given  in  Appendix  IV  show  that  cooling  produces  negative  buoyancy  of 
order  10“2  mm  s-2.  Because  the  velocity  or  f'  profiles  show  no  reversal, 
the  combined  molecular  and  turbulent  diffusion  of  positive  momentum  into 
the  region  of  negative  buoyancy  overwhelms  this  small  amount  of  negative 
buoyancy. 

A comparison  of  the  f*  profiles,  given  in  Figure  21,  to  either 
the  4>  profiles  or  9 profiles,  given  in  Figure  24  shows  that  the 
momentum  boundary  layer  is  always  thicker.  Unlike  <sm>  6^,  and  6^ 
increase  in  thickness  for  increasing  Trf  and  this  thickening  manifests 
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itself  as  longer  tails  in  each  profile. 

Figure  21  also  shows  that  the  maximum  value  of  f ’ , f'  , 

fnaX 

increases  with  increasing  T.  up  to  T.  = 8°C.  Above  T , = 8°C,  f' 

3 d d d max 

remains  almost  constant  at  2.5  x 10"2.  Table  7 gives  f'  for 

nJaX 

2°C  * < 12°C  for  three  cases;  namely,  density  as  a function  of  both 

S and  T;  density  is  only  a function  of  S,  and  density  is  a function  of 
S and  T with  no  blowing  at  the  interface.  In  the  case  of  no  blowing, 
for  Td  > 10°C,  the  numerical  scheme  failed  to  converge.  The  values  of 
fnax  for  p(S,  T)  increase  rapidly  with  Td  up  to  8°C  where  they  reach  an 
almost  constant  value  and  then  slightly  decrease  at  12°C.  The  increase 
results  from  greater  buoyancy  at  increasing  Td  but  the  available  buoy- 
ancy reaches  a maximum  as  S approaches  zero.  Figure  14  from  the 

W 

laboratory  results,  shows  that  S approaches  zero  in  the  neighborhood 

W 

of  Td  = 10°C;  therefore,  f^ax  should  approach  a constant  value  near  10°C. 

Above  Td  = 10°C,  the  increasing  thermal  effect  on  density  causes  the 

slight  decrease  of  f*  at  T.  = 12°C. 

3 max  d 

The  increase  of  f'  is  not  linearly  related  to  the  available  buoy- 
ancy  because  the  increased  mixing  at  higher  Td  also  increases  the  upward 
mass  transport.  A comparison  of  the  first  two  columns  of  Table  7 
shows  that  the  disproportionate  increase  in  f^  results  from  the 
increased  mixing  and  not  the  increasing  production  of  downward  thermal 
buoyancy.  The  second  column  shows  the  case  with  the  density  is  only  a 
function  of  salinity  and  the  same  disproportionate  increase  is  observed 
in  addition  to  higher  values  of  f'  . These  higher  values  result  from 
the  elimination  of  the  downward  thermally  produced  buoyancy  force  and  a 
later  section  of  this  chapter  discusses  the  effect  of  both  thermal 
buoyancy  and  blowing  at  the  interface  of  the  turbulent  boundary  layer. 


Table  7 
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Table  7 also  gives  the  nondimenslonal  upward  mass  transport, 
f(®),  for  the  same  three  situations  and  f(®)  exhibits  the  same  variation 
as 

max 

To  compute  the  upward  transport  per  unit  width,  (84a)  evaluated  at 
y = * gives 

~ ax1/4  f(°°) 

where  x is  the  vertical  length  of  the  turbulent  boundary  layer  and 
Table  6 gives  values  of  a for  2°C  - 5 12°C.  As  an  example,  in  the 

laboratory,  the  upward  transport  by  a 0.5  m long  boundary  layer  next  to  a 
block  of  ice  0.2  m wide  after  1 hour  is  0.13  m3.  During  an  experiment 
there  are  two  sides  of  the  ice  block  melting  at  the  same  rate;  therefore, 
the  transport  must  be  doubled  to  give  0.26  m3.  Because  the  experimental 
tank  has  a horizontal  cross  sectional  area  of  1.2  m x 0.45  m,  after  one 
hour  the  water  in  the  upper  0.48  m is  water  upwelled  by  the  boundary 
layer.  To  counter  this  vigorous  upflow,  I used  the  pumping  and  heating 
systems  previously  described  to  maintain  uniform  far-field  conditions. 

A final  example  applies  these  results  to  an  iceberg  with  a 100  meter 
vertical  length  scale  floating  in  seawater  at  a thermal  driving  of  2°C. 
The  resulting  upward  volume  transport  equals  294  m2  day-1  per  waterline 
meter.  If  the  iceberg  has  a circumference  of  1 km,  the  upwelling  would 
equal  2.9  x 105  m3  day. 

The  Effect  of  Blowing  at  the  Interface 

The  following  section  discusses  the  effects  of  blowing  at  the  ice- 
water  interface  which  are  greatest  on  the  salinity  distribution.  To 
determine  the  blowing  or  melt  water  injection  effects,  I set  the  blowing 
or  f ( 0)  equal  to  zero,  then  numerically  solved  the  differential 
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equations  (91,  92a, b)  at  = 2,  4,  and  8°C  with  the  turbulent 
parameters  previously  used  at  the  respective  T^'s. 

Figure  25  displays  the  effect  of  no  blowing  on  f'(0)  and  e'(0); 
the  circled  numbers  are  the  results  from  each  calculation  and  the  broken 
line  connects  the  points.  For  comparison,  the  solid  line  connects  the 
points  from  the  calculations  with  blowing  at  the  interface.  The  results 
show  that  as  increases,  4>‘(0)  and  o'(0)  become  more  negative  which 
corresponds  to  both  higher  salinity  and  heat  fluxes  to  the  ice;  however, 
the  experiments  show  that  with  increasing  the  heat  flux  increases 
while  the  salt  flux  decreases. 

In  summary,  the  reduction  in  both  S and  salt  flux  to  the  ice 

w 

results  from  blowing  at  the  interface  and  the  low  saline  molecular 
diffusivity.  The  blowing  indirectly  effects  the  momentum  distribution 
in  the  boundary  layer  and  the  heat  transfer  through  the  large  effect  of 
salinity  on  density.  As  an  interesting  note,  Radbill  and  McCue  (1970) 
show  that  blowing  is  used  to  protect  space  vehicles  from  thermal 
extremes  during  re-entry  into  the  atmosphere.  These  vehicles  are  coated 
with  a substance  that  vaporizes  upon  heating  which  creates  a blowing 
velocity  that  is  large  enough  to  inhibit  the  flow  of  heat  toward  the 
vehicle. 

The  Effect  of  Opposing  Buoyancy  Forces 

To  investigate  the  effect  of  opposing  buoyancy  forces  on  the  model, 
I solved  the  governing  equations  with  the  buoyancy  determined  by  a 
linear  dependence  of  density  on  salinity  only 
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Figure  25.  The  effects  of  blowing  and  opposing  buoyancy  forces  on 
<j>'(0)  and  e ' (0)  at  various  Tj;  results  including  both 

blowing  and  opposing  buoyancy  forces,  ; no 

blowing, ; and  density  only  a function  of 

salinity,  . 
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I retained  blowing  at  the  interface  and  used  the  same  turbulent 
parameters  that  gave  agreement  between  the  laboratory  and  theoretical 
work.  The  only  difference  between  this  model  and  the  first  model  is 
that  in  this  case  temperature  diffuses  as  a passive  contaminant. 

Figure  25  shows  the  effect  of  opposing  buoyancy  forces  on  ^'(0) 
and  o ' (0)  for  Td  from  2°C  to  10°C  which  are  the  points  in  squares 
connected  by  the  dashed  line.  At  all  T^,  the  values  of  ^'(O)  remain 
almost  unchanged  while  the  values  of  o'(0)  become  more  negative  which 
indicates  a melt  rate  higher  than  the  observed  melt  rates.  With 
buoyancy  independent  of  temperature,  the  upward  flow  becomes  more 
energetic  which  increases  the  heat  transfer  to  the  ice.  The  salt 
transfer  remains  unchanged  because  the  flow  near  the  ice  ' saline 
anomaly  is  high  remains  virtually  unchanged. 

Table  7 gives  the  effect  of  opposing  buoyancy  fora.  n and 

f(“>).  The  retention  of  the  temperature  effect  on  density  reduces  both 

quantities.  At  T,  = 2°C  the  omission  of  thermal  effects  increases  f' 

^ d max 

by  15%  and  f(>")  by  46%;  at  T.  = 10°C,  f‘  increases  by  44%  and  f (■>■ ) 
increases  by  32%.  These  large  changes  result  directly  from  the  thin 
salinity  boundary  layer  which  allows  the  thicker  temperature  boundary 
layer  to  reduce  the  net  buoyancy  driving  the  flow. 

In  conclusion  the  numerical  results  show  that  the  turbulent  model 
consisting  of  the  three  regions  described  in  the  previous  chapter 
accurately  models  the  observed  flow.  In  addition,  the  effect  of  blowing 
at  the  interface  and  opposing  buoyancy  forces  must  be  included. 

Previous  turbulent  studies,  either  experimental  or  theoretical,  have  not 
included  either  of  these  effects. 


The  High  Temperature  Reversal 


The  high  temperature  reversal  of  the  turbulent  boundary  layer 
results  from  three  factors.  First,  in  the  turbulent  boundary  layer 
the  model  shows  that  the  difference  in  molecular  diffusi vities  and  the 
lack  of  turbulent  mixing  near  the  wall  produces  a saline  boundary  layer 
that  is  thinner  than  the  thermal  boundary  layer.  Second,  the  large 
temperature  difference  across  the  boundary  layer  produces  more  dense 
water  in  the  thermal  layer.  Finally,  the  effect  of  temperature  on 
density  increases  as  the  temperature  increases.  These  factors  combine 
to  produce  a small  amount  of  much  loss  dense  fluid  inside  of  a large 
amount  of  slightly  more  dense  fluid.  For  Tx  less  than  the  reversal 
temperature,  T , enough  low  salinity  and  positive  momentum  water 
diffuses  away  from  the  ice  to  overwhelm  the  thermal  effect  on  density; 
as  Tv  increases  past  T , the  cooling  produces  a large  amount  of  slightly 
more  dense  water  that  overwhelms  the  small  amount  of  very  less  dense 
water. 

Integration  of  the  density  distribution  across  the  boundary  layer 
gives  the  net  density  anomaly,  B,  within  the  boundary  layer. 


B = 


CO 

0 


[p„  - p(x)] 

— dy 


(116) 


The  use  of  a linear  equation  of  state. 


p = pjl  + B(S-SJ  + «(T  - TJ]  , (117) 


with  a evaluated  at  T and  S simplifies  the  calculations  and  gives  the 
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maximum  effect  of  temperature  on  density  and  results  in 

oo 

B = | [B(S-SJ  + a(T-Tj]  dy  . (118) 

> o 

Elimination  of  T and  S from  (118)  through  substitution  of  (82a, b),  the 
definition  of  ♦ and  6,  into  (138)  and  setting  B equal  to  zero  yields 
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(119) 


At  high  temperatures  the  experiments  show  that  S and  T approach  zero, 

w w 

which  gives 
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In  terms  of  the  similarity  variable,  n,  (120)  becomes 
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(121) 
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after  eliminating  the  common  similarity  scales.  The  ratio  of  two 
integrals  taken  from  the  model  at  Td  = 10°C  and  12°C  is  0.145  and  0.157. 
Because  the  buoyant  driving  is  constant  above  10°C  T^,  the  profiles 
remain  unchanged  except  for  the  effect  of  the  increased  blowing  which 
is  greatest  near  the  wall;  so  I took  0.157  as  a typical  value  of  the 
ratio  to  use  in  the  computations.  I also  used  the  equation  of  state 
of  Gebhart  and  Mollendorf  (1977)  to  calculate  a^,  Ta,)  which  resulted 


in  a transcendental  equation  in  for  a given  Sot.  An  iterative 
technique  solved  (121)  for  Sai  up  to  40  °/00  and  Figure  26  shows  the 
position  of  the  zero  net  buoyancy  curve,  labeled  Z.  For  far-field 
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Figure  26.  A T-S  diagram  showing  the  far-field  conditions  of  the 
turbulent  experiments;  the  arrows  indicate  direction 
of  flow.  Flow  reversal  predicted  by  a simple  mixing 

model,  , and  by  the  buoyancy  integrated  across 

the  boundary  layer  which  is  the  curve  label  z. 
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conditions  below  Z,  the  net  buoyancy  is  upward  and  for  far-field  condi- 
tions above  Z,  the  net  buoyancy  is  downward.  The  accuracy  of  equation 
(121)  decreases  with  decreasing  S because  the  linear  equation  of  state 
becomes  increasingly  erroneous  at  low  salinities. 

The  points  on  Figure  26  locate  the  various  experiments  and  the 
arrow  at  each  point  indicates  the  direction  of  the  observed  turbulent 
flow.  In  all  cases,  the  direction  of  flow  agrees  with  the  direction 
of  the  net  buoyancy;  for  the  two  experiments  nearest  Z the  mean  flow 
took  much  longer  to  establish  itself  than  in  any  of  the  other  experi- 
ments. For  all  the  upward  flowing  cases  the  ice  melted  smoothly  and 
there  was  a laminar  bidirectional  at  the  bottom  of  the  ice.  For  all  of 
the  downward  flowing  cases  cusps  formed  in  the  ice  and  there  was  a 
laminar  bidirectional  flow  at  the  top  of  the  ice.  The  open  circle  points 
are  experiments  performed  with  Dr.  Herbert  Huppert  while  he  was  visiting 
the  laboratory  during  September  1978. 

The  dashed  line  in  Figure  26  locates  the  far-field  conditions  for 
zero  net  buoyancy  assuming  1 gm  of  ice  melts  in  a volume  of  salt  water 
that  contains  just  enough  heat  to  melt  the  gram  of  ice.  Conservation 
of  heat,  mass,  and  salt  along  with  the  freezing  condition  gives  a 
closed  problem  that  yields  the  final  salinity  and  temperature  which 
determines  the  final  density.  Comparison  of  the  initial  density  to  the 
final  density  then  determines  the  buoyancy.  Conservation  of  heat  gives 

L'  * <T.  - VV-  022) 

where  l1  = L/Cp^^,  and  the  subscripts  f,  ® refer  to  final  and  initial 


149 


conditions.  Conservation  of  salt  gives 


Conservation  of  mass  gives 


(123) 


V 


f 


1 + V 

00 


(124) 


The  linear  freezing  condition  is  = -mS^  where  m is  a constant 
equal  to  0.057.  Solving  for  gives  the  quadratic 


S2 

5f 


(T  + L1) 


L'S 


= 0 


(125) 


The  solution  to  equation  (125)  is  the  dashed  line  on  Figure  26  which 
does  not  fit  the  data.  This  discrepancy  results  from  the  assumption 
that  all  of  the  cooling  and  diluting  takes  place  in  one  parcel  of  water 
and  the  effect  of  the  large  Lewis  number  is  neglected.  In  this  case, 
the  large  Lewis  number  allows  one  parcel  of  water  to  be  diluted  and 
another  parcel  of  water  to  be  cooled.  The  relative  sizes  of  these 
parcels  becomes  important  in  determining  the  net  buoyancy  across  the 
boundary  layer. 


Conclusions 


In  summary  these  experiments  show  that  the  flow  generated  by  a 
melting  vertical  ice  wall  immersed  in  salt  water  of  oceanic  temperatures 
and  salinities  consists  of  a laminar  bidirectional  boundary  layer  flow 
adjacent  to  the  bottom  portion  of  the  ice  which  becomes  a turbulent 
upward  flow  further  up  the  ice.  The  bidirectional  laminar  flow  results 
from  the  high  Lewis  number,  Le  = 200,  which  separates  the  cooling  and 
dilution  effects  within  the  boundary  layer.  The  low  value  of  0, 

'10~  -1  m ' s'1,  contains  the  saline  anomaly  within  3 mm  of  the  ice  where 
it  generates  an  upward  flowing  boundary  layer.  The  comparatively  high 
value  of  k.-,  10"  ’ m2  s'1,  allows  the  thermal  anomaly  to  diffuse  past  the 

saline  anomaly  where  it  generates  a downward  flowing  boundary  layer 
10  mm-20  mm  thick. 

The  large  effect  of  salinity  on  density  allows  the  inner  saline 
layer  to  determine  the  length  of  the  laminar  regime.  The  inner  layer 
becomes  turbulent  when  the  Grashof  number  reached  a critical  value  of 
108  which,  for  oceanic  cases,  corresponds  to  a maximum  laminar  length 
scale  of  0.5  m. 

When  the  inner  layer  becomes  turbulent,  the  low  salinity  water  near 
the  ice  diffuses  further  out  where  it  overwhelms  the  thermal  effect  on 
density  and  results  in  a unidirectional  upward  flowing  turbulent  bound- 
ary layer.  The  divergence  between  the  downward  flowing  laminar  thermal 
layer  and  the  upward  flowing  turbulent  boundary  layer  produces  a hori- 
zontal jet  of  ambient  water  flowing  towards  the  ice.  The  jet  marks  the 
top  of  the  laminar  region  and  the  origin  of  the  turbulent  flow. 
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Measurements  of  the  ice-water  interface  temperature,  made  by 
thermistors  frozen  in  the  ice,  show  a vertically  varying  temperature 
in  the  laminar  region  and  a uniform  temperature  in  the  turbulent  region. 
In  the  laminar  region  the  interface  temperatures  are  coldest  at  the 
bottom  of  the  ice  and  slowly  increase  with  height  until  the  turbulent 
region  is  reached.  In  the  turbulent  region  the  interface  temperatures 
are  uniform  over  the  ice  face.  For  both  regions,  the  interface  tempera- 
ture approaches  0°C  as  the  thermal  driving  increases  past  10°C. 

Measurements  of  the  melt  rate  in  the  turbulent  region  by  a measur- 
ing transit  show  the  ablation  rate  increases  as  the  thermal  driving  to 
the  3/2  power  below  Td  = 10°C,  while  above  Td  = 10°C  the  ablation  rate 
increases  linearly  with  thermal  driving.  In  the  region  between  the 
laminar  and  turbulent  flow,  the  ice  ablates  the  fastest  because  of  the 
jet  of  ambient  water  flowing  towards  the  ice  at  this  level. 

To  analyze  the  turbulent  flow,  I solved  the  boundary  layer  equa- 
tions of  motion  for  a turbulent  flow  under  the  assumption  of  a self- 
similarity. I modeled  the  turbulent  transfer  of  salt,  heat,  and 
momentum  with  an  eddy  diffusivity  that  varies  across  the  boundary  layer 
in  the  following  fashion:  Near  the  ice  the  eddy  diffusivity  equals  zero 
then  increases  to  a constant  value. 

I numerically  solved  the  similarity  equations  and  used  the  smoothed 
laboratory  data  to  determine  the  turbulent  parameters.  The  numerical 
results  show  that  the  maximum  value  of  the  eddy  diffusivity  is  of  order 
10"6  m2  s-1.  The  turbulent  diffusion  and  this  low  value  for  the  eddy 
diffusivity  further  supports  the  similarity  assumption.  The  main  direct 
effect  of  the  turbulence  is  the  increased  diffusivity  of  salt  that. 
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through  the  large  effect  of  salt  on  density,  indirectly  affects  the 

distribution  of  momentum  within  the  boundary  layer. 

T he  numerical  results  also  show  that  the  blowing  at  the  interface 

and  the  opposing  buoyancy  forces  have  an  important  effect  on  the 

turbulent  flow.  The  blowing  reduces  S which  increases  the  buoyant 

w 

driving  of  the  boundary  layer  as  Td  increases.  The  downward  thermal 
buoyancy  remains  important  because  the  high  Le  allows  for  the  production 
of  a small  amount  of  positively  buoyant  water  and  a large  amount  of 
negatively  buoyant  water  within  the  boundary  layer.  In  addition,  the 
increasing  effect  of  temperature  on  density  at  higher  temperature 
combines  with  the  Lewis  number  effect  to  reduce  the  upward  velocities 
and  transports  and  eventually  produces  the  flow  reversal  at  a Td  of  25°C. 
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APPENDIX  I 

The  Equation  of  State 


In  this  study  I used  the  equation  of  state  given  by  Gebhart  and 
Mollendorf  (1977)  at  1 atmosphere  of  pressure  because  of  its  compactness 
and  ease  of  manipulation.  The  authors  fitted  the  experimental  density 
data  to 

(M) 


p(S,  T)  = Pm(S){ 1 - >(S) (T-Tm(S) (ql 


where  r(S,T)  represents  density,  Pm(S)  is  the  maximum  density  at  a given 


salinity,  y(S)  is  the  thermal  coefficient,  Tm(S)  is  the  temperature 


of  maximum  density,  and  q is  an  unknown  coefficient.  They 


took  both  a linear  and  a quadratic  for  the  dependence  of  p , v,  and  T 

m m 


on  salinity  and  retained  pressure  effects  in  each  case. 


I used  the  linear  version  of  p , > , and  T at  1 atmosphere  to  com- 

m m r 


pute  p(S,  T)  and  a(S,  T)  in  the  numerical  calculations.  The  linear 
forms  are 

( I -2a ) 


Pm(S)  = + a:S  , 


y(S)  = a 3 + a4S  , 


( 1 -2b ) 


and 


Tm(S)  = a5  + a&S 


(l-2c) 


In  this  case 


1 .894 


( I -2d ) 


and  Table  A-l  gives  the  values  of  the  coefficients  at  through  af,  for  S 


It 
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In  parts  per  thousand. 

This  formulation  gives  accurate  densities  and  allows  for  easy 
computation  of  | and  | ||-,  two  quantities  needed  in  the  numerical 
computations.  For  oceanic  conditions, the  density  dependence  on  salinity 
is  nearly  linear  so  that 


6 = 8.05  x 10-“  °/00-i  , 


oo 


(1-3) 


while  the  density  dependence  on  temperature  is  more  complicated. 
Differentiation  of  (A-l)  with  respect  to  T gives 


a(S.  T) 

= 1 
p 

ff  , ±y(S)q | T-T 

a < 0 and 

for 

T < Tm,  a > 0. 

TABLE  A-l 

1.00 

a2 

8.046  x 

10'“ 

a 3 

9.30  x 

lO"6 

a„ 

-2.646  x 

10‘8 

a$ 

4.007 

a* 

-2.124  x 

10'1 

Values  for 

the 

coefficients  of 

the 

q-'  . 


equations  for  p.  y,  and  T . 

m m 


(1-4) 
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APPENDIX  II 


Molecular  Diffusivities  of  Salt,  Heat,  and  Momentum 

In  seawater  the  molecular  diffusivities  of  salt,  heat,  and 
momentum  differ  by  three  orders-of-magnitude;  the  large  difference 
between  the  thermal  and  saline  diffusivities  produces  a bidirectional 
flow  in  the  laminar  case  and  retains  a strong  influence  in  the 
turbulent  case. 

Caldwell  (1974)  gives  the  dependence  of  the  saline  diffusivity,  D, 
on  temperature,  T,  as 

D = 10” 1 1 (62.5  + 3.63T)m-’  s'1  . (11-1) 

Neumann  and  Pierson  (1966)  show  that  the  thermal  diffusivity,  ►,  has 
only  a slight  variation  with  temperature  such  that  k can  be  taken  as 
constant : 

s = 1.44  x 10"7  m-’  s‘l  . (II-:’) 

The  ratio  of  ►.  to  D defines  the  Lewis  number,  Le, 

Le  = jj  . (H-3) 

For  seawater  Le  varies  from  245  to  170  as  the  temperature  changes  from 
-1°C  to  +6°C.  This  study  shows  that  the  large  Le  produces  the 
bidirectionality  of  the  laminar  flow  and  it  is  partially  responsible 
for  the  high  temperature  reversal  of  the  turbulent  flow. 

Neumann  and  Pierson  (1966)  values  of  the  kinematic  viscosity  of 
momentum,  v,  in  seawater  for  a temperature  range  of  0°C  to  20°C;  they 


I 


! 


show  that  the  viscosity  dependence  on  salinity  is  negligible.  In  this 
study  I took  v at  0°C, 

v = 1.8  x 10-6  m2  s"1  . ( 1 1-4) 

I assumed  a constant  viscosity  because  ice  is  found  in  cold  water 
and  the  convective  motions  occur  in  water  cooled  by  the  ice. 
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APPENDIX  III 
The  Field  Study 

During  June  1977,  in  a cooperative  study  with  the  International  Ice 
Patrol,  1 investigated  iceberg  melting  processes  at  an  iceberg  in  the 
Labrador  Sea  at  50°  49'N,  52°  21 'W  aboard  the  USCG  Cutter  Ei v.*yr.v>:. 
Plate  XIII  shows  the  iceberg;  it  had  a waterline  length  of  100  m,  a 
height  of  30  m and  a draft  of  100  m.  The  draft  was  determined  by  side- 
scanning sonar.  The  arrow  A in  Plate  XIII  indicates  one  of  many  old 
waterlines  on  the  iceberg;  these  features  are  a result  of  increased 
melting  due  to  the  action  of  surface  waves  on  the  side  of  the  iceberg. 

To  therightof  this  waterline  the  subsurface  portion  of  the  iceberg  has 
been  exposed  to  view  by  rolling  of  the  iceberg.  The  iceberg  probably 
rolled  to  this  position  in  the  previous  two  or  three  days  because  the 
present  waterline  has  not  been  deeply  eroded.  Comparison  of  the  ice 
above  and  below  the  old  waterline  indicates  that  the  ice  melts  faster  in 
the  water  than  in  the  air  and  the  features  that  were  underwater  are 
more  rounded  and  polished  than  the  aerial  features. 

The  water  surrounding  the  iceberg  was  well  mixed  and  turbulent  as  a 
result  of  surface  waves  breaking  against  the  iceberg.  This  effect  was 
especially  pronounced  if  the  waves  were  incident  on  one  of  the  many 
sloping  ice  beaches.  The  presence  of  a swell  with  a period  of  approxi- 
mately 5 s and  an  amplitude  of  1 m made  extended  observations  next  to 
the  iceberg  impossible.  I observed  numerous  seabirds  and  small  blocks 
of  ice  near  the  iceberg. 

Because  the  laboratory  and  theoretical  results  show  that  the  upward 


flowing  turbulent  boundary  layer  is  only  a few  meters  thick  for  this 
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iceberg,  direct  measurements  in  the  boundary  layer  are  impractical  in 
the  field.  Therefore,  I measured  the  salinity  and  temperature  of  the 
surrounding  water  to  determine  the  alteration  of  these  properties  that 
results  from  melting  ice.  The  field  program  consisted  of  one  synoptic 
surface  survey  of  the  temperature  and  salinity  around  the  iceberg  and 
two  STD  sections  taken  radially  out  from  the  iceberg. 

The  synoptic  survey  consisted  of  continuously  measuring  the  temper- 
ature and  salinity  at  a depth  of  2 m while  the  ship  steamed  in  a circle 
around  the  iceberg.  The  ship  made  two  circuits  around  the  iceberg,  one 
at  a range  of  100  m to  200  m and  the  other  at  a range  of  500  m to  600  m. 
The  wind  during  the  surveys  blew  from  a direction  of  196°T  with  a speed 
of  10  m s'1 . The  inner  survey  showed  that  the  temperature  dropped  from 
4.2°C  to  3.4°C  in  a 150  meter  wide  plume  centered  at  17°  to  the  right  of 
the  downwind  direction.  I found  no  comparable  salinity  plume  in  this 
survey.  The  outer  survey  did  not  show  any  changes  that  could  be  attri- 
buted to  the  presence  of  a melt  plume.  Although  this  melt  plume  did  not 
extend  from  the  iceberg  more  than  0.5  km,  Bajzak  (1974)  found  a melt 
plume  from  a 300  m long  iceberg  in  the  Labrador  Sea  extending  for  a 
distance  of  over  1 km. 

As  a result  of  the  Coriolis  force  the  melt  plume  is  not  found 
directly  downwind  of  the  iceberg.  Neuman  and  Pierson  (1966)  show  that 
the  surface  current  is  at  an  angle  of  45°  to  the  right  of  the  wind  and 
the  current  veers  to  the  right  with  increasing  depth.  Therefore,  the 
iceberg  moves  to  the  right  of  the  wind  at  an  angle  determined  by  the 
depth-integrated  water  drag  while  the  surface  waters  move  to  the  right 
at  an  angle  of  45°  to  give  a melt  plume  location  of  approximately  17° 
for  this  case. 
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The  STD  sections,  which  determined  the  spatial  distribution  of 
temperature  and  salinity  with  depth,  indicate  upwelllng  may  occur  as  a 
result  of  the  turbulent  natural  convection  next  to  a melting  iceberg. 

Figure  27 shows  the  positions  of  the  11  STD  casts  taken  around  the  ice- 
berg (the  location  of  the  surface  melt  plume),  and  the  direction  of  the 
wind,  which  had  a speed  of  5 m s-1.  The  first  section,  taken  on  10  June 
1977,  consists  of  casts  2 through  6;  cast  4 determined  the  undisturbed 
water  column  properties  3.3  km  upwind  of  the  iceberg.  The  second 
section,  taken  on  11  June  1977,  consists  of  casts  7 through  12;  cast  7 
determined  the  undisturbed  water  column  properties  3.3  km  upwind  of  the 
iceberg.  I took  casts  2,  6,  and  12  within  2 m of  the  iceberg  from  the 
bow  of  the  ship  using  a portable  Inter-Ocean  STD  that  I hand-lowered  to 
depths  of  20  m.  I used  the  shipboard  Bisset-Berman  CTD  for  the  rest  of 
the  casts.  At  each  cast,  I used  deep-sea  reversing  thermometers  and  a 
Nansen  bottle  to  obtain  calibration  values  of  temperature  and  salinity 
for  the  CTD;  I took  these  measurements  at  the  greatest  depth  of  each 
cast,  usually  200  m.  I then  calibrated  the  Inter-Ocean  STD  against  the 
shipboard  CTD  by  taking  simultaneous  readings  during  cast  4. 

i 

Figure  28 shows  the  temperature  field  for  section  1 to  a depth  of 
100  m,  the  level  of  the  bottom  of  the  iceberg.  Below  100  m I found  no 
significant  differences  in  either  the  salinity  or  the  temperature  for 
all  of  the  casts.  Casts  within  100  m of  the  iceberg  showed  water  colder 
by  1°C  than  the  water  3.3  km  upwind.  Casts  2 and  3 showed  isothermal 
water  of  3.0°C  to  a depth  of  15  m out  to  100  m downwind  of  the  iceberg. 

Casts  5 and  6 showed  an  unexpected  result:  The  water  next  to  the  iceberg 
is  generally  warmer  than  the  water  30  m away.  This  results  from  the 
advectlon  of  warm  surface  water  towards  the  iceberg  by  the  turbulent 
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Figure  27.  Location  of  the  STD  casts  and  the  observed  melt  plume 
See  text  for  further  explanation. 
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motions  generated  by  the  incident  surface  waves.  Casts  3 and  5 showed 
colder  water  between  40  m and  70  m that  may  also  be  the  result  of  ice 
melting. 

Figure  29a  shows  the  salinity  field  for  the  first  section.  The 
salinity  at  the  iceberg  is  lowered  by  0.5  °/00  at  the  surface  and  then 
increases  to  the  far-field  value  of  33.1  °/00  at  5 m.  Cast  6 shows  the 
salinity  decreases  by  0.3  °/0o  at  18  m;  this  corresponds  to  the  1.0°C 
decrease  in  temperature  at  this  depth.  These  changes  occurred  when  the 
handline  of  the  portable  STD  went  slack,  indicating  that  the  STD  was  in 
contact  with  an  underwater  prow  of  ice. 

Figures  29b  and  30  i llustrate  the  temperature  and  salinity  fields 
for  section  2.  Near  the  iceberg  the  surface  water  was  2.9°C  colder  than 
the  surface  water  at  3.3  km.  The  water  temperature  to  20  m deep  was 
2.3°C  colder  than  those  found  in  the  far-field.  At  30  m away  from  the 
iceberg  the  water  temperature  is  2.1°C  colder  than  the  far-field  water 
at  4 m depth;  this  difference  decreases  to  zero  at  the  40  m level. 

Figure  29b  shows  the  salinity  for  section  2 and  the  far-field  salinity 
distribution  for  this  section  was  identical  to  the  distribution  found 
in  section  1.  In  the  upper  20  m,  at  distances  of  2 m and  30  m,  the 
salinity  had  perturbations  of  0.6  °/00. 

The  two  STD  sections  suggest  that  upwelling  occurs  as  a result  of 
subsurface  ice  melting.  The  temperature  gradient  in  the  upper  75  m 
makes  the  advection  of  cold  water  from  depth  visible  in  the  upper  layers. 
In  a 100  m region  surrounding  the  iceberg,  anomalously  cold  water  is 
found  down  to  40  m.  This  colder  water  results  from  the  combination  of 
cooling  due  to  melting  and  advection  of  cold  water  from  below,  through 
the  action  of  the  convective  boundary  layer  next  to  the  ice.  As  the 


6 CO 


169 

salinity  is  uniformly  almost  distributed  with  depth,  upward  advectlon 
cannot  be  seen  In  the  salinity  sections;  therefore,  any  salinity  changes 
must  be  due  to  dilution  by  Ice  melting.  The  observation  of  salinity 
fluctuations  in  the  upper  20  m and  temperature  fluctuations  in  the  upper 
40  m suggests  that  the  cold  less  saline  water  at  the  surface  Is  a result 
of  melting  and  advectlon,  while  the  cooling  between  20  m and  40  m is 
solely  due  to  advectlon  of  cold  water  from  below. 


APPENDIX  IV 


This  appendix  contains  the  tables  of  f,  <j> * , <j>,  e'  and  e 

generated  by  the  numerical  solution  of  the  governing  equations  with  the 
turbulent  parameters  determined  by  (114).  Tables  IV-1  through  IV-6 
contain  the  solutions  for  Td  from  2°C  to  12°C  in  2°C  increments 
respectively.  These  solution  retain  both  blowing  at  the  interface  and 
opposing  buoyancy  forces. 

Tables  IV-7  through  I V- 1 2 contain  the  solutions  for  Td  from  2°C  to 
12°C  in  2°C  increments  respectively.  These  solutions  retain  blowing  at 
the  interface  but  density  is  only  a function  of  salinity. 

Tables  I V- 1 3 , 14  and  15  contain  the  solutions  for  Td  = 2,  4,  and 
8°C  with  opposing  buoyancy  forces  but  no  blowing  at  the  interface. 
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A laboratory  study  of  melting  vertical  ice  sheets  immersed  in  sodium 
chloride  solutions  of  oceanic  temperatures  and  salinities  found  that  the 
convective  motions  generated  by  the  melting  consist  of  two  regions  confined 
to  a boundary  layer  adjacent  to  the  ice.  The  first  region,  at  the  bottom  of 
the  ice,  is  laminar  and  bidirectional,  with  a 2 to  3 rrn  thick  upward  flowing 
layer  next  to  the  ice  inside  of  a 10  inn  thick  outer  downward  flowing  layer. 

This  bidirectionality  results  from  the  large  difference  between  the  thermal  and 
, saline  diffusivitles  which  contain  dilute  water  near  the  ice  and  allows  the 
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20.  (cont.) 

- cooling  to  diffuse  further  from  the  ice.  Near  the  ice,  dilute  cold  water 
rises;  while  away  from  the  ice,  cold  saline  water  sinks.  Second,  further  up 
the  ice,  the  inner  layer  goes  through  a region  of  transition  and  becomes 
fully  turbulent  when  the  saline  Grashof  number  reaches  2 x 10s.  The  turbu- 
lent flow  is  upward,  unidirectional  and  the  dominant  flow  in  the  oceans. 

Ice-water  interface  temperature  measurements  over  the  ice  surface  show 
a vertically  varying  temperature  in  the  laminar  region  and  a uniform  tempera- 
ture in  the  turbulent  region.  The  measured  melt  rates  were  highest  in  the 
transition  region,  lowest  in  the  laminar  region,  and  had  a slow  vertical 
variation  in  the  turbulent  region.  For  water  temperatures  greater  than  25°C, 
the  turbulent  flow  reversed  direction  and  a laminar  bidirectional  region 
formed  at  the  top  of  the  ice. 

In  a theoretical  study  of  the  turbulent  flow,  an  eddy  diffusivity  models 
the  turbulent  transport  of  momentum,  salt  and  heat  in  the  conservation  equa- 
tions and  the  laboratory  measurements  determine  the  model  parameters.  The 
turbulent  model  consists  of  three  regions:  a laminar  inner  region,  an  inter- 
mediate region  where  the  eddy  diffusivity  increases  linearly  with  distance 
from  the  ice  and  an  outer  region  of  constant  eddy  diffusivity.  A similarity 
transformation  reduces  the  governing  partial  differential  equations  to  a 
system  of  ordinary  differential  equations  which  are  then  numerically  integrated. 

The  numerical  results  show  the  following:  The  eddy  diffusivity  depends 
parabol ical ly  on  the  density  difference  across  the  boundary  layer.  The  major 
effect  of  the  turbulence  is  to  increase  the  transport  of  dilute  water  away 
from  the  interface.  The  dilution  effect  overwhelms  the  thermal  effect  on 
density  to  produce  a unidirectional  upward  flow.  Finally,  the  injection  of 
melt  water  into  the  boundary  layer  reduces  the  salt  flux  to  the  ice  which 
lowers  the  interface  salinity  as  the  water  temperature  increases.  Use  of  the 
results  to  form  a zero  net  buoyancy  condition  on  the  buoyancy  integrated 
across  the  boundary  layer  predicts  the  reversal  of  the  turbulent  flow  which 
accurately  fits  all  of  the  known  data. 

Salinity  and  temperature  measurements  with  depth  in  the  vicinity  of  an 
iceberg  in  the  North  Atlantic  Ocean  show  variations  in  these  quantities  near 
the  iceberg  that  may  be  the  result  of  the  upward  flowing  turbulent  boundary 
layer. 


